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Abstract 


To investigate the influence of lightning strokes on surfaces with geologically non- 
homogenous structures, experiments were carried out in a model test tank filled with 
dried sand; and in this a vein of iron ore or other materials. Artificial lightning dis- 
charges against the surface of the test tank were undertaken by the aid of an impulse 
generator at voltages of one million volts and spark lengths of 75 and 170 cm respect- 
ively. It was found that a discharge from a positive electrode was pronouncedly in- 
fluenced by a geological discontinuity as an iron ore vein. On the other hand, negative 
strokes hit the sand surface at random more independent of geological discontinuities. 
The difference in effects from the two polarities is related to-the mechanism of the dis- 
charges and the development of leaders initiating the main discharges in the form of 


lightning ones. 


Several authors! have maintained that the 


geologic structure of the earth exerts an in- 


fluence on the path of lightning from cloud 
to earth and that therefore lightning will strike 
in some places more often than in others. 

One of the earliest statements in this re- 
gard was made by the Swedish author Olaus 
Magnus in his famous publication, »Historia 
de Gentibus Septentrionalibus”, Rome, 1555. In 
the first chapter of his sixth book, Olaus Mag- 
nus wrote: 

“Nature itself shows by many symptoms 
where these [ore deposits] are located. They 
will mainly be found on mountains which have 


1 See references at the end. 


rounded summits without ravines. When 
lightning strikes such mountains at the summit, 
on the side, or at the foot, it opens up the silver 
vein glittering in the clefts.” 

Following this same line of reasoning, pro- 
spectors have looked for ore deposits in areas 
where the local population reports that light- 
ning strikes frequently. This method has been 
used especially in certain regions of Africa. 

The danger of lightning is especially great 
over a non-homogeneous surface covered by 
soil that is a poor conductor. On the other 
hand, a surface that possesses the same electric 
qualities everywhere is not apt to produce 
concentrations of lightning, sometimes termed 


eels Ld , 
lightning nests”. 
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To confirm the above statements by obser- 
vations in nature is very laborious and largely 
a matter of accident. As we shall see, only the 
rarely occurring positive lightning is influence 
by geologic subsoil. This makes it even more 
difficult to obtain direct corroboration in 
nature. It is even doubtful whether direct ob- 
servations can yield ultimate proof. The geo- 
electrical qualities of the subsoil are variable, 
all circumstances cannot be defined with cer- 
tainty, and the influencing factors cannot 
always be identified. Therefore it seems desi- 
rable to test the validity of the hypothesis by 
model experiments in the laboratory. Attempts 
to obtain confirmation in nature should of 
course follow. 

Such experiments were made at the Institute 
of High Tension Research in Uppsala. This 
report is preliminary; a more exhaustive paper 
will be published later in the new Arkiv för 
geofysik of Kungl. Svenska Veten:kapsakade- 
mien, Stockholm. 


Experimental equipment 


The model consisted of a wooden box with a 
sheet iron bottom (2% 3m) which was set u 
horizontally and grounded. The box was filled 
with dry sand to an average depth of 22 cm. 
A streak or vein of dry iron ore or loam pow- 
der about 7 cm wide was arranged in the 
sand so as to extend across the box and to 
reach to the bottom. Granite stones were also 
placed here and there in the sand and small 
metal disks were put on top of it in an attempt 
to represent a geologically unequal subsoil. 

Vertically above the sand surface, at a 
stroke distance of 75 to 170 cm, a rod or ball 
(diameter 25 cm) was placed to act as a voltage 
electrode. The electrodes were placed some- 
times directly above the vein and sometimes 
displaced sideways. The impulse generator 
used had a nominal impulse wave of 1/50 usec. 
One pole of the generator and the bottom 
plate of the box were grounded. All experi- 
ments were conducted with minimal break- 
down voltage. 

There were made 100 discharges on the 
model surface for statistical evaluation of the 
paths. These paths were drawn on tracing 
paper. It was assumed that 100 discharges 
would suffice to exclude chance and to bring 
out any tendencies which might exist. 


Positive Discharges 


A number of preliminary experiments were 
made with discharges to a metal sheet covering 
the sand surface. When this was done, the 
sparks were distributed over the whole soil 
surface in accordance with the law of proba- 
bility. It became evident that both polarities 
have the samé properties. Thus we can con- 
clude that it is impossible to predict where 
lightning will strike on a surface that is a good 
conductor. All parts of the surfaces have the 
same electrical properties, and choices of light- 
ning path are quite accidental. 


When experiments were made on the sand 
surface containing the iron ore vein, results 
were completely different. In case of a geolo- 
gically unequal surface composed of matter of 
different electric properties the (conductivity 
of the iron ore vein is 300 times larger than 
that of the sand surrounding it), polarity has a 
decisive significance. The better conductors 
used in the model strongly influence the path 
of the lightning — but only in the case of posi- 
tive polarity. In such cases, “lightning nests” 
can even develop in a homogeneous sand 
layer due to small local inequalities of the elec- 
trical property. 

If it were possible to apply the results of the 
laboratory experiments to nature, we would 
observe the following properties of lightning. 
Positive lightning is influenced by the vein 
up to an angle of incidence of 45°. If the cloud 
is situated exactly above the vein, all the strokes 
will hit the vein, with preference for the 
edges. Even with an angle of 55°, veins will 
attract about 80% of the sparks. In this 
case, the so-called “creeping-paths”, form. 
The stroke at first hits the sand in front of the 
‘vein, then jumps to the vein, creeping along 
the surface. (Fig. 1.) 

This indicates that positive lightning prefers 
to strike superior conductors such as veins, 
faults, humid crevices, etc. Positive lightning 
therefore tends to be grounded with prefer- 
ence in such areas and sparks creeping on 
the surface will endanger objects located there. 
Thus, when there is positive lightning, the 
location of the cloud with respect to the vein 
determines whether the vein will protect or 
endanger objects. Not only veins, but also 
good surface conductors, such as moist cre- 
vices in mountains can have a moderate in- 
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Fig. 1. Positive discharge. The spark has hit at a distance 
of 13 cm in front of the vein and has then crept to 
the vein on the surface. Stroke distance, 75 cm. 


fluence on positive lightning. Strokes of light- 
ning have been observed which “ploughed” 
long paths along the surface, terminating in a 
crevice. 


Negative Discharges 


During the laboratory experiments, nega- 
tive lightning showed different properties than 
positive. Negative discharges when made 
to the homogeneous sand surface hit at ran- 
dom, as in the case of the metal plate. There 
is no tendency toward steering. This random 
distribution of strikes persists even with a 
non-homogeneous sand surface, for instance 
with vein. In correspondence with this re- 
sult one may conclude that if a negative cloud 
lies directly above the vein, one cannot expect 
that it will strike the vein as the positive light- 
ning does. Discharges will occur near the vein 
and sometimes the stroke will hit the vein. 
At other times, however, it will strike in the 


1*—904613 


surrounding sand surface. A tendency toward 
creeping paths was seldom observed in the 
model experiments with negative discharges. 

Thus, objects near the vein are not protected 
by its presence from negative strokes but are 
always endangered since the vein has some 
capacity for attraction, During lateral displace- 
ment of the electrodes, the vein loses its effect 
on the path of lightning more rapidly than with 
positive lightning. Under these circumstances, 
therefore, the existence of a well-conducting 
vein is less dangerous to objects during nega- 
tive than during positive lightning. The diffe- 
rent properties of both polarities are shown in 
tig 2) andes 


Fig. 2. The cloud is above the vein. A—E—arbitrary 
paths of lightning with negative polarity (left). No 
place, even at a distance from the vein, is safe from 
being hit. 
With positive polarity (right) the probable path of 
lightning lies between F and H. The other part of 
the surface is free from discharges. 


FA 


Fig. 3. The cloud is displaced laterally from the vein. 
A—D arbitrary paths of lightning with negative pol- 
arity. The whole surface is endangered. The vein in- 
creases the danger only by a small amount. 

Between E and F—preferred path;of lightning with 
positive polarity. The line e—f is especially endangered. 
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Positive discharge : 


Fig. 4. 
short lower (negative) leader. The place of juncture is 
given by the loop. 


Long upper (positive) and 


The different behavior of the positive and 
negative lightning is presumably due to differ- 
ent mechanisms of development. In order 
to answer this question, we have begun spe- 
cial examinations. Preliminary experiments 
have shown that upper and lower leaders are 
present where there is a poor conducting, non- 
homogeneous soil plate. Their meeting initia- 
tes the stroke. During positive discharges, the 
upper (positive) leader covers nearly the whole 
distance. The lower (negative) leaders are short. 
We suppose that “lightning nests” form where 
the ground can attain greater field strength 
because of stronger electric properties fav oring 
development of the lower (negative) leader. 
Fig. 4 shows a positive discharge with long 
upper and short lower leaders. The loop in 
the lightning trace indicates their meeting 
point. 

During negative leaders 


discharge, both 


Fig. 5. 

leaders (loop) lies about at the middle of the distance. 

Note especially the dying lower subleader at the right 
hand. 


Negative discharge: The meeting point of both 


have always been observed. Both have about 
equal length and junction occurs roughly in 
the middle of the distance (fig. 5). Lower 
leaders are found more frequently in this case. 
It is a question of chance which of the many 
lower leaders first joins the upper one to deter- 
mine the direction of the main discharge. In 
consequence, lightning hits all over the sur- 
face; the vein loses the character of concen- 
trator; even the portion of the sand surface 
that closely adjoins the vein is hit occasionally. 


Conclusion 


In summary, it appears that the model ex- 
periments have answered the question of the 
existence of “lightning nests” and shown that 
their formation is closely connected with the 
polarity of lightning. Positive lightning pre- 
ferably strikes places with good electric pro 
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perties. There is a marked tendency toward 
formation of “lightning nests”; and protec- 
tive measures are more effective. 

It is more difficult to predict where negative 
lightning will strike. It does not exclusively hit 
places with well-conducting subsoil; it also 
hits surrounding soil which is a poorer con- 
ductor. Protection therefore does not have the 
desired effect. The zone protected by an in- 
stallation is narrow; the formation of “light- 
ning nests’ is poorly developed. 


It is necessary to take note of these facts 
since natural lightning comes from clouds 
charged either positively or negatively. Per- 
centage distributions have often shown that 
about 90 % of all lightning comes from nega- 
tively charged thunderstorm clouds. There 
are, however, exceptions, depending on the 
conditions in the atmosphere. From this 
viewpoint, the question of the electrical polar- 
ity of the thundercloud has the greatest signi- 
ficance. 
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Principal Results of a Comprehensive Investigation of the 


Structure and Dynamics of the Thunderstorm 


By HORACE R. BYERS, University of Chicago 


(Manuscript received 26 August, 1949) 


Abstract 


Certain intensive investigations into the nature of thunderstorms in the United States 
have resulted in a new model of thunderstorm structure, life cycle and thermodyna- 
mics. Thunderstorms and their effects have been investigated by airplanes, radar, surface 
micronetworks and swarms of balloons. Some of the more significant aspects of these 
measurements are described. A comprehensive final report will be published soon by the 


U. S. Government. 


Introduction 


In 1945 the United States Weather Bureau, 
Air Force, Navy and National Advisory 
Committee for Aeronautics organized a pro- 
ject to investigate the structure and dynamics 
of the induvidual thunderstorm by means of 
a concentrated program of observations and 
measurements. A description of the observa- 
tion system has beenpublished in various ar- 
ticles (WHITE, 1946, 1948; U. S. WEATHER 
BUREAU, 1946, 1947). It will suffice here to 
state that it has been the largest project ever 
directed at the study of a single atmospheric 
phenomenon. Airplanes, balloons, radar, sur- 
face micronetworks and certain specialized 
equipment were used to investigate the 
structure and dynamics of thunderstorms 
during two summers of pronounced thun- 
derstorm activity. In the summer of 1946 
observations were carried’ out over the 
Florida Peninsula, and in 1947 south-western 
Ohio was the scene of operations. 

Analysis of the data has been completed to 
the point of making a general report of the 
scientific results. This report, prepared with 
the collaboration of Roscor R. BRrAHAM, 
JR, Louis J. BATTAN, Harry L. HAMILTON, 
JR, Harry M. Moses and others of our 


Weather Bureau group, is expected to appear 
in printed form early in 1950. It is being 
printed by the U.S. Government Printing 
Office under the direction of the U. S. Weather 
Bureau, which organization had charge of 
the scientific investigations. The report, un- 
der the title “The Thunderstorm” will be 
available from the Superintendent of Docu- 
ments, Washington, D.C. 

At the suggestion of the Editor of Terrus, 
a summary of the findings of the Project is 
presented on the following pages in order to 
acquaint meteorologists and geophysicists in 
general with the more important new facts 
about thunderstorms without waiting for the 
full report. The difficulty of this task can be 
measured by considering that the printed 
report will consist of some 300 pages of a 
size somewhat larger than this journal and 
including 254 illustrations. This summary, 
therefore, can only state conclusions, so the 
reader is referred to the full report to ascertain 
the nature of the data and reasoning leading 
to these conclusions. The results will be 
summarized here in conformance with and 
in the order of the ten chapters that form the 
report, with the exception that the last two, 
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covering thunderstorm conditions affecting 
flight operations, will be omitted. The preli- 
minary results of these two have already 
been published by the U. S. WEATHER BUREAU 
(1948) and the Aır WEATHER SERVICE (1949). 


1. Thunderstorm Structure and Circulation 


The BÉNARD cell (BENARD 1901) is rejected 
as a model of convective circulation in the 
cumulonimbus cloud partly because steady- 
state convection cells rarely, if ever, are 
duplicated in deep atmospheric layers. The 
atmosphere with strong convection is in a 
very disturbed state. Radar scans over a large 
area having thunderstorms show either a 
random distribution of the storms or a con- 
centration in squall lines. Downflow spreads 
over areas many times greater than the areas 
of upflow. The large areas of return settling are 
cloudless or have only growth-arrested cumuli. 
The energy released by condensation exceeds 
the energy of the simple BÉNARD-cell type 
convective motion. 

As indicated by previous investigators, it 
is found that the thunderstorm consists of 
several units of convection, called cells. Many 
of these units are initially detectable as separate 
echoes appearing on a radar ’scope, but even 
after they have formed a conglomerate 
thundercloud mass they can be distinguished 
through a study of the data measured in air- 
planes flying through them. The patterns of 
vertical motions serve to identify the cells. 
Between the cells is found ‘‘connective tissue” 
of inactive, i.e., nonturbulent, cloud. 

As described in a previous article (Byers and 
BRAHAM, 1948), the life cycle of the thunder- 
storm cell is divided naturally into three stages 
depending upon the direction and magnitude 
of the predominating vertical flow. They are: 


1. Cumulus stage—characterized by updraft 


through the cell. 


2. Mature stage—characterized by the pre- 
sence of both updrafts and down- 
drafts, at least in the lower half of the cell. 


3. Dissipating stage—characterized by weak 
downdrafts prevailing throughout the 
cell. 


Cumulus Stage—Measurements made in cells 
in this stage show that throughout the entire 


cell there is an updraft. The strongest updrafts 
are found at the higher altitudes and toward the 
end of this stage. The surface air converges to 
feed the updraft, but air also flows into the 
cloud from the unsaturated environment. This 
mixing of surrounding air or “entraining” is 
made possible through the evaporation of 
some of the liquid water carried in the updraft. 
The entraining is characteristic of all stages of 
the thunderstorm. 


Within the developing cloud the tempera- 
tures are everywhere higher than in the 
environment at corresponding heights. The 
greatest differences between the in-cloud and 
environment temperatures are found in areas 
of greatest updraft speeds, at the upper levels, 
reaching a maximum toward the end of the 
stage. 

A distinguishing feature of the cumulus 
stage as contrasted with the other stages is 
the absence of precipitation falling out of the 
bottom of the cloud. The flight crews flying 
within the clouds reported rain or snow, 
particularly near the end of the stage, but these 
hydrometeors appeared to be suspended, since 
no precipitation was observed at that time at 
the ground. The greatest concentrations of 
hydrometeors—liquid, solid, or both—occur 
at the freezing level and somewhat above it. 
The depth of the transition zone in which the 
water passes from a liquid to a solid state 
varies, depending upon the speed of the 
updraft. 

Mature Stage—TIhe initial occurrence of 
rain at the ground has been chosen as the 
beginning of the mature stage. Apparently 
this happens when the size and concentration 
of the drops or ice particles becomes so great 
that they can no longer be supported by the 
existing updraft and begin to fall relative to 
the earth. The drag exerted by the precipita- 
tion upon the ascending air is an important 
factor in forming a downdraft and, as a matter 
of observed fact, the beginning of rain at 
the surface and the initial appearance of a 
downdraft are coincidental. The downdraft, — 
which at first is found only in the lowest 
15,000 feet, gradually increases in horizontal 
and vertical extent. (Fig. 1.) 

Throughout the mature stage, an updraft 
still persists. Its magnitude is usually greatest 
at the highest levels flown (25,000 feet). It is 
also strongest in the early part of the mature 


8 , HORACE R. BYERS 


FEET 


| Vile et: 


* 


= RER EU 


ss 
NN 


abat ET eis Peat SE 
u lai ara oad 
tess ae oder nr 


ee 
7" 


era NS us pe 
Eh TU ya 
pe eee 
Reese 777, 7 ety née 
2 
RARE +280 


RAIN 


Horizontal Scole 
Draft Vector Scoe ID 3° trrsec Ice Crystals 


Fig. 1. Vertical cross section sketch of a thunderstorm 

cell in the mature stage with updrafts and downdrafts 

represented vectorially, typical hydrometeors and tem- 

perature distribution for summer thunderstorm condi- 
tions in the eastern United States. 


stage when speeds may locally exceed 100 
ft/sec. The downdraft is usually weaker than 
the updraft and has its greatest speeds in the 
lower part of the cloud. At the ground it 
is forced to spread out. Thus areas of down- 
draft, rain and lateral divergence coincide. 

Low temperatures are associated with the 
downdraft, with the greatest negative tempe- 
rature anomalies in the lower levels. The 
updraft is still warmer than the clear air 
surrounding the cloud and, of course, has 
considerably higher temperatures than the 
downdraft at corresponding heights. The 
warmest regions at each level are found in the 
area of maximum updraft speed, and the 
coldest areas occur where the downdraft is 
the strongest. 

As discussed in the 1948 Journal article, the 
updraft air follows an adiabatic reduction in 
temperature with ascent which is modified 


both by condensation and by mixing with 
entrained air, two processes which act in 
opposite senses. The result is that the updraft 
or “entraining-wet-adiabatic” rate is greater 
than the simple wet adiabatic. Therefore the 
updraft air is permitted to become only 
slightly warmer than the environment air. In 
the downdraft, evaporation is the principal 
process besides compression, and this appears 
to occur essentially reversibly and wet- 
adiabatically. At the low levels, then, the 
downdraft appears as a small cold-air mass. 
It forms the well-known pseudo-cold front 
as it advances along the surface of the earth. 

The greatest intensities in all aspects of the 
thunderstorm are reached in the mature stage. 
Heavy rain and strong winds are observed at 
the ground and severe turbulence, including 
short, intense gusts as well as updrafts and 
downdrafts are found in the clouds. It is during 
this stage that hail occurs, although this 
hydrometeor is not found in every storm. 
Above 26,000 feet no direct observations have 
been made of hydrometeors, but from radar 
echoes returned from cells in the mature stage, 
it is concluded that the size and number of 
particles are large to heights as great as 60,000 
feet in some cases. With extraordinarily 
strong updrafts (up to 100 ft/sec), it is possible 
for liquid water to be carried well above the 
freezing line. One case of heavy rain at 26,000 
feet, nearly 10,000 feet above the freezing 
level, was reported. 

Dissipating Stage.—The downdraft spreads 
over the entire area of the cell to start the 
dissipating stage. As the updraft is cut off, the 
Eka a of liquid water released becomes 
smaller, finally reducing the mass of water 
available to accelerate the descending air. 
Then complete dissipation occurs or only 
stratified clouds remain. The entire cell is 
colder than the environment as long as down- 
drafts and rain persist. With decreasing down- 
draft speeds the temperature within the cell 
rises to a value approximately equal to that 
of the environment at corresponding levels. 
The disturbance of the wind field at the ground 
also disappears. 


II. Thermodynamics of Entraining 


Reference is made to the work of STOMMEL 
(1946) and AUSTIN (1948) in showing the 
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thermodynamic results of entraining of en- 
vironment air in the updraft of a convective 
cloud. Computations of inflow rates from the 
swarms of balloons sent aloft around thunder- 
storms show that the cumulus cloud which 
develops into a thunderstorm entrains at a 
rate of approximately 100 per cent per 500 
mb of ascent, that is, it koche its mass as it 
rises through a pressure decrease of soo mb. 
This is considerably less than the entrainment 
rates of tradewind cumuli found by SToMMEL 
and perhaps is related to the lack of develop- 
ment of the latter. 

It is shown that when entraining takes 
place the properties of the updraft are consider- 
ably different from those to be expected from 
the conventional considerations of parcel 
ascent. In general, it may be said that an up- 
draft undergoing entraining has less liquid 
water and a lapse rate steeper than the moist 
adiabatic. The rate of temperature decrease 
with ascent and entraining depends on the 
temperature lapse rate and water-vapor content 
of the environment, as well as on the rate of 
entraining. The drier the environment or the 
more rapid the entraining, the steeper will 
be the updraft entraining-moist-adiabatic lapse 
rate. With high rates of entraining the updraft 
might theoretically be colder than the environ- 
ment. This was not found to be the case 
except in isolated or transitory conditions in 
the Thunderstorm Project observations. 

The conditions of free convection based on 
temperatures higher in the updraft than in the 
surroundings, place limitations on the rates of 
entraining that can exist. A series of charts, of 
which fig. 2 is an example, was constructed 
showing the difference in temperature between 
updraft and environment for a group of diffe- 
rent environments and different entraining 
rates. For each environment is shown the 
effective temperature difference between the 
updraft and the environment that would 
result from each of the several rates of entrain- 
ing. The effective temperature difference is 
the difference between updraft and environ- 
ment reduced by the amount which is requi- 
red to produce a buoyancy force equi- 
valent to the drag force exerted by the con- 
tained water. It is therefore the temperature 
difference available for overcoming other 
retarding forces and for accelerating the air. 
The examples show that entraining rates 
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Fig. 2. Plot of distribution with height of computed 


effective temperature difference between updraft air and 
environment air under the specified rates of entrain- 
ment. In this case the environment lapse rate is assu- 
med to be midway between the dry and saturated 
adiabatic rates with a convective condensation level at 
900 mb and 23°C. Relative humidity in per cent in 
the environment is considered to be 1/10 of the pres- 
sure in mb, i. e., 90 % at 900 mb., 80 % at 800 mb., etc. 


greater than 100 per cent per 250 mb are 
unfavorable for the growth of cumulus 
clouds to heights great enough to produce 
thunderstorms in American maritime tropical 
air. 

The lapse rate in a typical, well-developed, 
bulging cumulus cloud is represented by line 
A’ B’ in fig. 3, and the lapse rate of the environ- 
ment by line AB. If a saturated parcel is 
displaced downward from point C, it warms 
at the moist-adiabatic rate CD, if no environ- 
ment air is entrained. If environment air is 
entrained, the parcel warms at some other rate, 
such as CE. If the parcel is displaced downward 
beyond D or E, it will be colder than the en- 
vironment and sink. The mass of liquid water 
provides the means whereby a parcel in a 
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Fig. 3. Graphihal representation on a section of an 
adiabatic chart of the thermodynamics of the downdraft. 
AB is the environment lapse rate, while A’B’ is the 
updraft lapse rate. CD is a moist adiabat and CE repre- 
sents the downdraft lapse rate. If the downdraft reaches 
below E, it is colder than the environment. In U. S. 
thunderstorms, the pressure interval from C to E is 
usually less than 100 mb. 


thunderstorm can thus be forced below point 
D or point E, whence it continues as the 
thunderstorm downdraft. With sufficient liquid 
water available to maintain saturation, the 
parcel continues to the ground, reaching the 
surface with a temperature several degrees 
lower than the surface environment wet- 
bulb temperature. 


II. Thunderstorm Weather Near the Surface 


Rainfall—It was found that the rainfall 
pattern follows closely the arrangement of 
the cells and reflects to a considerable extent 
their stages of development. A single, isolated 
thunderstorm cell is comparatively rare in 
the regions studied and when it does occur, 
it is generally weak and not representative of 
the average thunderstorms. Usually a thunder- 
storm consists of a group of three or more 
cells adjacent to each other, each manifesting 
itself in the rainfall pattern. The first rain 
reaching the ground from a newly-developed 


cell is limited in area to a few square miles. 
Later, as the cell develops further, the rain 
area expands along with the increase in 
horizontal extent of the downdraft with 
which the rain is associated. The cold air of 
the downdraft spreads out after reaching 
the surface, while the rain falls directly to 
the ground, so an ever-expanding outer area 
of cold air without rain develops. In the 
dissipating stage of the cell, the rain area 
contracts while the cold air continues to 
spread. 

As would be expected, under a moving cell 
the area of rain advances approximately at the 
speed of movement of the cell itself. If the 
rainfall is considered with respect to the moving 
cell, it is found that the duration of moderate 
or heavy rain from a single cell may vary 
from a few minutes in the case of weak, short- 
lived cells, to almost an hour in the large, 
active ones. However, the duration of the 
rain period at a surface station depends upon 
such factors as the number, size and longevity 
of the cells passing over the station, the posi- 
tion of the station with respect to each passing 
cell and the rate of translation of the cell. 
Without separating the above factors, one finds 
that the average duration of rainfall at Florida 
as well as Ohio stations is about 25 minutes. 

The maximum rainfall rates occur at the 
location of a cell core within 2 or 3 min. 
after the first measurable rain from the cell 
reaches the ground and the rain usually remains 
very heavy for a period of 5 to 15 min. There- 
after the rainfall rate decreases, but much 
more slowly than it increased. Lesser rates 
occur around the edges of a cell. 

Wind Field—When a cell is in the early 
building stage there is a gentle inward turning 
of the surface wind, forming an area of weak 
horizontal convergence under the updraft. 
Except for changes in the electric field, this 
change in wind is the first surface manifesta- 
tion of the usual type of thunderstorm. As 
the cell grows and a downdraft develops, the 
surface winds take on an entirely different 
character, becoming strong and gusty as they 
flow outward from the downdraft region. 

In the mature stage, the outward-flowing 
cold air underruns the warmer air which it 
er, and a discontinuity in the wind 
and temperature fields is established. This 
discontinuity or discontinuity zone moves 
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outward with the spreading cold air of the 
downdraft. Strong horizontal divergence is 
measured in the surface winds of the down- 
draft ‘outflow. Values as high as 29 hr 
(8.06 - 1073 sec) have been evaluated from 
the surface micronetworks. 

In the slowest-moving storms, the outflow 
is almost radial, but in most cases the outflow 
field is asymmetrical, with wind speeds sub- 
stantially higher on the downwind side. This 
is due to the reinforcement or cancellation, 
as the case may be, of the radial flow by the 
prevailing air movement in the lower layers. 
Thus the wind discontinuity in most cases is 
easily detected only in the forward portions 
of the storm, where it appears as a micro-cold 
front. 

The cold dome of outflowing downdraft 
air has a form illustrated schematically in 
fig. 4. In this sketch the thunderstorm cell is 
considered to be in the mature stage and is 
moving from left to right. The cold air is 
represented as having spread out considerably 
farther on the downwind side of the cell 
than on the upwind side, as would be expected 
in a moving system. 

After the outflow has been spreading for 
15 to 20 minutes, the discontinuity zone will 
have traveled about s or 6 miles from the cell 
center. The surface winds near the discon- 
tinuity surface are still strong and gusty, but 
farther behind, inside the cold-air dome, the 
surface wind speeds have decreased so that the 
strongest winds are no longer underneath the 
cell itself. This increase in the wind speed as 
one approaches the discontinuity zone from 
the cold air side results from a continued 
settling of the outflow air, bringing high 
horizontal momentum downward. 

The discontinuity zone, which everywhere 
marks the limit of the cold downdraft air, 
is one of the most interesting of all surface 
weather features associated with the thunder- 
storm passage. With the passage of this zone 
at a station directly underneath the center 
of a cell in the early mature stage, there is a 
sharp increase in the wind speed, occasionally 
to De Neve force, and a marked reduction 
in temperature. The sharp increase in wind 
has been termed the “first gust”, since it 
often appears as the first major gust of a 
period of high, gusty winds. After the cold 
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Fig. 4. Schematic vertical cross section through the 

lower part of a thunderstorm cell in the mature stage 

showing the spreading cold dome and typical pattern 
of air flow. 


the boundary of this air decrease. It is found 
that, at least during some stages of the storm, 
the discontinuity zone at the ground moves 
faster than the normal component of the 
winds immediately behind it. This may be 
explained by the fact that there is rapid down- 
flow of colder air toward the ground along 
the boundary zone which thrusts it forward. 

In direction, the wind shows clockwise 
shifts in most cases. This was especially true 
in Ohio where the winds were usually south- 
west or south and shifted to west or north- 
west at the discontinuity. 

Temperature—The “first gust” and the 
“temperature break”, i.e., the point on the 
thermograph where the temperature suddenly 
starts its drop, are two of the most pronounced 
features observed at the surface, and they occur 
essentially together. At the time of formation 
of summer afternoon thunderstorms in the 
areas studied the temperature is usually above 
85°F, often in the 90’s. As a result of the rain 
and the downdraft, the temperature may 
reach a value as low as 65°F. without a change 
in air-mass. As the thunderstorm activity 
dies down after sunset, the temperature 
usually has recovered to an intermediate value 
representing a mixture of strongly cooled 
and less cooled or uncooled portions of the 
air-mass. 

The area affected by the cooling is many 
times greater than that over which rain falls, 
but the temperature change is, of course, most 
marked in the rain core (downdraft center). 
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Cooling may be detected as much as 15 to 
20 mi. downstream. Near the center of a 
mature cell the temperature reaches a minimum 
within 10 to 15 minutes after the temperature 
break; farther from the cell the temperature 
drop is much slower. The amount of the 
temperature decrease observed in any given 
storm varies inversely with the distance of 
the observation point from a cell core. The 
temperature discontinuity is sharp and well 
defined under the core but is less pronounced 
farther away. Since the downdraft is only 
a few miles in diameter, the area over which 
the first and most rapid temperature fall 
occurs is relatively small. As a result, strong 
horizontal temperature gradients are created 
after the downdraft first reaches the ground. 
Gradients exceeding 20°F per mile have been 
observed. As the storm ages, the cold air 
spreads out and the magnitude of the horizon- 
tal temperature gradient decreases. Regardless 
of the spread of the cold air, the area of mini- 
mum temperature remains in the general 
location where the cold downdraft made its 
first appearance at the surface, except in those 
cases where a new cell with its own downdraft 
and rain core develops over another part of 
the cold dome. 

Pressure.—Early in the cumulus stage a fall 
in surface pressure almost invariably occurs. 
This fall is observed before the radar echo 
forms, and it is recorded over an area several 
times the maximum horizontal extent of 
the echo. When the radar echo appears, the 
pressure trace levels off in the region directly 
underneath it, but continues to fall and fre- 
quently ata more rapid rate in the surrounding 
areas. The pressure drops in the cumulus 
stage are usually small in magnitude — less 
than 0.02 in. (0.67 mb) below the diurnal 
trend and take place over a period of 5 to 
t5 minutes. Following the fall, the pressure 
trace remains steady for as long as 30 minutes. 

The pressure falls appear to be caused by 
the combined effects of vertically accelerated 
air motions, the expansion of the air due to 
the release of the latent heat of condensation 
and the failure of the convergence near the 
surface to compensate fully the expansion or 
divergence aloft. Wind patterns in the vicinity 
of thunderstorms showed velocity convergence 
in the developing stages with divergence 
above 20,000 feet, suggesting a mass balance. 


In the mature stage, two features of the 
pressure trace—the “dome” and the ‘‘nose”— 
are recognized. The dome is registered at 
all stations to which the cold outflow air 
penetrates. The pressure nose, the abrupt, 
sensational rise that some meteorologists 
regard as typical of the thunderstorm, really 
occurs only at stations that happen to be 
passed by the main rain and downdraft just 
after they have first reached the earth in the 
beginning of the mature stage. It is superim- 
posed upon or may mark the start of the 
pressure dome. 

The displacement of the warmer air by 
the cold outflowing air from the downdraft 
results in the pressure rise initiating the pressure 
dome. A study of 206 thunderstorm pressure 
records from the surface micronetwork showed 
that in 182 of the traces there was a pressure 
rise associated with the arrival of the cold 
outflowing air. Since the rate and total amount 
of pressure rise depend on the slope of the 
cold air mass, the temperature difference 
between the cold air and the displaced warm 
air, the depth of the cold air itself and the 
speed with which the system travels, the 
most marked pressure changes are found near 
the cell core and decrease with distance from 
it. The areal extent of the pressure dome is 
similar to that covered by the cold air. There- 
fore the pressure remains high for a period of 
from one-half hour to several hours, depending 
on the amount of cold air involved. 

The distinction between the pressure nose 
and the pressure dome can be made only 
with difficulty on the conventional week-long 
barograph traces, but on the twelve-hour 
recording drums used on the Thunderstorm 
Project the distinction is clear. 

At any given time the pressure nose was 
usually detected by but one station per cell 
in Ohio, indicating that the maximum diame- 
ter of the area in which it occurred was less 
than 5 mi. It is thus apparent that it is caused 
by a transitory process which exists for only 
a brief period, namely the time of com- 
mencement of the mature stage when rain 
and downdraft reach the ground. Two 
effects which are thought to be important 
in creating the pressure nose are the weight 
of the suspended and falling water, and a lack 
of balance between convergence and di- 
vergence. Both of these effects are at a 
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maximum at this time. An accumulation of 
an average of 1 gm of liquid water per cubic 
meter from the cloud base to 35,000 feet, a 
conservative value, would increase the surface 
pressure by about ı mb, other factors being 
equal. When the downdraft becomes 
established, it reverses at least part of the ver- 
tical circulation pattern and reduces the 
divergence above 20,000 feet, or perhaps 
even reverses it to convergence. Should this 
happen before the divergence (outflow) in 
the surface layers becomes well established, 
a net increase in mass convergence could 
result -in a substantial surface pressure rise. 
For example, an uncompensated convergence 
of ı hr’ (a conservative value for thunder- 
storm convergence) within a layer only 1,000 
feet thick at 20,000 feet would increase the 
surface pressure at a rate exceeding 1 mb in 
3 min. Relationships of pressure changes to 
vertical motions in thunderstorms which have 
been a favorite subject for theoretical treat- 
ments by meteorologists could not be con- 
firmed by data of the Project. One of the 
difficulties is that a downdraft at the lower 
levels is under an updraft at a higher level, so 
that the vertical accelerations are opposed 
and tend toward compensating each other. 

After the brief pressure nose, the pressure 
remains at the value of the pressure dome 
which prevails for the particular thunder- 
storm. The dome persists through the dissipat- 
ing stage of the cell, after which the pressure 
returns to the trend prevailing before the 
passage of the storm. In the case of a thunder- 
storm associated with a cold front or a fast- 
moving squall line, the pressure remains high 
or even continues to rise as a result of cold-air 
advection or the passage of a wave in the 
pressure and wind field. 

Humidity.—An extraordinary phenomenon 
is noted in the hygrograph traces from the 
surface micronetworks in the form of a sharp 
decrease in the relative humidity in the midst 
of the heaviest downpour of rain. With the 
onset of rain, the relative humidity rapidly 
approaches but usually does not quite reach 
100 per cent, then, in many cases, drops 
suddenly to values as low as 60 or 70 per cent, 
returning to near saturation again in a few 
minutes. Such a fluctuation, termed the 
humidity dip, and illustrated in fig. 5, is 


associated with the rain area. In practically 
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Fip. 5. Example of humidity »dip» as shown on hygro- 

graph during heavy rain under thunderstorm of August 

14, 1947 near Wilmington, Ohio. The time divisions 

indicated on this hygrograph trace are for $-minute 
intervals. 


all cases the humidity dip occurred in a region 
of divergence in the surface wind, therefore 
in the outflowing downdraft. The temperature 
was usually decreasing, but if it had already 
reached its lowest point, as was sometimes 
the case, a rise in temperature of 2 or 3°F 
would sometimes accompany the humidity 
dip. 

The failure of surface air to become saturated 
during the heavy rain, together with the fre- 
quently observed relative-humidity dip, sug- 
gest that the downdraft air fails to maintain 
saturation as it descends, even in the presence 
of large concentrations of liquid water. Two 
processes are suggested: first, dessication of 
the air by condensation on cold precipitation 
particles; second, a time lag in the evapora- 
tion of the water drops so that there is in- 
sufficient accomodation to the increase in the 
saturation mixing ratio as the air descends 
to lower levels. 

The only measurements of raindrop tem- 
peratures were made by the Project at 3 feet 
above the ground. There it was found that in 
afternoon thunderstorms the rain temperature 
in the first few minutes of the rain period is 
several degrees lower than that of the ambient 
air. This would be sufficient to keep the rela- 
tive humidity from reaching 100 per cent. If 
the second process mentioned above is taking 
place, one would expect that the air of the 
thunderstorm downdraft would be heated at 
a rate between the moist- and dry-adiabatic 
and would reach the ground in an unsaturated 
state. 
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Development of New Cells—From a study 
of numerous cases of new cell development, 
the action of the cold outflowing air appears 
unmistakably as causing or contributing to 
the new growth. When two cells in the mature 
stage are within a few miles of each other, the 
cold outflows collide. The greatest frequency 
of new cell development was in the area 
between two existing cells whose edges were 
three or less miles apart. A three-mile band 
downwind was next in importance, then the 
lateral edges and least frequently the upwind 
or rear side. 

In some cases the time interval between the 
beginning of the outflow and the appearance 
of the new cell on the radar ’scope was too 
short to permit explanation of the new de- 
velopment as a result of the underrunning 
cold air or similar time-consuming process. 
There are cases in which one or a cluster of 
new cells came into existence, as indicated 
by the radar echoes, almost simultaneously 
with the initial or parent cell, which suggests 
that a preferred region of updrafts favors the 
development. 


IV. Electric Fields in Thunderstorms 


Owing to a lack of suitable equipment and 
qualified personnel, the electrical data of the 
Thunderstorm Project left much to be desired. 
However, by means of recording electric 
field meters mounted on top and bottom of 
the fuselage of the airplanes which operated 
satisfactorily on a few occasions, some data 
on the electrical configuration of thunder- 
storms were obtained and some crude ground 
measurements gave further information. 

The predominant charge distribution was 
that reported by previous investigators, namely 
a positive charge in the upper part of the cloud 
and a negative charge in the lower part, with 
the field reversed under the cloud to indicate 
the ground as positive relative to the cloud 
base. In five low-level airplane traverses under 
a thunderstorm in Ohio, interesting changes 
in the electrical field were measured. The 
first traverse showed an almost neutral field. 
Only eight minutes later, with increased cell 
height and heavy rain, a strong positive 
gradient (directed earthward) was recorded 
in the rain area and a strong negative gradient 
outside the rain area. Then, with dissipation 


of the cell and decreased rain intensity, the 
later traverses showed the positive field 
vanishing and being replaced by a negative 
one. The finding of the positive charge in the 
rain core on the second and third traverses 
appears to bear out the findings of G. C. 
SIMPSON and his co-workers (1927, 1941) 
that a positively-charged region may occur 
near the base of the thundercloud. The more 
commonly observed charge distribution with 
the cloud base negative with respect to the 
ground, was found under the updraft of the 
cell and under the entire cell after the time of 
maximum rain. The first result is in con- 
formance with well-known measurements by 
a number of investigators showing that rain- 
drops falling to the ground usually carry a net 
positive charge. 

In a study of conditions in 7 areas of marked 
vertical velocities in thunderstorms occurring 
on three different days in Ohio, negative 
gradients were found in all except the two 
traverses mentioned above. These negative 
gradients were characteristic of both updrafts 
and downdrafts. The two traverses showing 
positive gradients, noted above, recorded this 
sign only in the downdraft portion. 

From a study of point collectors located on 
the ground and a height-finding radar scanning 
the thunderclouds from a point about 20 miles 
away, the following conclusions concerning 
lightning were reached: 


1. The cloud tops (echo tops) reached a 
height where the environment tempera- 
ture was around —20° C before the first 
lightning occurred. 

. The maximum lightning frequency 
occurs at the same time as the cell reaches 
its maximum height. 

3. As the height of the cell top decreases, 

the frequency of lightning also decreases. 

4. It appears that a greater cell height (or 

lower temperature of the cloud top) is 
necessary to initiate lightning than is 
required to maintain it once it has started. 

5. The maximum frequency of lightning 

precedes the time of the maximum 
s-minute rainfall. 
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V. The Thunderstorm as Disclosed by Radar 


From the photographic records of the range- 
height-indicating (RHI) radar ’scope it was 
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possible to obtain data on the altitude of first 
formation of 66 radar echoes from convective 
clouds. A graph was made of the frequency 
distribution of the differences between alti- 
tudes of the tops of the initial echoes and the 
concurrent heights of the freezing level. 
The distribution showed a pronounced mode 
at 1,000 feet above the freezing level but the 
mean was at 42,200 feet. In convective clouds 
the echoes appear abruptly after the cloud 
has been visible to field observers for some 
time, suggesting a sudden release of large 
quantities of large water drops after penetra- 
tion above the freezing level, in accordance 
with the Bergeron theory. 

It was found that the rate of vertical growth 
of the top of a radar cloud echo agrees closely 
with the updraft velocity measured in that 
portion of the cloud if a correction factor of 
about 3 feet per second is added to the radar 
growth rate to account for the relative free 
fall of the attenuating snowflakes or other 
particles. The use of radar in this manner for 
measuring updrafts appeared to have a prac- 
tical application in detecting hail possibilities, 
since all observed cases of hail accompanied 
strong updrafts. 

It was noted that the thunderstorm tops 
ascended in a series of steps, appearing as the 
growth of new protuberances or »turrets» 
from the cloud top. Large clouds usually 
exhibited several turrets during their growth 
and, as a matter of fact, in 32 cases studied, 
the maximum height reached was correlated 
with the number of turrets formed by a 
coefficient of + 0.67 + 0.10. Each successive 
turret was higher than the preceding one and 
the mean lapse of time between successive 
turret peaks was 17.8 minutes. It is believed 
that each turret makes it easier for the following 
ones by increasing the moisture content in the 
cloud-top environment air which must be 
entrained in further growth. Less heat of 
condensation is robbed from the new turret 
by entraining than would be the case if it 
were standing alone in a dry environment. 
Each turret underwent a growth period 
followed by an interval of subsidence. The 
growth period averaged about 16 minutes in 
duration. The average vertical growth rate 
was about 18 feet per sec. and the subsidence 
rate about 12 feet per sec. 

The mean of the maximum heights reached 


by 199 Ohio storms as indicated from the 
radar was 37,500 feet. There was no significant 
difference in this respect between different 
types of thunderstorms, such as air-mass, 
squall-line or frontal. 

The boundaries of the radar cloud echoes 
were found to agree closely with the visible 
cloud limits except in the levels near the cloud 
base where non-echo-producing “outrider” 
clouds were common. Also the anvils and 
other layer-type lateral extensions usually 
did not return a radar echo. The horizontal 
and vertical extents of thunderstorms were 
found to be correlated, that is, the taller the 
cloud the broader it was. The greatest areal 
coverage was at an altitude around 10,000 
feet. The cross-sectional area was slightly less 
at the lower level and tapered at the top, 
forming a total cloud echo shaped like a 
rosebud. 

From radar scans covering an area of over 
55,000 square miles, it was found that on 
average thunderstorm days in Ohio 10 per 
cent of the area would be covered by cloud 
echo and on a day of maximum thunderstorm 
activity, 40 per cent of the area would be 
covered. It is found from indirect comparison 
that at 20,000 feet the in-cloud areas would 
average § per cent, and have a maximum of 
22 per cent, indicating the better chance of 
avoiding thunderstorms by flying high. (Other 
measurements showed that these high levels 
are the worst for flights within thunderstorms). 


VI. Effect of Environment Wind Field 


The causes and effects of vertical shear, 
including effects of entraining and momentum 
transfer in the drafts, are investigated. Some 
of the results have already been published by 
Byers and BATTAN (1949) and Markus (1949), 
so it is not necessary to summarize them here, 
although much new material is contained in 
the final report. 

With the aid of long-range radar and 
abundant upper-air wind data, the movements 
of thunderstorms are studied in relation to 
the wind fields in which they are embedded. 
A method was devised whereby the transla- 
tional component of the motion could be 
separated from the growth or dissipative 
component. 

It was found that the motion of the storm 
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corresponded most closely to the mean vector 
wind between the gradient level and 20,000 
feet. The correlations were better in Ohio 
than in Florida owing to a stronger, more 
consistent wind flow in the former region. 
In direction, the correlations were 0.95 or 
better in both regions. It was found that the 
speed of the cloud was less than that of the 
vector mean wind from gradient to 20,000 
feet. When the two speeds are plotted, a nearly 
straight line is formed, but the relationship is 
better represented as 


Or — 1.01 0.05) Wa.) 05020 


where Uy is the vector mean wind speed and 
U. is the cloud speed. 


VII. Preferred Areas of Thunderstrom 
Development 


Both areas studied were essentially flat, 
although the plain to the west and southeast 
of the Ohio network was sharply dissected. 
The radar made it possible to find out why 
thunderstorms formed where they did. 

In Florida the only topographical effects 
were the alternation of land, small lakes and 
swamps which are characteristic of the Central 
Florida landscape. By making proper allow- 
ance for displacement of a cumulus cloud 
downwind before it produced a radar echo, 
it was found that the number of new echoes 
per square mile of land was one-third more 
than the number per square mile of water. A 
qualitative estimate of the effect of swampy areas 
was made by plotting a map showing isolines 
of number of new echoes, again allowing for 
an appropriate displacement. Although there 
are some exceptions, the map shows a definite 
tendency for the preponderance of new echoes 
to be over the dry land. 

A study of the diurnal distribution of 
thunderstorms in the Florida area showed a 
distinct maximum of temperature and rate 
of thundercloud development from 1230 to 
1630, 75th-meridian time. The effect of the 
sea breezes from the two sides of the peninsula 
in causing the thunderstorms has already been 
covered in an article by Byers and RODEBUSH 
(1949). 

In Ohio, a study of the initial appearance of 
584 radar clouds on 21 days showed a larger 
number of echoes from day to day over the 


rougher parts of the terrain, especially to the 
northwest and southeast of the surface mi- 
cronetwork area. However, plenty of new 
echoes formed over the flat area of the network. 
Neither in Florida nor in Ohio was it 
possible to find a well-marked region from 
which the first thunderstorm of a series 
would develop. The thunderstorm hearth 
(Gewitterherde) so aptly described by von 
FICKER (1932) was not readily apparent. 


VIII. Squall Lines 


During the 127-day period from May 17 
to September 21, 1947 there were in Ohio 56 
thunderstorm days on which extensive radar 
data were obtained. On 32 of these days lines 
of thunderstorms were observed. The lines on 
6 days occurred along surface fronts, on 19 
days were ahead of surface cold fronts and on 
the remaining 7 days apparently had no con- 
nection with the surface fronts. As is usual 
throughout the eastern United States, the 
pre-cold-front squall line was the predominat- 
ing scene of thunderstorm activity. 

It was found that the pre-cold-front squall 
lines are not quite parallel to the cold front 
but have an orientation averaging 13 deg. 
clockwise from that of the front. Their 
orientation is somewhat counterclockwise 
from the orientation of the overlying 700-mb 
contour lines. The movement of the squall 
line, as distinct from the movement of the 
individual storms composing the line, is 
not very closely related to the speed of the 
following cold front; in many cases it moves 
faster than the front. The individual elements 
of the line, or storms, move with the prevailing 
upper winds as described in the preceding 
section. 

The pre-cold-front squall line on the synop- 
tic chart is frequently observed to last for 
periods as long as 24 hours, during which it 
may travel several hundred miles. When 
viewed on the radar, it was found that there 
is usually a zone of convective activity com- 
prising several lines and, although the zone 
may persist for a considerable period of time, 
the lines within it are constantly undergoing 
change—new lines forming and older ones 
dissipating. 

The elements in a line range from a single 
isolated echo 1 to 2 mi. in diameter to a large 
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aggregate of storms appearing on the "scope 
as a single, more or less homogeneous echo 
having a maximum dimension of over 30 mi. 
The number of separate elements in a line 
varied from 4 or § small echoes in the early 
stages to as many as 40 or 50 at the time 
of maximum echo intensity. For the most 
part, elements were found to form and dissipate 
on a given squall line, suggesting the existence 
of a preferred line formation. 

The squall line or, more properly, the squall 
line zone was found to be at an average distance 
of 170 mi. ahead of the cold front, with 
extreme values of 80 and 325 mi. Generally 
there was an area of relative inactivity between 
the cold front and the trailing edge of the 
squall zone; but on one or two occasions the 
convective activity finally extended back to 
the front itself. 

Generally the thunderstorms associated with 
the squall lines were little different from any 
others except for a tendency for them to be 
more severe, especially in producing effects at 


the surface, such as strong, gusty winds. 
Since the storms come in a line, there is a 
more or less connected discontinuity line from 
the cold outflowing downdrafts which is 
often mistaken for a true front. The return 
to tropical air-mass conditions afterward 
shows the local character of the discontinuity. 
The pressure “dome” of the combined storms 
often results in the creation of a long pressure 
ridge along the squall line. Sometimes a 
squall line of this character may after a time 
produce a true frontal discontinuity involving 
more cold air than can be accounted for by 
the downdrafts. 

Of the seven squall lines that were not 
associated with fronts, three were related to 
some type of large-scale cyclonic disturbance, 
two were on the west side of a warm anti- 
cyclone and the other two were in rather 
flat, uniform pressure fields. These lines were 
usually short compared with the prefrontal 
types and seemed to be less well defined. 
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Die aerologischen Gewittertendenzkarten 


Von A. SIMILÄ 


Schwedisches Institut für Meteorologie und Hydrologie, Stockholm 


(Manuscript received 22 June, 1949) 


Abstract 


A procedure is outlined for estimating numerically the factors which favor or oppose 
the development of thunderstorms. The resulting “thunderstorm tendencies” may be 
represented on synoptic charts; examples are given of the use of these charts in fore- 


casting. 


Das Problem der Gewitterprognose be- 
schaftigt sich mit dem Teil der Labilität, der 
zu den stärkeren Umlagerungsprozessen, be- 
gleitet mit elektrischen Entladungen, führt. 
Nach statistischen Arbeiten in Finland (Rossi 
1940) variieren die tiefsten Temperaturen in 
den Schauerwolken zwischen — 8° und — 15° 
C, aber in den Gewitterwolken zwischen — 15° 
und — 30°C. Das bedeutet, dass in den 
meisten Fällen in Nord- und Mitteleuropa die 
Schauerwolken nicht die soo-mb-Fläche er- 
reichen, während sich die Gewitterwolken 
über dieses Niveau erstrecken. 

Das Gewitterproblem umfasst hauptsäch- 
lich zwei Teile: 


1. Labilität unter dem Kondensationsniveau. 
2. Labilität in der Wolkenbildungsschicht 
(über dem Kondensationsniveau). 


Der erste Teil sucht eine Antwort auf die 
Frage: Wird eine Wolkenbildung durch 
Konvektion oder durch horizontal konver- 
gierende Luftströme stattfinden? Dieses Pro- 
blem, das in den meisten meteorologischen 
Lehrbüchern behandelt ist, wird hier nicht 
berührt. 

Der zweite Teil umfasst die Labilität in der 
freien Atmosphäre und soll klarlegen, ob die 
neugebildeten Konvektionswolken die Mög- 
lichkeit haben, so stark weiterzuwachsen, dass 
sie zur Gewitterbildung führen können. Die 


Labilität in der freien Atmosphäre kann man 
mit Hilfe der aerologischen Messungen be- 
stimmen. Das aerologische Stationsnetz ge- 
stattet jetzt auch, zweimal täglich eine synop- 
tische Karte über die atmosphärische Labilität 
zu machen, wenn man nur diese Aufgabe 
zweckmässig methodisch lösen kann. Das 
Ziel dieser Arbeit ist eine solche Methode 
anzugeben. : 

Unter Wolkenbildung durch Konvektion 
versteht man Wolkenbildung in einem ge- 
sättigten Luftstrom, der in einer ungesättigten 
Umgebung aufsteigt. Bei diesem Prozess 
kann man mit verschiedenen Bedingungen 
rechnen. Man kann voraussetzen, dass die 
Umgebung in Ruhe ist und kein Wärme- 
transport zwischen dem aufsteigenden Luft- 
quantum und der Umgebung stattfindet 
(NORMAND 1938). Eine andere Methode von 
J. BJERKNES (1938) und S. PETTERSEN (1939) 
setzt folgendes voraus: 

Durch eine horizontale Luftschicht steigt 
ein Luftstrom M’ mit einer Geschwindigkeit w’ 
auf und gleichzeitig sinkt ein Luftstrom M 
mit einer Geschwindigkeit w ab. Wenn diese 
Luftschicht so gross ist, dass die horizontale 
Luftströmung keine Ein- oder Ausströmung 
verursacht, müssen die aufsteigenden Luft- 
quanten gleich mit den absteigenden sein: 


M’w = — Mw w’ = positiv (1) 


w = negativ 


DIE AEROLOGISCHEN GEWITTERTENDENZKARTEN 19 


Weiter muss im Anfang ein barotroper 
Zustand herrschen und alle Bewegungen müs- 
sen adiabatisch sein. 

Die Temperaturdifferenz zwischen der auf- 
steigenden Wolkenluft und der absteigenden 
Luft in der Umgebung ergibt sich aus 


’ 


ar-12 y- rn] (2) 


ar=-12 |y—r+ nn] (3) 


y = gemessener Temperaturgradient 
I’ = trockenadiabatischer Temperaturgra- 
dient 


I” = feuchtadiabatischer Temperaturgradient 


Wir sehen also, dass ein feuchtadiabatischer 
Temperaturgradient nicht gross genug ist, um 
konvektive Wolken endlichen Querschnitts zu 
unterhalten. Es ist notwendig, eine grössere La- 
bilität zu haben, die gleichzeitig die grösstmög- 
liche Menge der konvektiven Wolken be- 
stimmt. Als maximale Menge konvektiver 
Wolken kann 5/10 angesehen werden, was 


Br 


und w = — w voraussetzt. 


Die Bedingung für indifferente Schichtung 
ist durch Gleichung (4) bestimmt: 


’ 


AT = 0, ARE nel 


ae 


(4) 


Man kann dann folgende Berechnung machen: 


M’ ; 
Pies ie gr (irons (6) 
PT NEST) (6) 
vr = N» = Mengekonvektiver Wolken. 


Bei labilem Zustand muss der Temperatur- 
gradient der Schichtungskurve 


eae AN AP (7 


In jedem Labilitätsfall ist die Labilität am 
grösstem im Anfangsstadium: N; =o. Die 
Labilität nimmt dann im Zusammenhang mit 
der Bildung konvektiver Wolken ab und 
erreicht das indifferente Stadium, wenn Glei- 
chung (6) erfüllt ist. 


Tk~Taso 


Graphische Bestimmung der AT-Skala für 
Gewitterinstabilitat. 


Fig. 1. 


Wenn man gemäss fig. 1. die Tem- 
peratur der Schichtungskurve im Konden- 
sationsniveau mit JT, bezeichnet, und die 
Temperatur, bei der eine Feuchtadiabate 
durch T, die soo-mb-Fläche schneidet mit 
Tr, und die Temperatur, bei der eine 
Trockenadiabate durch T; die soo-mb-Fläche 
schneidet mit Tp’, kan man approximativ 
die Gleichung (7) für den labilen Zustand fol- 
gendermassen schreiben: 


TR aan, Te TONE) 


Hier gibt die linke Seite die observierte 
Temperatur im soo-mb-Niveau. 


Wenn man eine Karte über die Gewitter- 
labilität über grösseren Gebieten machen will, 
kann man, um Arbeit zu sparen, statt Ty, 
welches die Berechnung des Konvektions- 
kondensationsniveaus für jeden Aufstieg er- 
fordern würde, die Temperatur in einer 
gewissen Druckfläche, die sich dem Konden- 
sationsniveau nahe anschliesst, nehmen. 

Nach Rossıs statistischer Arbeit ist in Finland 
bei konvektiven Wolken H; im Mittel = 1 270 

m. Da die Verhältnisse in Skandinavien 
ziemlich gleichartig sind, scheint es günstig, 
Teso mb Statt IT, zu wählen. Die Fehler, die 
dabei entstehen, sind nicht gross. Die Höhe 
der 8so-mb-Fläche bei einem Druck von 
1000 mb an der Erdoberfläche und einer 
Mitteltemperatur T,, von + 10° ist 1 400 gm. 
Im übrigen sind die Fehler folgende: 


| | 
N 
+5 


Le => 0 +5 +10 +15 +20 BEE 


Fig. 2. Diagramm für Gewitterinstabilitit Ar als 
Funktion der Temperatur der 8so-mb-Fläche (T0) 
und der soo-mb-Fläche (T;oo)- 


CN te Berechnete Instabilität 
Ho His << Tr zu gross 
Hoss < He Yes ZT, "zaiklen 
Hl) Berechnete Instabilität 
Hosp HB; Too Tu (ze klein 
Heso < Hr Tss0 > Tr zu gross 
EN Berechnete Instabilität 
He50 > Hi T 50 = Ib; richtig 
He50 = JE T 350 = iis richtig 


Da in den meisten Fällen in der Nähe der 
8so-mb-Fliche y ~ J” ist, kann diese Verein- 
fachung eingeführt werden, ohne dass grössere 
Fehler entstehen. 

Wenn man nun mit der Instabilität zwischen 
den 8so-mb- und soo-mb-Flächen rechnet, 
so erhält man, wie schon erwähnt, in Europa 
die Fälle, wo die Labilität so stark ist, dass 
Gewitterbildung wahrscheinlich ist. Wenn 
man statt dessen die 8so-mb- und 600-mb- 
Flächen als Grenzen für die Instabilitätsberech- 
nung wählen würde, würde man in den Re- 
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sultaten auch die labilen Fälle erhalten, die 
Regenschauer aber nicht Gewitter verur- 
sachen. Es gilt also hier, geeignete Grenzen zu 
wählen, um in den Endresultaten nur die 
Gewitterfälle zu erhalten. In anderen klima- 
tischen Verhältnissen muss man vermutlich 
andere Grenzen wählen. 

Man bestimmt also die Instabilität, die Ge- 
witter hervorrufen kann, auf folgende Weise: 


500 500 
Tr RES I 
Dy a 500 500 x 10 
IR Zu Tr 


Ar, das die Instabilitat der Temperaturschich- 
tung zwischen den 850-mb- und soo-mb-Flä- 
chen angibt, erhält man aus einer Skala, die 
bei den verschiedenen in Frage kommenden 
Temperaturen verwendbar ist. Wenn Ar posi- 
tiv ist, gibt es einfach an, wieviele Zehntel 
konvektive Wolken sich entwickeln können: 


Ar ZEN OC 
DATE Fo mit der Temperatur variiert, 
muss man die Ar-Skala für den Teil der Tem- 
peraturskala berechnen, der bei Gewitterlagen 
in Frage kommt. 

Um eine Vergleichsskala für Stabilität zu 
zu erhalten, erstreckt man dieselbe Ar-Skala 
auch über die Fälle, wo Tp? — À Be aot 
Fig. 1 zeigt, wie die Ar-Skala für jedes Ts50 
bestimmt wird. Ar =f (T0, 760) ist “als 
Diagramm dargestellt (siehe fig. 2). 


Die Einwirkung der Feuchtigkeit auf die 
Instabilität 


Sowohl in der Luftpaket- wie in der »Slice- 
Methode» wird nur die Temperaturstabilität 
ohne Rücksicht auf die Feuchtigkeitsverhält- 
nisse behandelt. Man kann zwar virtuelle 
Temperatur statt gewöhnlicher Temperatur 
anwenden, aber eine solche Korrektion ist 
ziemlich klein und auch nicht berechtigt, 
wenn man Rücksicht auf die Gewichtszu- 
nahme nimmt, die die Wassertropfen in einer 
konvektiven Wolke verursachen. 

Eine hohe relative Feuchtigkeit kann eine 
mitgeführte Eigenschaft einer Luftmasse, z. B. 
einer maritimen Luftmasse, sein, aber sie 
kann auch das Resultat einer allgemeinen auf- 
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steigenden Bewegung sein. Ebenso kann eine 
geringe relative Feuchtigkeit eine mitge- 
führte Eigenschaft sein (kontinentale Luft- 
masse) oder das Resultat eines Herabsinkens 
der Luftmasse. Es scheint mir auch, dass die 
Temperaturinstabilität im Durchschnitt zu 
klein ist im Vergleich mit der horizontalen 
Ausdehnung N; der konvektiven Wolken. 
Hier kommt zwar in den synoptischen Beo- 
bachtungen häufig ein Fehler vor, der durch 
die Perspektive und die Höhenausdehnung 
der konvektiven Wolken verursacht wird. 
Nach S. PETTERSEN 1940 werden s/ro kon- 
vektive Wolken, die ihre Basis in 1000 m 
Höhe und eine vertikale Ausdehnung von 
1 000 m haben, von der Erde aus als 8/10—9/10 
angegeben. Es ist jedoch wahrscheinlich, dass 
konvektive Vertikalbewegungen nach oben 
nicht vollständig gemäss der Gleichung M’ w’ = 
— — Mw kompensiert werden, so dass der 
Luftmassentransport durch eine horizontale 
Lamelle am Tage einen Überschuss der auf- 
steigenden Bewegung zeigt. Dieses wird teils 
am Abend und in der Nacht durch ein all- 
gemeines Herabsinken, teils durch Herabsinken 
in einem anderen, nahebelegenen Gebiet 
kompensiert. Wenn die Instabilität stark ist, 
muss die letztgenannte Kompensation die 
grösste Rolle spielen, und dann entstehen ab- 
wechselnd »Wellen» feuchter und trockener 
Luft. Natürlich spielt bei der Kompensation 
der aufsteigenden Luft, sowie bei der Feuchtig- 
keit der Luft die Orographie, besonders die 
Verteilung der Gewässer, eine gewisse lokal- 
bedingte Rolle, die in der Gewitterstatistik 
hervortritt. 

Gewitteruntersuchungen in den Vereinigten 
Staaten (ByERs und BRAHAM 1948) zeigen, dass 
bei Gewitterbildung Luft von der Umgebung 
mit der Wolkenluft durchmischt wird und 
der Feuchtigkeitsgehalt der Umgebung auch 
auf diese Weise auf die Gewitterkonvektion 
Einfluss hat. 

Um die relative Feuchtigkeit fiir Gewitter- 
prognosen verwenden zu können, muss man 
versuchen, die Grenzen zwischen »feuchten» 
und »trockenen» Gebieten festzulegen. Nach 
den aerologischen Messungen in Schweden 
scheint die Feuchtigkeitsgrenze zwischen Ge- 
witterlagen und Fällen ohne Gewitter un- 
gefähr bei einer relativen Feuchtigkeit von 
60 % zu liegen. Auf Grund von Rossıs statisti- 
scher Arbeit habe ich dann einen Mittelwert 


2— 904613 


für die vertikale Verteilung der relativen Feuch- 
tigkeit berechnet aus dem Ausdruck 


tr a Ui,5 + U, + U U, 
m 4 > 


wo Ui 5, U; u.s. w. die respektiven relativen 
Feuchtigkeitswerte in der Höhe von 1,5 km, 
3 km u.s. w. sind. Diese Mittelwerte wur- 
den mit dem Vorkommen von Gewitter ver- 
glichen. 

Daraus ergab sich folgendes Resultat: 


Gruppe ohne Regenschauer und Gewitter 
Um = 47 % 

Gruppe Gewitter > 15 % 
tungen U, = 67% 

Gruppe Gewitter < 15 % aller Beobach- 
tungen "U, = 50% 

Gruppe Regenschauer ohne Gewitter U, = 
= 58 % 


aller Beobach- 


Es scheint nahe zu liegen, Um = 60 % als 
neutralen fiir Wert Gewitterprognosen zu 
wahlen. Wenn man die Zunahme der Labilität 
durch die relative Feuchtigkeit mit 4, be- 
zeichnet, muss man versuchen, eine geeignete 
Skala für A, zu finden. Vorschlagsweise kön- 
nen wir 5 % relative Feuchtigkeit als Ein- 
heit in der A,-Skala nehmen. Dann erhält 
man folgende Tabelle. 


na Al, 

TOO. + 
95 a 
90 ale 
85 + 
80 == 
WS = 
70 ir 
65 Er 
60 
5» == 
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Fig. 3. Au 


als Funktion der Temperatur (T) und 
Taupunktstemperatur (Ta). 


Es ist jedoch wahrscheinlich, dass man den 
relativen Feuchtigkeitswerten des soo-mb- 
Niveaus grösseres Gewicht zumessen sollte als 
denen des 850-m- und 700-mb-Niveaus. 


A, ist in fig. 3. als Funktion der Tempera- 
tur (T) und der Taupunktstemperatur (Td) 
dargestellt. 


Gewittertendenzkarten 


Für jeden aerologischen Aufstieg werden 
Ar-Werte und die A,-Werte für 850, 700 
und soo mb graphisch bestimmt. Die Höhen- 
winde für dieselben Höhen werden auch 
eingetragen. Sie dienen zur Bestimmung der 
Bewegung der verschiedenen Gewitterten- 
denzgebilde und zeigen gleichzeitig die advek- 
tiven Anderungen des vertikalen Tempera- 
turgradientens. 


Bei der Analyse der Gewittertendenzkarten 
werden nur solche Gebiete, wo sowohl Ar wie 
A, positiv sind, als gewitterriskant betrachtet 
und mit Hilfe der Zahlenwerte dieser Fak- 
toren weiter untergeteilt. 


Gewittertendenzkarten zeigen nur die Labi- 
litätsverhältnisse in der Wolkenbildungsschicht. 
Die negative Gewittertendenz muss ja mit der 
negativen Gewitterprognose identisch sein. 
Aber die positive Gewittertendenz führt zur 
Gewitterbildung nur, wenn eine Labilität 
unter dem Kondensationsniveau existiert, die 


zur Wolkenbildung führt. 


Im Sommer, wenn die Einstrahlung sehr 
stark ist, sind tatsächlich die meisten positiven 
Gewittertendenzgebiete mit Eintreffen des Ge- 
witters verbunden. 


Gewittertendenzkarten vom 4. bis 7. Mai 
1949 bieten ein gutes Beispiel dafür, wie ein 
positives Gewittertendenzgebiet wandert und 
auch seine Intensität ändern kann. 

Am 4. Mai 1949 (fig. 4a) gibt es positive Ar- 
Werte über dem grössten Teil Westeuropas, 
aber die Luft ist so trocken, dass 4,-Werte 
allgemein negativ sind. Nur in West-Frank- 
reich gibt es feuchte Luft, die auf der Gewitter- 
tendenzkarte ein positives Tendenzgebiet ver- 
ursacht. Die Höhenwinde über West-Europa 
sind vorwiegend östliche, weiter im Osten 
südliche. Synoptische Wetterkarten weisen 
eine Gewitterbeobachtung in Nord-Frankreich 
und eine andere in Belgien auf. 

Am nächsten Tage (fig. 4 b) hat das positive 
Gewittertendenzgebiet sich nach Norden ver- 
breitet und die Intensität hat markant zu- 


_ genommen. In demselben Gebiet gibt es auf 


den synoptischen Karten beinahe hundert 
Gewitterbeobachtungen. Ein negatives Ge- 
wittertendenzgebiet mit grosser Gewitterstabi- 
lität nähert sich vom Westen dem europäi- 
schen Kontinent. 


Am 6. Mai (fig. 4c) befindet sich das posi- 
tive Tendenzgebiet weiter im Nordosten und 
ist etwas schwächer geworden. Die Anzahl 
der Gewitterbeobachtungen in demselben 
Raum steigt bis c:a 5o. Das stabile Gebiet vom 
Westen hat sich auch weiter nach Osten ver- 
breitet. 


Am nächsten Tage (fig. 4 d) gibt es einen 
schwachen Rest des positiven Tendenzgebietes 
mit einzelnen Gewittern über dem östlichen 
Teil des Kontinents. Im Übrigen herrscht 
eine grosse Gewitterstabilität über West- och 
Nordwest-Europa. 


Die Höhenwetterkarten zeigen, dass diese 
Gewitterinstabilität sich in einer degenerier- 
ten Kaltluftmasse, die über dem Kontinent 
in den untersten Schichten rasch aufgewärmt 
wird, entwickelt. 

Umgekehrt wird sie auf dem weiteren 
Wege nach Nordosten wieder schwächer teils 
wegen einer langsamen Erwärmung der Kalt- 
luft in der Höhe und teils wegen der Abküh- 
lung der untersten Schichten bei der Passage 
über die im Frühjahr kalte Ostsee. 
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Fig. 4 a—d zeigen die Gewittertendenzkarten vom 4. bis 7. Mai 1949. Das Stationsmodell ist folgendes: un- 
mittelbar neben dem Stationspunkt sind von unten nach oben die Au-Werte für 850, 700 und soo mb 
eingetragen, neben dem Au-Wert für 700 mb der AT-Wert. Soweit möglich sind auch die Winde für 
die drei Höhenlagen durch Pfeile angegeben und zwar punktiert für 850, gestreichelt für 700 und voll ge- 
zeichnet für soo mb. Die Ziffern an den Pfeilspitzen geben die Windgeschwindigkeit in Knoten an. 
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Seismic Sounding of Shallow Depths 


By H. HEDSTRÖM and R. KOLLERT, Stockholm! 


(Manuscript received 14 October 1949) 


Abstract 


After a short review of the development of “artificial seismology’’ a description is given 
of the modern field procedure when using the seismic refraction method for measuring of 
small depths, of the order of only a few meters. A small portable instrument equipment used 
for this purpose is described. The theoretical and physical basis for the seismic exploration 
method is briefly discussed, and by an example it is shown that laboratory tests on rock sam- 
ples can be used to calculate the transmission velocity of seismic waves. Some examples are 
given on the use of the seismic refraction method for civil engineering work. It is also rela- 
ted how the small portable instrument equipment previously described was used to measure 
the thickness of a glacier by means of the seismic reflection method. Repeated reflections 
(reverberations) were obtained and the smallest depth accurately determined in this way was 


96 m, 


The science of seismology, which dates back 
to the middle of the eighteenth century, devel- 
oped into a science in the modern sense during 
the last decades of last century, when instru- 
ments of sufficient sensitivity, damping, and 
so on, had been developed for registration of 
earth tremors. About the turn of the century 
a complete theory had been developed for 
the transmission of seismic waves through 
the earth. (WIECHERT and ZÖPPRITZ 1907, 
WIECHERT 1910.) 

From the studies of seismograms obtained 
during earthquakes the seismologists have 
generally assumed the existence of three main 
horizontal discontinuity surfaces in the earth’s 
interior, at depths of about 30, at 1,000 and 
2,900 km, respectively. (Attempts have also 
been made to subdivide the region between 
30 and 2900 km depth in 6 “concentric shells”, 
each with its characteristic transmission velo- 
city of the seismic waves.) The first disconti- 
nuity surface, known as the Mohorovitié 
Discontinuity, is considered to be the contact 
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surface between the “crystalline”” crust of the 
earth and the underlying “vitreous” silicate 
shell, which has a thickness of some 950 km. 
At the depth 1000 km one has the upper 
surface of a thick sulphide-oxide shell ch 
reaches down to a depth of 2,900 km, where 
it rests on the nickel-iron core of the earth. 

In a recent paper by L. Mintrop (1949) the 
author says: “It has been much easier for the 
seismologists to determine the core of the 
earth at a depth of about 2,900 km than any 
interface between two formations with differ- 
ent wave velocities within the crust of the 
earth.” However, in 1947 the depth to the 
“Peridotite layer”, with a wave velocity of 
8.1 km/sec was determined to be 27 + 3 km, 
from seismograms obtained at the great 
explosion of 4,000 tons of munitions on 
Heligoland (18.4.1947), in 19 stations with 
distances up to 608 km fromi the explosion. 
No deeper discontinuities were reported at 
the time. From a study of seismograms from 
this explosion and from a number of earth- 
quakes and also from previous observations 
by leading geophysicists, MINTROP now arrives 


| 


LR nds ll elk dé 


SEISMIC SOUNDING OF SHALLOW DEPTHS 


at the following “stratigraphy” for the upper- 
most part of the earth’s crust. (V-longitudi- 
nal = velocity of dilatational waves; V-trans- 
versal = velocity of distorsional waves.) 


V-long. V-transv. 


km km/sec km/sec Rock 
O— 4 3.5 2.0 Sediments 
4—13 s2 2.9 Granite 
13—28 6.5 3.6 Basalts 
Mohoroviéié Discontinuity 
28 — $7 8.1 4.5 Peridotite 
$7—I10 9.2 5.0 Pyroxenite? 
110—118 7.0 - (Plastic layer) 
118— 170 RO 6.2 2 
170—183 16.3 8.8 ? 
283-230 0232006 Ca. 5.0 Alternatingly 
plasticand sol- 


id formations 


Experimental study on “artificial” earth- 
quakes were first taken up about 100 years 
ago (MALLET 1851). Buried gunpowder charges 
were detonated, and the resulting elastic 
vibrations set up in ground were picked up at 
a distance by means of rather primitive “seis- 
moscopes”. It was observed that different rock 
formations were characterized by different 
velocities for the seismic waves, although 
accurate determinations of the velocities could 
not be obtained. H. L. ABBOT was the first 
observer to take advantage, in 1876, of the 
exploding of a large dynamite charge (25 tons) 
to measure the transmission velocities of seis- 
mic waves near the earth’s surface (ABBOT 
1878). Seismoscopes were set up at distances 
from 8 to 19 km. Unfortunately the results 
of the measurements were rather poor be- 
cause the sensitivity of the seismoscopes was 
too low. Some more experimental work was 
done up to about 1900, when it was first 
recognized that the seismometers would have 
to be much more sensitive than previously 
assumed in order to register the first weak 
pulse to arrive, before the arrival of later 
pulses of high amplitude. 

In the first decade of the present century it 
gradually became clear that there might be 
some practical and commercial use fot “aru- 


ficial” seismology. WIECHERT’s above men- 


tioned paper of 1910 marked a turning point. 
In an article the same year MINTROP made refer- 
ences to practical uses of the method, and in 
1914 J. Wırıp in an article wrote: 
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“Up to the present it has not been possible 
to quite accurately determine the hodograph 
(travelling time curve) of the seismic rays 
for shorter distances, and the topmost layers 
of the earth’s crust. However, it is the very 
characteristics of the vibrations of the top- 
most layers of the earth’s crust which could 
be of great practical importance in that they 
may give us information about those parts 
of the earth’s crust in which mining opera- 
tions are carried on” (Wie 1914). 


In 1919 Mintrop applied for a patent in 
Germany covering the use of “refraction pro- 
filing” for determining the depth and type 
of subsurface formations. Exploration work 
with this “refraction method” was first carried 
out in North Sweden by Mrntrop and his 
associates in 1922. In the following year the 
method was introduced for oil prospecting 
work in Mexico and 1924 in U.S.A., on the 
Gulf Coast, where it promptly led to the 
discovery of a number of oil fields during the 
next few years. In 1929 the number of “re- 
fraction seismographs” at work in U.S.A. had 
already risen to 175. 


The seismic “reflection method”, which is in 
principle similar to the marine echo sounding, 
has been worked out, from the beginning, 
in U.S.A., where a patent for the method was 
applied for in 1919 by J. C. KARCHER. At first 
comparatively little work was done with this 
method and technical progress was slow. 
About 1930, however, the method was defi- 
nitely established and by 1937 there were be- 
tween 225 and 250 “reflection crews’ working 
in U.S.A., as against 4 in 1929. 

Last year (1948) the number of seismic crews 
operating in U.S.A. was estimated at about 
450, corresponding to an annual expenditure 
for ‘artificial seismology” of the order of 100 
million dollars. 


In Sweden the refraction method was first 
used in Vasterbotten (ore prospecting purposes) 
in 1922, later on, 1938—40, in Scania for 
large-scale structure mapping of the Cretaceous 
sediments in that region, and from 1941 also 
for civil engineering purposes. For the seismic 
refraction work in Scania the distances between 
shot point and seismometers were of the order 
300 m to 17 km and the largest depths deter- 
mined generally under s—6 overlying layers 
with successively higher wave velocity were 
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Fig. 1. Places in Sweden where the seismic method has 
been used for determination of depth-to-bedrock. In the 
two years 1947 and 1948 alone §1 contract jobs have been 
completed, resulting in 4,381 depth determinations, each 
one a substitute for a bore hole down to solid bedrock. 


between 2,500—3,000 m, corresponding to the 
depth to the crystalline basement under the 
thick sedimentary cover in this part of the 
country. In this work, as well as in seismic 
reflection work carried out in the same area 
in the following years, the depths measured 
were thus still of the orders of kilometers or 
parts of kilometers. This is the case in all 
seismic exploration work in oil field terrain — 
the smallest depth that is said to have been 
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measured by means of the seismic reflection 
method, is about 200 m, and generally the 
depths measured are of the order of 1,000— 
2,000 m. 

In the above mentioned application of the 
refraction method for civil engineering pur- 
poses on the other hand, the depths involved 
are very often of the order of a few (2—3) 
meters only, and the depths will often change 
very quickly along a profile on the ground, 
necessitating a lot of careful detail work. 
Generally, the problem is to determine, ra- 
pidly, cheaply and accurately, the contour of 
a buried hard-rock surface under a cover of 
soil, sand and clay, glacial drift and hard mo- 
raine, which may vary in thickness from 1—8o 
m. In the last few years this work has reached 
considerable proportions, for Swedish con- 
ditions. The number of complete depth de- 
terminations of this work now totals about 
8,000 (average depth to bedrock 12.1 m) divid- 
ed on a number of contract jobs within Swe- 
den, mostly in connection with water power 


projects. (Fig. 1.) 


Field Procedure and Apparatus 


In the above mentioned paper of 1910 
WIECHERT described the principle of “refrac- 
tion profiling” in the following way: 


“The farther from the origin the earth- 
quake waves are observed at the earth’s 
surface, the deeper the rays have penetrated 
into the earth and the more of the assumed 
layers they have therefore traversed. If one 
therefore proceeds from the origin gradually 
in the distance, then initially only the first, 
later also the second, and then the third, 
and so on of the layers will come into action. 
One proceeds correspondingly with the 
computation. From the observations in the 
vicinity of the origin one draws conclusions 
first regarding the constants ... of the first 
layer ... The paths of all rays through this 
layer can then be calculated. From the paths 
of those rays which penetrate the first layer 
but remain close to it one draws conclusions 
regarding ... the second layer, and so on, 
in this manner.” 


This description fits very closely to the 
present practice of refraction profiling which 
is in civil engineering work carried out in the 
following manner, 
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Fig. 2. Seismic registration along a proposed tunnel line in central Sweden. 


A profile line is staked out on the ground 
and measured off from one end, where a shot 
point is planned for. Further shot points are 
placed in the points 30 m, 55 m, 70 m a.s.f. 
from the first shot point. The seismometers 
are then placed on the ground along the profile 
at intervals of 5 m. For best efficiency and 
good portability under average field conditions 
one uses generally 6 seismometers. The first 
seismometer set-up will therefore cover the 
section o—30 m of the profile, with 2.5 m 
distance from the shot points to the nearest 
seismometers. Each seismometer is connected 
to the recording equipment by means of a two 
conductor cable of 200 m length. The recording 
equipment consists of a set of 6 amplifiers, 4 


recorder with oscillographs, timing device 
and camera, and also equipment like shooting 
box, storage battery etc. From the recording 
apparatus two-conductor cables are layed out 
to the shot points at o and 30, but also to the 
shot points at 55, 80 and 105, and so on. After 
these preparations the shots (with $o—250 gr 
ne in the shot points 0, 30, 55 and 80 
are fired in rapid succession. Then the seismo- 
meters I—$ are moved into the section be- 
tween the shot points at 30 and 55 meters. 
Seismometer number 6, however, remains in 
its place and becomes number 1 in the new 
set up. The double record obtained in this 
station delivers the correction which has to 
be applied to the travel-times if the conditions 
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Fig. 3. Rope-way across a rapid, for the placing of 
explosive charges at the river bottom. 


in a shot point should change between two 
successive shots. (The preceding shot might 
have aerated the soil, driven out the pore 
water etc.) Then the shots in o, 30, 55, 80 and 
105 are shot successively. The seismometers 
are now moved into the third profile section, 
except the one at 52.5 m, which stays where 


it is. The firing line to shot point 30 is moved 
to shot point 130 and shots are fıred in 30, 55, 
80, 105 and 130. This procedure is continued 
throughout the entire profile length. 

For the planning of water power projects 
seismic depth-to-bedrock determinations often 
have to be carried out under lakes and rivers, 
and across rapids. 

When one is working across fairly quiet 
waters, a long (up to 1 km) wire, also used as 
measuring tape, is stretched out along the 
profile line from bank to bank or between 
anchored boats. It is supported by small floats 
which keep it near the surface. Along this wire 
the work is carried out from small boats. In 
this case it is preferable to have seismometers 
in fixed positions on the bottom at the points 
0, 30, 55 2.5.0. and to shoot in the points (at 
the lake bottom) corresponding to the seismo- 
meter positions of the above example. 

When profiling across rapids one makes use 
of “travelling shotpoint”, while the seismo- 
meters remain fixed at their stations on each 
bank of the river. Two wires are stretched 


Fig. 4. Modern seismic equipment for depth-to-bedrock determinations. The heaviest 

piece weighs 53 pounds (24 kg), the total weight is 330 pounds (150 kg.). — In front 

8 seismometers, at left behind seismometers the amplifier (for 6 ‘“‘channels’’), at right 

behind seismometers the recorder with oscillograph, timing device and camera, at extreme 

right the “fool-proofing’” device for firing of shots and simultaneous starting of film 
in the recorder. In the background an accumulator battery, cable reels etc. 
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Fig. 5. Seismogram from a shot with a charge of ı kg dynamite, recorded in stations which 
were 220 m to 320 m apart from the shot point. 


across the river at a distance of about two 
metres from each other. One wire serves to 
suspend the cables which connect the seismo- 
meters on the other side of the river with the 
recording equipment. On the other wire, 
which is actually an aerial rope-way, a device 
is mounted which can be drawn across the 
river by means of a light wire (fig. 3). The 
explosive charge (0.2—0.5 kg dynamite), 
weighted with a heavy piece of scrap iron or 
a stone, is fixed to this device and pulled out 
to a certain point along the profile. By a 
sudden pull in a “release” wire, a clamping 
device which holds the explosive charge lets 
go and the charge drops to the river bottom 
without appreciable off-drift. The charge is 
then detonated and the procedure is repeated 
in a point 5 m farther away along the traverse. 

The seismometers referred to above pick up 
the ground vibrations caused by the explosions 
and transform them into electrical pulses. A 
seismometer belonging to the equipment 
shown in fig. 4 consists in principle of a per- 
manent magnet, rigidly connected to the 
casing, and a moving armature which is sus- 
pended by springs. The vibrations cause 
changes of the air gap between magnet and 
armature, causing changes in the magnetic 
flux due to which a voltage is induced in the 
coil around the armature. This voltage is as 
high as s mV, measured over a load resistance 
of soo ohms for vibrations of 40 c.p.s. of 


0.001 mm amplitude. The seismometers are 
shaped like hemispheres with the center of 
gravity near the large flat bottom. Due to the 
“streamlined” shape the effect of wind noise 
is reduced to a minimum. The construction 
is waterproof and very solid, and clamping 
(arresting) of the seismometers is never 
required. 

The amplifier mentioned above and shown 
in fig. 4 has 6 “channels”, one for each working 
seismometer, in which the weak pulses from 
the seismometers are amplified (some 60,000 
times) before they are delivered to the re- 
corder. The amplification for each channel 
can easily be adjusted so that a suitable de- 
flection is obtained on the corresponding 
oscillograph galvanometer (as can be observed 
on a ground glass window on the recorder.) 
The filament voltage for the amplifier unit 
— 6 volts — is delivered from a storage 
battery, the plate voltage from dry batteries. 

The amplified pulses are photographically 
recorded on a fast moving photographic paper 
strip of 100 mm width (fig. 5). The recorder 
contains 7 galvanometers and all other means 
for photographic registration like light source 
and optical system, motor for paper transport 
and camera. It houses also a timing device of 
the shutter type (tuning fork with slotted 
diaphragms on the prongs) which projects 
time lines at intervals of 0.01 sec. on the 
record. One of the galvanometers is connected 
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Fig. 6. Travel-time graph, ray path and wave front 
diagram calculated for a three 3-layer assumption. 


to the shooting circuit and records the shot 
instant, the other 6 galvanometers register the 
vibrations picked up by the seismometers. 
The recording paper is transported by a D.C. 
motor, the speed of which is kept constant by 
a governor, which can be adjusted to a feed 
of so to 100 cm paper per second. Thus the 
time difference between shot instant and 
arrival of the impulses can be determined with 
an accuracy of better than one millisecond, 
depending upon the shape of the onset of the 
impulse. As a rule the travel time is read off 
with an accuracy of 0.2 milliseconds or better. 
The over-all amplification factor of the 
equipment, i.e. the ratio of the amplitude on 
the record to the amplitude of the ground 
vibrations is of the order of 4 x 10°. 

The sensitivity of the equipment is more 
than sufficient for refraction work on a large 
scale, where usually profiles of several kilo- 
meters’ length are investigated. (The ampli- 
fication will always have to be decreased until 
“the normal ground unrest” does not show 
up too much on the seismograms, and the size 
of the explosive charges used is then determined 


from this “permissible” amplification.) For 
work at which many short profiles have to be 
investigated, like at depth-to-bedrock deter- 
minations, several special features have been 
incorporated which facilitate the handling of 
the equipment and increase the efficiency. 
One of these arrangements is that the shot is 
automatically released by the camera at the 
instant when the recording paper has the 
proper speed and that recording is automati- 
cally stopped after a certain time. Another 
feature is the balanced shot circuit, resulting 
in transmission of the shot instant only, though 
the shot is fired via the same transmission line. 
Moreover the equipment is extremely portable. 
No part of it weighs more than 24 kg and the 
total weight, including cables and standard 
storage battery is about 150 kgs. The equip- 
ment can thus be transported by the normal 
working team generally consisting of 6 men, 
either on back or on simple transporting 
means like bicycle trailers or on sledges, 
depending on topography and season. 


Interpretation and calculation work 


If energy is suddenly imparted to the ground, 
e.g. by a blast of explosives, a small fraction 
of the energy is transmitted away from the 
shot point in the form of elastic waves of 
longitudinal and transversal types. Their ve- 
locity of propagation is widely different in 
different geological formations, and seismic 
exploration methods as described above can 
therefore be used to identify geological forma- 
tions hidden in the subsurface. Hereby in 
general use is made only of the longitudinal 
wave, the velocity of propagation of which in 
the encountered formations is found in the 
following way. 

The time elapsed between the shot-instant 
and the arrival of the first pulse at each seismo- 
meter station is read off from the seismogram. 
This time is plotted as ordinate to the corres- 
ponding distance of the seismometer station 
from the shot-point. When these “plots” are 
connected, a broken line, a so-called travel- 
time graph is obtained (fig. 6). The first part, 
beginning in the origin, corresponds to pulses 
from direct waves through the first (= surface) 
layer and therefore the velocity of propagation 
in this layer Vi = dx/dy. In the points on the 
next part of the travel-time graph, those waves 


SEISMIC SOUNDING OF SHALLOW DEPTHS ZT 


which have penetrated to larger depth and 
travelled with high speed along the surface of 
the second layer and back to the surface have 
overtaken the waves in the first layer and are 
the first ones to arrive. The ratio dx/dy for 
this part of the travel-time graph gives the 
true value for the velocity of wave propagation 
V: in the second layer, provided that the 
interface is parallel to the surface (and an 
apparent velocity when this is not the case). 
For the determination of the true value V2 
one will therefore generally need travel-time 
graphs shot from both ends of the profile. 
Some velocities of wave propagation which 
are usually encountered for soil types and 
rocks in Sweden are given in the following 


table. 


Soil types and rocks V in m/sec. 


eh eth ae A 
glacial drift (boulders embedded 

inberispsande idlay)* is. en 1 
hard Moraine. 5.2.20. vs 


300— 500 


I,200—1I,700 
2,200—3,000 
$,000—$,$00 


The travel-time curves yield, however, also 
data from which the depth to the interfaces 
of the different layers can be calculated. As 
the geometry of the “ray path”, corresponding 
to the shortest time of travel, is determined by 
the same principle which FERMAT developed 
for the propagation of light, the formulas for 
depth calculations are easily derived. For the 
simple case of 3 layers with interfaces parallel 
to the surface, as illustrated in fig. 6, the 
following relations give the travel-time for 
any station of the travel-time graph: 
Section I: tt =A/Vı 
Section 2: t2 = T2 + A/V2 = 2hr cos ir2/Va 


Section 3: t, = T; + 4/V,; =2hr cos ir3/Vr 

+ 2 hz cos i,,/V2+ A/V; a.s.o. 
where fr, f2 and f3 are the traveltimes for points 
in the three sections of the graph at any 
distance A from the shot point, T; and'T, 
the time intercept on the ordinate in the shot 
point, Vr, V2 and V, the true velocities in the 
corresponding layers, fr and h: the thickness 
of the first and second layer respectively and 
i the respective angles of total reflection. From 
these expressions, formulas for the thickness of 
the layers can be derived. Thus for the first 
layer we may write 


IES = 2 hy COS i12/Vs 
and therefore 


hy = Jada etc: 


2 COS er 


An other type of formulas for depth calcula- 
tions uses the distance x of the breaking points 
in the travel-time graph from the shot point. 
The corresponding expression of the depth hı 
to the first interface (= thickness of the first 


layer) is 

Re x ve — VW 

2 Ve 4- Vi 

The depth formulas become rather complicated 
when dealing with several layers, especially 
when the interfaces are not parallel to the 
surface nor to each other. For these cases each 
computor usually has diagrams of his own 
construction in order to simplify the calcula- 
tion operations. In many cases it will be still 
more convenient to evaluate the depths 
graphically by means of wave-front construc- 
tions according to HUYGENS principle. A 
description of this procedure will, however, 
carry too far for this paper and reference must 
be made to seismic literature. (ANSEL 1931, 
THORNBURGH 1930.) 

The velocity of propagation of the elastic 
waves in the ground is directly related to the 
elastic constants of the geological formations 
penetrated by the waves. It can be shown 
mathematically that this velocity, for the 
longitudinal and the transversal waves, is given 


by the formulas: 


a 
PRE ES 2 
DIN Bee 0 


and Virans = \/ £ 


where 

E = Young’s modulus (stress-strain ratio for 
elongation or compression) 

rigidity modulus (stress-strain ratio for 
change of shape e.g. twisting) 


G 


I 


| 


= compressibility (fractional change of 


volume per unit increase of hydrostatic 
pressure) 
oe = specific gravity (volume weight) 
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Fig. 7. Section showing bedrock profile and overburden layers as determined through a seismic investigation. 


Thus from the velocities measured in the 
field the elastic constants for the material of 
the “‘subsurface’” can be calculated and con- 
versely the velocities can be calculated when 
these constants are known. The elastic con- 
stants can be determined by laboratory tests 
on samples either by static or by dynamic 
methods. Static methods usually do not give 
accurate results on samples taken near the 
surface, (weathered, fractured), due to perma- 
nent “setting”, or hysteresis phenomena, after 
the first application of the load, while the 
dynamic methods furnish satisfactory and 
reproducable results in most cases. Through 
Ipe’s method (1935, 1936), Young’s modulus 
is determined from the natural frequency of 
longitudinal vibrations in the length direction 
of a bar cut out of the material to be investi- 
gated, and for determination of the rigidity 
modulus the bar is brought into torsional 
vibrations. A comparison of the elastic con- 
stants as obtained by this method and as cal- 
culated from the velocities found in seismic 
field work will generally show good agree- 
ment. 

As example may be cited an instance where 
drill cores from an exploration drill-hole in 
the southernmost part of Sweden were subject- 
ed to laboratory investigations. The method 
used for the tests was similar to IpE’s above 
mentioned method, i.e. the test samples were 
subjected to vibrations of varying frequency, 
and a “‘spectrum’’ of resonance frequences was 
determined. In all 47 drill cores were tested, 
from depths between 328 and 1,097 m. From 
the elastic constants determined, and from the 
volume weights, the “velocity of sound” in 
each sample was determined. From these 
figures an average velocity for the longitudinal 
seismic wave between 330 and 1,100 m depth 
was determined. Further, the corresponding 
velocity down to a depth of 330 m had pre- 


viously been determined through seismic work 
with the refraction method. The combination 
of these results gave an average velocity down 
to a depth of 1,100 m of 3,350 m/sec. Reflection 
work in the same area, with the reflection 
profiles extended comparatively far out from 
the shot point, had indicated a “vertical 
velocity” from a depth of about 1,250 m of 
3,250 m/sec. In this same area reflections had 
been obtained at 0.772 sec., at 1.092 sec. and 
at 1.307 sec. Using the average velocity of 
3,350 m/sec., the reflection time 0.772 sec. gave 
a depth to the reflecting surface of 1292 m. 
At a depth of 1,285 m the drill hole passed from 
a thick, sandy and porous formation into a 
clayey, harder formation of considerable 
thickness. No doubt this interface was the one 
that had caused the reflection at 0.772 sec. 
This goes to show that very good agreement 
can be obtained between field results and dy- 
namic test results on samples in the laboratory. 


Practical Results 


The field procedure described above has 
been developed on the base of experience from 
the interpretation of tens of thousands of 
seismograms. It was soon found that travel- 
time graphs, obtained from shots at the end 
points of a profile of no greater length than 
necessary in order to reach the depth required, 
do not furnish sufficient data for a correct 
interpretation. One reason for this is that the 
velocities that are characteristic of the upper- 
most layers of the overburden often vary 
irregularly and quickly along the profile, due 
to variations in weathering and aerating. Other 
reasons are irregularly varying thickness, along 
the profile, of the different layers in the over- 
burden, (fig. 7), and rough configuration of 
the bedrock surface (quickly changing depth 


along the profile). Therefore the investigation 
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must furnish sufficient data to give a continuous 
record of the changes of velocity in the over- 
burden layers and especially in the surface 
layer—where it may in places be as low as 
about 150 m/sec. due to aeration. It must also 
give continuous information regarding the 
configuration of the rock surface all along the 
profile, and especially underneath the shot- 
points. The field technique described above 
is the result of these considerations. From this 
shooting scheme it follows that for each shot 
point a travel-time graph is obtained in each 
direction and of such a length that impulses 
from the bedrock surface are obtained over 
a sufficiently long distance. In that way the 
“bedrock-sections” of the travel-time graphs, 
shot in the same direction from two adjacent 
shot points, overlap each other for a certain 
distance. The same holds also for the sections 
of the travel-time graphs belonging to the 
deeper layers of the overburden. In this manner 
continuous observation of all layers is obtained 
all along the profile, thus facilitating the inter- 
pretation of the travel-time graphs. This 
shooting technique has proved to be satis- 
factory with regard to interpretation at the 
same time being very efficient. 

On land a seismic crew, consisting of one 
observer and 6 helpers, covers on the average 
a profile length of about 150 m per day, 
corresponding to 750—900 m length of travel- 
time graphs per day. This involves recording 
of 25—30 shots for 5 to 6 seismometer set-ups. 
Under favourable conditions an efficiency has 
been reached of 400 m of profile, corresponding 
to 2000 m travel-time graphs per day, over 
several days in succession, with up to 60 shots 
per day. On lakes, about the same efficiency 
is usually reached at the cost of 1 or 2 more 
men who prepare the wire arrangement. In 
work across rapids usually one profile (of say 
8o—100 m length) is surveyed per day. In one 
such instance a seismic survey across a river 
was made with the object of determining the 
width of a shear zone that was expected to run 
along the center of the river. The investigation 
gave depths to bedrock in the center of the 
shallow river of some so feet, where it had 
previously been assumed that the bedrock was 
“outcropping”. Subsequent drilling fully con- 
firmed the picture obtained from the seismic 
work. The first two of these bore holes took 
five weeks to complete, each hole naturally 


giving the depth at one single point only. At 
the seismic investigation each traverse took 
only one day to complete, resulting in a 
complete bedrock profile across the river, and 
s such traverses could be measured at the cost 
of one single bore hole. 

The accuracy of these seismic depth-to- 
bedrock determinations is largely dependent 
on local conditions in regard to the composition 
of the overburden and the configuration of the 
bedrock surface. The data collected from the 
field work described above are sufficient for 
an accurate interpretation. However, it goes 
without saying that a good deal of experience, 
and a sufficient knowledge of physics and 
geology, is required from the “‘interpreter” 
of the field results, if the best accuracy possible 
is to be obtained. In this calculation work it 
is necessary that the interpreter can at all times 
feel sure about the measurements made with 
the instrument equipment. With the equip- 
ment described above this is positively the case 
—as a result of criticism and suggestions from 
field crews it has been continually improved 
during 8 years of field work under partly very 
difficult conditions, until it is now considered 
“fool-proof”. 

From excavations made and drill holes put 
down after seismic depth-to-bedrock determi- 
nations in Sweden it has been found that the 
accuracy of the seismic depth determination 
can be put at about + 10% of the depth, or 
to +1 m at depths smaller than ro m. In 
places with “normal” and fairly favourable 
conditions the difference between the seismic 
depths and the results of subsequent drilling 
usually is only a few decimeters. Where the 
bedrock surface has very strong relief, which 
often happens in Sweden, it must be kept in 
mind that the depth obtained from a seismic 
survey represents the average along a certain 
length of profile, while a drill hole gives the 
depth in one point only. In such a case a drill 
hole 10—15 m off the profile may show a 
depth that appreciably differs from the seismic 
figures along the profile. 

Fig. 8 shows a typical case of a seismic 
depth-to-bedrock determination at a proposed 
dam site. The greatly increased demand for 
this type of seismic work in Sweden in the 
last few years shows that the method is now 
looked upon by Swedish civil engineers as a 
convenient and reliable tool for quick recon- 
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Fig. 8. Plan and section showing results of seismic investigation and of subsequent boring operations. Proposed 
dam site in northern Sweden. 


naissance surveys of proposed sites for dams 
and tunnels. 

The above equipment is now also being used 
for a very special purpose, namely for deter- 
mination of the thickness of glaciers. The 
amplifier set shown in fig. 4, however, has 


then been replaced by an amplifier set with 
reflection amplifiers, including automatic gain 
control, variable filters and “mixing” (fig. 9). 

This equipment has been put together for 
the Norwegian-Swedish-British Antarctic Ex- 
pedition, which has just now started out. At 
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Fig. 9. Equipment used at the measurement on the main 
glacier of the Kebnekaise. At left the amplifier set with 
6 reflection amplifiers, in the center the recorder with 
oscillograph, timing device and camera, at right the 
power unit and in the background an accumulator 
battery and the »fool-proofing» device for firing of shots 
with simultaneous starting of film in the recorder. 


first a test was made on Hardangerjökul in 
Norway. From the experiences made during 
this test some modifications and special features 
have been included in the equipment and a 
final test has been made on the main glacier 
in the Kebnekaise mountains in North Sweden. 
One of the main objects with this test was to 
train the personnel of the above expedition in 
the handling of the instruments. 

The equipment was sent by rail to Kiruna. 
From there part of the equipment was taken 
by airplane to Kebnekaise where it was 
dropped by means of parachutes. Other parts 
of the equipment were transported by motor 
cat, boat and pack horse to the rim of the 
glacier, and from there transported on back 
and by sledge to the centre of the glacier. 
(Fig. 10.) 

Where the ice was covered by 1—2 m of 
firn it was found that charges of 3—5 gr dyna- 
mite, exploded in drill holes at a depth of 1—2 
m would give good reflections. Where the 
firn cover was very thick, charges of up to 
60 gr of dynamite had to be used, fired in 
drill holes of 2—3 m depth. For refraction 
work, with profile lengths up to 600 m, it was 
found that charges of up to 1 kg dynamite had 
to be used, but that transmission of elastic 
vibrations to the firn could then easily be 
obtained if the charge was placed at the centre 
of a semi-spherical cavity in the firn surface, 
of about 1 m diameter. 

From the refraction and reflection seismo- 


grams obtained from these measurements it 
was found that the velocity of transmission 
was 3450 m/sec. for the longitudinal waves 
and 1620 m/sec. for the transversal waves. 
Reflections were obtained for both types of 
waves. The optimal distance for these reflec- 
tions seems to correspond to the largest possible 
angle of incidence, which of course must be 
smaller than the angle for total reflection for 
the interface between the ice and the bedrock 
underneath. 

In one of the places tested, where the firn 
cover was fairly thin, reflections were obtained 
from both longitudinal and transversal waves. 
A calculation of the depth of the ice at this 
point gives for both reflections the same thick- 
ness of 96 m. One pronounced reverberation 
(repeated reflection) appears on all of the best 
records. 

At another testing place, where the firn 
cover over the ice was evidently very thick, 
no transversal waves and no reverberations 
could be observed. Good reflections of longi- 
tudinal waves were obtained, however, and 
the depth evaluated from these reflections, and 
from the velocities determined by the refrac- 
tion method, was 190 m. The probable reason 
for the absence of transversal waves in this 
case is that the firn has too small rigidity 
modulus to carry appreciable amplitudes of 
distorsional waves. 

Although these tests on the main glacier of 
the Kebnekaise had only been planned as 


instrument check-ups and not as measurements 


Fig. 10. Transport of the seismic equipment on the main 
glacier of the Kebnekaise. 
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of the thickness of the ice, their results are of 
considerable interest. They show that, by the 
use of modern equipment, ice thicknesses of 
less than 100 m can be measured by reflection 
shooting. As far as we know, these depth de- 


terminations (by R. TEGHOLM) on the Kebne- 
kaise main glacier represent a new record in 
the determination of small depths by means of 
the reflection method. 
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Abstract 


This article gives a description of a new type of recorder for ionospheric sounding. 
The recorder covers the frequency spectrum 1 to 20 Mc/s in 30 seconds and is for that 
reason of great value for investigations in the polar regions, where ionospheric conditions 
are very fluctuating. The use of a wide band transmitter and receiver has made the short 
sweep time possible. Only one variable capacitor rotating at low speed is necessary. Auto- 
matic tracking between receiver and transmitter is secured through the use of a heterodyne 
system. Samples of records obtained at Kiruna since the beginning of observations in July 
1948 are shown. 


Ionospheric research has been carried out by through the frequency sweep. On the other 
the Research Laboratory of Electronics since 1942 hand, if the transmitter is followed manually 
at our Askim Observatory, 6 miles south of through the frequency range on a receiver, a 
Gothenburg. Among other projects we have frequency band of say 15 Mc/s can not be 
been interested in sweep frequency recording covered ina shorter time than about 10 minutes. 
and have used both automatical and manual The use of tuning capacitors in the transmitter 
recording apparatus (RYDBECK, 1942 a, b). also determines an upper limit of the output 
A corner of the recording room is shown in pulse power of the transmitter. One method 
fig. 1. The picture shows equipment for both to increase this upper limit is to use oil 
fixed frequency and sweep frequency re- as dielectric but this will be a complicated 
cording. solution. 

In fig. 2 is further shown an automatic To avoid these difficulties, the most natural 
mechanically tuned sweep frequency recording solution is to use wide band stages instead of 
receiver. Although this equipment has rendered tuned stages. The panoramic recorder em- 
quite satisfactory service, we found already at ploys both a wide band transmitter and an 
the beginning of the last war that the develop- untuned receiver. The transmitter and the 
ment of a high speed “electronic” sweep receiver are controlled through a heterodyne 
panoramic recorder was highly desirable for system. Thus the only rotating capacitor in 
the further development of ionospheric re- the receiver is the tuning capacitor in a com- 
search. Due to the war situation however it mon oscillator. As this capacitor is on low- 
was very difficult to obtain quality parts for voltage level and consists of only a few blades 
the new equipment, the completion of which there are no mechanical limits for the rotating 
therefore was considerably delayed. speed. The only speed limit of the frequency 

In general the usual disadvantage with sweep is the maximum permissible detuning 
ordinary sweep frequency recorders is their between transmitter and receiver when the 
mechanical tuning elements. These make it most delayed echoes arrive. Allowing a maxi- 
comparatively difficult to construct receivers mum permissible detuning of 8 kc/s when 
that automatically follow the transmitter echoes arrive from a virtual height of 600 km 
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Bios. 


A corner of the recording room at the Askim 
Observatory. 


a minimum time of 10 seconds is due for 
each frequency sweep. 


However, the single remaining rotating 
capacitor could also be omitted. A reactance 
tube could be used as a variable common 
oscillator. But this would probably call for the 
use of an intermediate frequency that would be 
unnecessarily high. 


The panoramic system we have used was 
suggested already in 1939 by WıLLıam ©. REED 
and ©. E. H. Rypsecx, both at that time at 
Harvard University. The recorder covers the 
frequency spectrum from 1 to 20 Mc/s in 
30 seconds, with an output pulse power va- 
rying from 16 kW on the lowest frequency 
to 6 kW on the highest frequency. The result 
is presented on the screen of a 12” cathode 
ray tube and recorded with a 16 mm ciné 
camera. This recording equipment has been 
in continuous operation at the Kiruna Ob- 
servatory since the summer of 1948. 


During the course of the development of 


our recording equipment information was 
received from the United States regarding 
a very similar outfit developed by P. G. SUL- 
ZER and associates (2948). Surzer’s work pro- 
ved to be very valuable and stimulating com- 
parison in our efforts to develop a high quality 
recording equipment. 


The Recorder 


The block diagram is shown in fig. 3. From 
this diagram is seen that the frequency deter- 


~ 


mining part in the equipment is a common 
oscillator, variable from 31 Mc/s to so Mc/s. 
The output of this oscillator is mixed with the 
output from a fixed oscillator operating on 
30 Mc/s, the result being an output voltage 
with a frequency continuously varying from 
1 Mc/s to 20 Mc/s. This radio-frequency is 
amplified in an untuned wideband amplifier 
and fed to the transmitting aerial. 


The common method when making iono- 
spheric investigations is to transmit hae pulses 
of radio- frequency and to receive these after 
they have been reflected by the ionosphere. 
By measuring the time delay of the reflected 
pulses one can determine the virtual height of 
the reflecting layer. In this recorder the 30 Mc/s 
oscillator, the driver, and final stages are pulsed 
by a 50 ee pulse. The pulse repeti- 
tion rate is 50, synchronized with the power 
line. The output from the variable common 
oscillator is also fed to the receiver input stage, 
where it is mixed with the incoming reflected 
signals. The difference frequency remains 
constant, 30 Mc/s, throughout the entire 
frequency band. The first intermediate fre- 
quency of the receiver is therefore chosen as 
30 Mc/s, which means that the transmitter and 


An automatic mechanically tuned receiver for 
sweep frequency recording. 
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the receiver always are tuned to the same 
frequency. 

The output voltage from the receiver, which 
contains both the transmitted ground wave 
and the echoes is fed to the indicator and 
applied as intensity modulation of a 12” cath- 
ode ray tube. By chosing a suitable vertical 
sweep, synchronized with the pulse frequency, 
and a horizontal sweep, proportional to the 
frequency of the transmitter, the result will be 
presented on the screen of the indicator tube 
in the form that is common for sweep fre- 
quency records. The time it takes to form the 
image of each sweep is 30 seconds, and as the 
indicator tube has an afterglow time of 30 
seconds each sweep can be observed visually. 
The various pulse and gate voltages that are 
needed in the recorder are generated in the 
time-base unit. 

A short description of the different units 
that form the panoramic ionospheric recorder 
will be given below. 


Common oscillator 


The tuning capacitor in the common oscilla- 
tor is made with semi-circular plates and is 
slowly rotated by a synchronous motor at a 
speed of one rpm. Thus, during one rotation 
of the capacitor, the frequency increases from 
its lowest value, 31 Mc/s, to the highest, so 
Mc/s, and then back to the lowest again. This 
calls for a very high degree of precision in the 
construction of the capacitor, the same fre- 
quency is obtained with two different settings 
of the capacitor, and those two settings must 
make the same angle with the lowest frequency 
setting. This is due to the fact that on the same 
shaft as the capacitor is mounted a precision 
linear potentiometer, over which is developed 
a dc-voltage proportional to the position of 
the tuning shaft. This voltage is zero when 
the oscillator frequency is 31 Mc/s and in- 
creases to + so volts on the highest fre- 
quency. The horizontal sweep of the indicator 
tube is driven from this source, the horizontal 
position of the electron beam will therefore 
be proportional to the frequency of the trans- 
mitter. 

The oscillator circuit is a Hartley oscillator, 
and the radio frequency voltages to the trans- 
mitter and the receiver are separated by a 
buffer stage. This was necessary in order to 
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avoid feed-back between the transmitter and 
the receiver through a common connection. 

In the same unit as the common oscillator 
is also mounted the input mixer tube for the 
transmitter. This was necessary in order to 
decrease the frequency that had to be fed to 
the transmitter unit. The local oscillator for this 
mixer is a pulse-operated miniature triode, 
6C4, working on 30 Mc/s. The output of radio 
frequency pulses is fed to the injector grid of a 
6K8 mixer tube and the signal grid of the same 
tube is fed by the common oscillator. Over the 
plate load resistor will therefore appear a pulsed 
radio frequency, varying from I to 20 Mc/s 
as the tuning capacitor covers the frequency 
band. These pulses are then fed to the wide- 
band amplifier. 


Transmitter and antenna system 


The transmitter is a broad-band amplifier 
with a band-width of about 20 Mc/s, using the 
miniature pentode 6AKs in the first three stages 
(FINK 1940). These are followed by another 
two stages of amplification, using a 6AG7 and 
a 807. The next stage is a phase-inverter stage 
driving a 4—250 A push-pull stage. This stage 
is operated with a plate voltage of 1,500 volts 
and a grid bias of — 200 volts. During the time 
when the pulsed 30 Mc/s oscillator operates, 
this bias is decreased to — 25 volts. Each tube 
draws about 0.4 A in plate pulse current. From 
this stage, the radio frequency is fed to the 
final stage, two 304TH connected in push-pull. 
These tubes operate with a plate supply voltage 
of 5,000 volts and a bias of — 1,000 volts. 
They are also pulsed, during the pulse the bias 
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The transmitting and receiving D-antenna 
used with the panoramic ionospheric recorder. 


Biss 74). 


is reduced to —- 20 volts, corresponding to a 
plate current of 2.5 A. Continuous operation 
under these conditions would of course be out 
of question. The maximum pulse output power 
is 16 kW, but a new version of this recorder 
is planned with an output power of 100 kW. 
The use of a wide-band amplifier requires 
more stages than an ordinary tuned amplifier, 
but on the other hand the problem of ob- 
taining high output power is much more 
simplified. 

One of the more important problems with 
a transmitter covering such a wide band as in 
this case is that of finding a suitable antenna 
system. The impedance of an ordinary antenna 
depends very much on the frequency, showing 
extreme values at the resonant frequencies. It 
is not possible to use such an antenna because 
the load characteristics should vary too much 
within the frequency band. We have therefore 
tried combined antennas of different kinds, 
and the antenna that seemed to give the best 
result was a construction used by the Depart- 
ment of Terrestrial Magnetism, Carnegie Institution 
of Washington (GEORGE 1945). This so called 
D-antenna has an average impedance of 600 
ohms, and the impedance does not change 
more in amplitude than 2 : r-within the en- 
tire frequency band. A drawing of the D- 
antenna is shown in fig. 4. As seen in this 
figure there are two antennas mounted at right 
angles. One of them is used for the transmitter 
and one for the receiver. We have not 
completed our tests to find an even better 
antenna. Further improvements are therefore 
to be expected. 


Receiver 


The input stage of the receiver 1s a miniature 
double triode, 6J6, with a 1,000 ohms resistor 
in the common cathode return. The radio 
frequency from the common oscillator is 
injected over this resistor and the grids are fed 
out of phase from the antenna. The plate circuit 
is tuned to the intermediate frequency, 30 Mc/s. 
After the mixer stage follow two transformer 
coupled stages of amplification. In the next 
stage this intermediate frequency is changed to 
a new intermediate frequency, 1,600 kc/s. The 
local oscillator is a Hartley oscillator, but in 
order to improve the stability this oscillator 
will be replaced by a crystal oscillator. After 
one stage of amplification the signals are recti- 
fied and passed through a limiter circuit in 
order to prevent the transmitted pulse and 
strong echoes to overload the cathode ray tube. 
In order still to improve this matter the receiver 
will be redesigned to have a logarithmic vol- 
tage response. 


A receiver for reception of these short pulses 
must have adequate band-width to pass the 
pulses without distortion. Special circuits have 
therefore been incorporated to fulfill this 
requirement. Another thing of importance is 
that the different parts of the receiver must be 
very well screened from each other. For 
example, even a small voltage leaking from the 
local oscillator to the input stage causes cross- 
modulation with strong broadcasting stations 
resulting in unnecessary interference lines on 
the recorded picture. 

The sensitivity of the receiver is better than 
4 microvolts over the whole band except for 
the lowest frequency where it decreases to 10 
microvolts. The receiver assembly is shown in 
Lig. 5. 
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Time-base and indicator 


The time-base, fig. 6, serves as a generator of 
the various pulse and control voltages that are 
required. The functions of the time-base are 
synchronized with the frequency of the power 
supply. It delivers pulse voltages for the 30 
Mc/s oscillator, the driver and the power 
amplifier stages in the transmitter. Further it 
delivers a saw-tooth voltage for the vertical 
deflection in the indicator tube. The duration 
of this voltage is variable from o up to a time 
corresponding to the time it will take for the 
signals to travel up to a height of 1,000 km and 
back again. There is also delivered a negative 
gate pulse to the cathode of the indicator tube 
with a duration equal to that of the saw-tooth 
voltage. This voltage intensifies the spot on 
the screen during the sweep. 


During the vertical sweep there are also 
delivered short positive pulses to the grid of 
the indicator tube. These pulses occur with a 
frequency of 3,000 cps. They will therefore 
serve as markers for each 50 km in virtual 
height. If there is no deflecting voltage on the 
horizontal plates of the indicator tube these 
pulses and the pulses from the receiver will 
appear as dots on a vertical line on the screen. 


But during the frequency sweep the voltage : 


from the linear potentiometer in the common 
oscillator will make the beam move horizon- 
tally over the screen. Then the markers and 
the echoes will appear as horizontal lines on 
the screen. The broadcasting bands will appear 
as vertical lines. 


In front of the screen of the indicator tube 
is placed a 16 mm Paillard-Bolex motion- 
picture camera. Dhisn cans besseensin fig... 7; 
which shows some units of the recorder. The 
shutter of this camera remains open during 
each frequency sweep. After the sweep is com- 


The time-base used for generating the various 
control voltages in the recorder. 
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pleted the shutter closes for a moment to allow 
the camera to advance a new frame of film. 
The record will therefore be a time exposure. 
The operation of the camera is made by means 
of a mechanical device driven by a solenoid, 
which is energized through a micro switch in 
the common oscillator. In addition to the 
frequency sweep a small watch showing local 
time is photographed on the same picture. 
This identification system will in the near 
future be completed with another marking 
showing the date of the day. 

The operation of the recorder is wholly 
automatical, a system of relays is incorporated 
to make the recorder start and stop at pre- 
determined times. A precision clock controls 
the starting pulses for this system of relays. The 
relay panel is shown in fig. 5. 

Usually the recorder starts every hour and 
records six pictures, but it is also possible to let 
the recorder run continuously day and night. 
Normally the ionosphere does not change so 
very much from picture to picture. Therefore 
it is possible to play the recording in an ordinary 
movie projector and study the gradual changes 
in the ionosphere. Especially during magnetic 
disturbances when ionospheric conditions are 
changing very rapidly, these films are of very 
great importance for investigations of the 
character of the disturbance. 


Records 


„In fig. 8 and fig. 9 two typical records are 
shown obtained with the panoramic recorder 
at Kiruna during tne spring 1949. 
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Fig. 8. Record made on an undisturbed day, April 


25, 1949, at 0900 GMT. 


Fig. 8 shows a typical record made on an 
undisturbed day. The echoes start at the lower 
frequency limit of the recorder, appearing at 
a height of roo km which is the normal height 
of the E-layer. At 2.1 Mc/s the virtual height 
seems to increase, and at 3.45 Mc/s the critical 
frequency of the E-layer is reached. At 3.3 
Mc/s echoes also arrive from the Fr-layer at a 
height of 200 km. The echoes are split into 
two components from 4.4 Mc/s. These are the 
ordinary and the extraordinary components 
of the wave due to the magneto-ionic effect 
of the terrestrial magnetic field. The critical 
frequency of the ordinary wave is 5.0 Mc/s. 

Echoes from the F2-layer appear at 5.7 Mc/s 
and these are also split in two components. The 
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Fig. 9. Record made on a disturbed day. 
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A sweep frequency Kerr-cell recording, 
July 18, 1945. 


Fig. 10. 


ordinary wave has a critical frequency of 6.4 
Mc/s and the extraordinary disappears at about 
7.1 Mc/s. The vertical lines seen on the picture 
are broadcasting stations within the frequency 
band. 

Fig. 9 on the other hand shows a record made 
on a magnetically disturbed day. There are no 
distinct layers present. The reflections from 
t.s Mc/s up to 9.1 Mc/s in the height interval 
90 km to 320 km are all probably due to ioniza- 
tions caused by corpuscular bombardment. 
The F2-layer also shows a very diffuse character 
with a critical frequency of approximately 
4.6 Mc/s. 

In fig. 10 is shown an example of a Kerr-cell 
recording (RYDBECK 1945, 1946), obtained 
on July 18, 1945 by the Chalmers’ Solar Eclipse 
Ionospheric Expedition. This record is obtained 
by a manual sweep frequency recorder. The 
time required for this frequency sweep was 
six minutes, the corresponding time for the 
panoramic recorder is 15 seconds. However, 
the example shows a very good proof of the 
quality of the Kerr-recorder. 

Finally in figs. 11 and 12 are shown two 
records obtained on May 5, 1949, fig. 11 at 
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Fig. 11. Record made on May $, 1949 at 2142 GMT. 
2142 GMT and fig. 12 at 2143 GMT. The 
time elapsed from the end of the first record 
to the start of the second is one minute. 

These two records give a very good example 
of how rapidly the conditions in the iono- 
sphere often vary north of the arctic circle. 

The scattered reflections at 2142 GMT 
within the height interval 400—550 km have 
disappeared at 2143. What remains in that 
interval is a multiple reflection of the F2- 
layer. Instead there appears several new echo 
patterns in the height interval 100—300 km. 
It is also interesting to note how the recorded 
signal strengths of the various traffic bands 
vary. 

These records are shown only as examples 
of the performance of the panoramic iono- 
spheric recorder, especially for investigations 
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Fig. 12. Record made on May $, 1949 at 2143 GMT. 
of rapidly changing ionospheric phenomena. 
Detailed analysis and discussions of the many 
records obtained at the observatory will appear 
shortly in special reports from the Research 
Laboratory of Electronics. 
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Abstract 


It is assumed that in horizontally uniform zonal currents certain long waves can grow in 
amplitude whereby the kinetic energy associated with the growth is drawn from the potential 
energy of the basic current system. This restriction as to the energy redistribution permits to 
treat the simultaneous physical changes of state as polytropic, which simplifies considerably the 
integration of the relevant equations of motion. It is shown that the amplitude growth, at a 
prescribed zonal dimension of the perturbation(wave-length), is primarily controlled by the 
value of a non-dimensional quantity, viz. the RICHARDSON number (= static stability divided 
by the square of the vertical shear of wind), in a similar manner as this number controls the 
time variation of turbulence to the extent of its being a result of energy redistribution. The 
derived criterion for instability predicts that all but the very long waves in the middle latitude 
westerlies are unstable. 

The intensity of the absolute circulation, in addition to the vertical shear, has an important 
bearing on the considered phenomenon of instability. This is discussed with regard to the non- 
equilibrium position of a cold anticyclonic vortex centered at the pole. The stability criterion 
is also scrutinized as to its validity at tropical latitudes. 


Controlling the Instability of Long Waves in Zonal Currents’ 


I. Introduction 


Considerable attention is at present given 
to the origin and life-history of the long, 
undulating upper level disturbances character- 
istic of the middle latitude belt of westerlies. 
To the extent that these disturbances may be 
considered as quasi-horizontal wave motions 
imposed on autobarotropic currents, their 
shape, wave-length and displacements are 
controlled primarily by the variation of the 
Coriolis parameter with latitude. 

In view of recent investigations into the 
breaking up of zonal motion (BERGGREN, Bo- 
LIN and Rosssy, 1949) one must accept the 
existence of some mechanism by which the 
planetary waves are dynamically linked to 
the minor scale disturbances of the frontal 


1 To be published shortly in the Communications, 
Series B, Sveriges Meteorologiska och Hydrologiska 
Institut. 


wave type. However, the breaking up of a 
quasi-zonal current into large scale vortices 
may be regarded, without reference to the 
existence of frontal waves, as a phenomenon 
of long-wave instability on the grounds of a 
rather marked difference in the horizontal 
scale of the respective disturbances. Theore- 
tical investigations — HAURWITZ (1940), JAW 
(1946) and CHARNEY (1947) — have dealt 
with the dynamics of horizontally uniform, 
barocline (vertically shearing) zonal currents. 
These theories demonstrate that the baroclinity 
and the static stability of the basic current 
system are the principal factors in the early 
stages of the life-history of the perturbation. 

The object of the investigations summarized 
below is an attempt to incorporate, in a 
general manner, the physical process into the 
theory for long waves (though restricted 
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by the same fundamental assumptions as 
those made by the above authors); to derive, 
by using relatively simple mathematical pro- 
cedure, such elementary instability criteria 
that can be more readily interpreted in the 
light of the energetical and dimensional aspect 
of the phenomenon treated here: large-scale 
turbulence appearing in a flow which initially 
is large-scale laminar; finally, to develop the 
elementary criteria with a view to attain 
more conformity with empirical evidence. 


2. A limiting case of energy redistribution] 
(Introduction of the polytropic invariant) | 


To the extent that the kinetic energy of 
the basic zonal current remains unchanged 
during the growth of the disturbance, the 
growth phenomenon is merely a manifesta- 
tion of a release, through the perturbation 
motion, of potential energy. Evidently, the 
existence of a barocline field (more specific- 
ally of solenoids in the meridional plane) in 
the basic system is a prerequisite for availab- 
ility of that part of potential energy which 
must be converted into kinetic energy of the 
perturbation motion. These specifications re- 
garding the energetic aspect of the process 
imply that the current system is thermally 
inactive and the physical process associated 
with the motion is characterized by polytropic 
changes of state. The polytropic process is a 
reversible thermodynamic path of constant 
heat capacity c’, so that the heat supplied to the 
system per unit of time to unit of mass is de- 
termined by the differential relationship dQ = 
=c’dT where T is the absolute temperature. 
To each of the polytropic classes (isobaric, 
isothermal, (dry) adiabatic and isopycnic) one 


can assign an individual invariant £ (conserva- 
tive property) or, in the mathematical formu- 


lation 3 = 0, so that e.g. in the case of dry- 


’ dt 
adiabatic process E is being identified with 


the potential temperature ©. It is, however, 
possible to treat, subject to subsidiary restric- 
tions as to the distribution of water-vapour, 
a reversible saturated-adiabatic change of 
state as a quasi-polytropic process so that 
c’ is also in this case very nearly constant. In 
doing so, allowance can be made for the 


4— 904613 


effect on energy conversion of the formation 
of extensive cloud masses which may accom- 
pany the growth of the baroclinic wave. 


3- Derivation of a basic frequency equation 


A first approximate solution for the phase 
velocity of plane waves, satisfying the hydro- 
dynamic and all other relevant differential 
equations as well as a minimum of boundary 
conditions, yet revealing a fundemental cri- 
terion for instability, is obtained in the follow- 
ing way. The theorem for the path variation of 
the vertical component of vorticity is trans- 
formed (and linearized) for “small perturba- 
tions” propagating with constant speed in a 
transversally uniform zonal current of infinite 
width. Making use of the theorem of geo- 
strophic meridional perturbation velocity!, a 
linear differential equation of effectively second 
order with the dependent variables p (pressure 
perturbation) and & (perturbation of the po- 
lytropic individual invariant &) is derived. 
Two solutions are found, one establishing a 
frequency relationship for short waves whose 
wave-length is of the order of magnitude 
equalling the depth of the atmosphere, the 
other yielding an equation for the phase 
velocity of long (planetary) waves. The waves 
of the first group are linked with the dis- 
turbance of hydrostatic equilibrium and are 
capable of growing, in certain conditions, if 
the static stability is negative. The long waves 
are linked with the disturbance of zonal 
balance. To understand and measure the 
extent of validity of the criteria for instability 
of the latter waves, the reader is asked to 
follow a few steps in the derivation of the 
principal equation for the phase velocity. 

In order to simplify the mathematical aspect, 
the effect of compressibility of the air will be 
disregarded throughout in this summary. 
Little generality is lost through this expedience 
as far as the main results are concerned. In 
view of the restrictions with regard to the 
energy redistribution, as mentioned earlier 


t The assumption that the meridional velocity is 
geostrophic or very nearly geostrophic, involves the 
exclusion of long gravitational waves from possible 
solutions of the equations of motion. RossBy (1939) has 
shown that in principle this must be so and later CHARNEY 
(1947) has given a mathematical proof of this circum- 
stance. 
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in this paper, it must be expected that the 
behaviour of the perturbation, i.e. its speed 
of propagation and its time-rate of growth (de- 


termined by the real and imaginary parts of 


the phase velocity c resp.) depends on a) 
steady state or “basic” quantities such as the 
intensity of the zonal circulation u or the 
meridional slope of the isobaric surfaces mp; 
the rate of basic baroclinity which in con- 
junction with m» is determined by the me- 


ridional slope of the constant &-surfaces, me; 
b) parameters constituting the structure of the 
perturbation such as wave-length, L, and 
phase angle ® as well as the amplitudes po and 


& of the respective perturbations in the 


pressure- and é-field; c) planetary and physical 
constants as e.g. the geographical latitude @ 
(central to the zonal current) and the radius of 
the earth, a; the acceleration of gravity, g. 
A general theory of atmospheric waves 
would have to allow for variations of all these 
quantities, except a and 8; with respect to 
latitude, longitude and elevation whilst within 
the limits set by the assumptions made already 
above, we have to allow for variations of u 
(or mp), me, po and & as well as the phase angle 
with elevation. There are no physical or 
kinematic boundary conditions which impose 


a restriction on the &-surfaces whose meridional 
slope may therefore without loss of generality 
assumed to be uniform. It is well known that 
the variations of the amplitudes po and & with 
elevation are subject to the boundery condi- 
tions and it is this circumstance which makes 
the integration of the relevant differential 
equation so difficult. The view is here put 
forward that among the factors that primarily 
govern the law of amplitude growth, the 
height-dependency of the perturbation am- 
plitude is of secondary importance. We shall 
therefore treat the amplitudes po and & as 
constants for the integration of the transformed 
vorticity equation.' We obtain 


' To further simplify the problem we shall, for the 
time being, disregard the tilt of the component waves 
in the vertical, in other words, treat also the phase-angle ® 
as independent of height. Later on the effect of the inclu- 
sion of the tilt will be discussed. 


~ 


where u* is thus the zonal wind velocity 
relative to the moving wave (i.e. the “velocity 
deficit”) and ß/k? its value in the case of non- 
divergent waves in a barotropic, homogenous 
zonal current (RossBy 1939). In this equation 
f denotes the Coriolisparameter, ß its rate of 
change northward, H the height of the homo- 
genous atmosphere whose surface temperature 
is equal to that at the reference level, k = 2/L 
the wave number and ¢ the lag of phase be- 


tween the &-wave and the pressure wave.‘ The 
second term on the r.h.s. gives the combined 
effect of horizontal divergence and horizontal 


solenoids in the p- and &-fields on the net acce- 
leration of relative vorticity — this effect 
expressed as a velocity —, whilst the preceding 
term, as is well known, gives in the same 
manner the effect of the planetary or “induced” 
variation of vorticity. 


Equation (r) is as yet underdetermined with 
respect to u* or c, since the amplitude ratio 


bo, M 

A is not expressed in terms of known quantities, 
and since the phase lag e is in general not 
determinable without regard to the vertical 
tilt of the two component waves, this tilt 
having been disregarded here so far. From 
general principles it can be shown that the 


or 

non-dimensional factor A = 2 ah cos ke 
po € 

depends upon the meridional horizontal 


; 0 & 
gradient of the conserved property, . , and 
Oy 


the orbital frequency »* and thus on wave- 
length L and velocity deficit u*, since always 


ERST 
pr = 


aL In particular, in a plane at which 


the motion is horizontal (nodal plane of verti- 
cal motion) the factor A is easily determined 
from the developed and transformed law of 


. * Û 
conservation, VIZ. — 


Ser which had been 


' The bar over the quantities p and & indicates that 
the value in the undisturbed state is taken. 
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discussed earlier in this paper. One finds at 
a nodal plane 


SH 
Ao = — = 5 
ur f (2) 
ht: f Ta 
where S* = 2 —,a quantity which, in analogy 
E Ô y _ 

| er tae 
to the static stability S = & 37 may be term- 

OZ 


ed “dynamical stability function”. Substituting 
(2) in (1) we obtain for the phase velocity c the 
equation 


CETTE 


p [4% (Se) 
> L2 LE \ Po B? Se | > (3) 


where io is the zonal velocity at the nodal 
plane. In the derivation of this equation we 
had assumed that S:* and S: are independent 
of height or, in other words, that the me- 


ridional slope of the iso- £ -surfaces is uniform 
throughout. This implies that e.g. in the me- 
teorologically relevant classes of polytropic 
change of state, viz. the dry- and saturated- 
adiabatic, (£ = © or ©, resp.) there is a free 
surface except if the absolute temperature 
does not vary or increases with height. Eq. (3) 
gives then the speed of propagation of simple 
harmonic waves in a zonal current with one 
nodal plane of vertical motion and one free 
surface provided the current is horizontally 
uniform, the vertical and (meridional) hori- 
zontal gradients of the quantity & are constant 
throughout. The adopted model does therefore 
not allow for the existence of a tropopause 
0€ 
oy 
tinuities. Neither does it allow for inho- 
mogeneity, e.g. within the troposphere, of the 
concentration of isentropes in horizontal 
surfaces which is frequently observed in even 
otherwise well organized zonal currents. And 


Ô k 
surface at which both — and = possess discon- 


ı For a barotropic current this relation, after some 
transformation, becomes identic with RossBy’s formula 
(1942) for the ratio between the amplitude of isotherms 
(iso-lines of conserved property) and the amplitude of 

; A u 
streamlines, et 


A yx 


the variations of the verticallapse-rate of tem- 
perature in both, the vertical as well as the 
(meridional) horizontal direction, which may 
be of importance, cannot be introduced as 
significant quantities in eq. (3). Now, disre- 
garding here the effect of the earth’s oro- 
graphy as well as that of frictional stresses, 
the earth’s surface is regarded as a nodal 
plane in the vertical motion, i.e. a horizontal 
nodal surface. In eq. (3) # represents then the 
surface geostrophic wind of the basic current 
system. 


4. The elementary instability criterion and 
its validity in the middle latitude westerlies 


An analysis of the equation (3) relating 
the phase velocity c to the intensity of the 
zonal circulation at the ground wo, the 
wavelength L (through k), the geographical 
latitude @ (through 8) and the stability func- 
tions S;* and Ss (henceforth in application 
to dry-adiabatic processes denoted by simply 
S* and S) shows briefly the following. A neces- 
sary condition for amplitude growth (ex- 
istence of an imaginary component c;) is that 
the static stability S = o, higher static stability 
values involving less rapid, lower stability 
values involving more rapid amplitude growth. 
This dependency on static stability appears 
physically quite in order, since an existing 
rate of static stability must damp out those 
vertical components of the wave motion 
which occur everywhere above the nodal 
plane and are a consequence of the requirement 
of mass continuity. In the limiting case when 
the static stability vanishes, i.e. at so-called 
“indifferent vertical stratification” of the fluid, 
the solution eq. (3) breaks down except 
when simultaneously the meridional in- 


equalities in the é-field (potential temperature- 
or entropy field) vanish as well.“ Because of 
this break-down at the singularity S =o 
(S* + 0), it cannot be expected that in the 
negative regime of S the relationship holds 
valid for one and the same phenomenon as it 
does in the positive regime. In the latter regime 
we have to take the first of the solutions of 


1 In this particular case we have homogeneity in 


£ and eq. (3) reduces to the well known Rossby formula 
for planetary waves in an autobarotropic zonal current. 
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the differential equation, mentioned in § 3 
which yields a relationship for very short 
waves only, when the relevant atmospheric 
values for S* and S are being introduced. 
We shall not here discuss the properties of 
these waves. However, it is worth mentioning 
that in the normal case of static stability being 
present (S>o), these short waves can be 
identified with the phenomenon of billowing 
clouds for which Haurwitz (1941) has given 
an adequate theory. One has then to consider 
the limiting case when the baroclinity (ver- 
tical shear) is concentrated in an infinitely thin 
horizontal layer of air and the static stability 
in the adjoining infinitely deep strata is 
disregarded. 

In the following it will be assumed that the 
static stability is positive. (It is, however, 
noteworthy that to a vertical lapse-rate of 
potential temperature or entropy which’ in 
meteorological practice is regarded as “very 
small”, corresponds a value of S such as to 
yield, at prescribed L (k) and S*, still values 
for & or c; of meteorologically reasonable 
dimension). With this specification the re- 
quirement for instability may be expressed 


by the inequality: 


[S*]_ 8 _ Qceosp L 
VS 2k a 


(4) 


2 IT 


or, if it is desirable to express this condition 
through a critical (maximum) wavelength L., 
by the equation 


27a RSA] 


Q cosp VS 


L, = 


(4 a) 


or, finally, if the critical (minimum) number 
n. of waves around the circumference is 
preferred, by the equation 


5 
| SX] 


== 


ne = 2 cos? p 


(4 b) 


In these relationships 2 is the angular velocity 
of the earth rotation, a the earth’s radius and 
@ the latitude which is representative for the 
zonal current under consideration. 

It is immediately seen that when the isen- 


~ 


tropic surfaces are horizontal (S* = 0), even 
the shortest waves would be stable, whilst 
when the static stability is reduced to zero, 
even the longest waves would be unstable. To 
which extent these elementary instability 
criteria conform with empirical data, could 
be checked (though only very roughly), 
since S* and S could be computed and mapped 
for discrete zonal belts over long periods of 
time. However, it is preferable to make first 
use (as is usually done) of the gradientwind 
equation which relates the dynamic stability 
function S* to the surface (zonal) wind, the 
vertical shear of the (zonal) wind, the so-called 
thermal wind, and to the static stability. 
Equation (4 b) may then be written as 


1—sing VS/|A| 
I HS 
g|4| 


(5) 


n= — 
2 sin y | 


where the symbol A denotes the vertical 
windshear. It is noteworthy that the non- 


dimensional quantity VS/[Alis the square 
root of the Richardson number (see BRUNT, 
1934) which plays a similar rôle in the theory 
of energy redistribution in turbulent currents 
in which shearing stresses are operating 
(RICHARDSON, 1920). According to various 
investigators, this number (R) is close to 
unity. The condition that turbulence should 
increase, i.e. the rate of supply of energy by 
the shearing stresses be at least as great as 
the work which has to be done to maintain 
the turbulence against gravity, is expressed 
by the inequality R < 1, or in analogy to 
eq. (5) 


Ray (6) 


In order to develop further the analogy 
between these criteria for turbulent redistribu- 
tion of energy — the breaking up of organized 
zonal currents into vortices has been since long 
regarded as a phenomenon of large scale 
turbulence — we shall for the time being 
neglect the term io S/g|A | in the denominator 
on the r.h.s. of eq. (5), as being small compared 
with unity. In fact this is permissible in all 
“normal” (statistical average) conditions in 
the belt of middle latitude westerlies. Im- 
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portant exceptions occur when normal or 
high static stability (S = 1.3 10-4 sec-2) is 
coupled with weak shear (A = 5 - 1074 sec-1) 
and/or normal to high zonal velocity at the 
ground (uo = 10 to so m/sec). The above 
criterion (eq. 5) may be given the form: 


qo ( an sin a 
I — sin? @ 


In the middle latitude belt of westerlies, be- 
tween 45° and 65° latitude, one observes 
frequently 4 to 6 deep troughs in the streamline 
pattern which — for the sake of providing a 
point of argument — may be looked upon 
as the result of a redistribution of energy in the 
limited sense as earlier specified. With the 
values n = 5 and 9 = 55° the quantity R is 
in the neighborhood of 600. In the diagram 
fig. 1 the relationship between R, n and y is 
represented. In “normal” conditions in the 
belt of westerlies the vertical shear is of the 
order of magnitude 1 m/sec. per kilometer 
elevation and the static stability about 1.3 . 1074 
sec": (corresponding to a temperature lapse- 
rate 0.65°/100 m). With these values we find 
R = 100. On the other hand, it is seen from 
the diagram that the critical Richardson 
number is for all but the very longest waves 
considerably above the value 100 (north of 
45°). Thus if the criterion, derived with the 
use of many simplifying assumptions, is 
expected to give some part of the truth, 
extremely long planetary waves only should, 
north of latitude 45°, exhibit a large degree 
of persistency with regard to shape. To some 
degree this persistency is confirmed by synoptic 
evidence but naturally other factors not having 
been considered here, must have a great and 
frequently overruling effect. 

In the diagram fig. 2 the above criterion 
of instability is given as a relation between 
the circumferential wave number, the vertical 


(5a) 


shear and the latitude for a prescribed “nor- » 


mal” rate of static stability. For a number 
n=5 (L=72° longitude) there is good 
agreement between this elementary criterion 
and that derived by CHARNEY (1947) who 
carried out a complete integration of the 
hypergeometric confluent differential equation 
which is obtained through taking into account 
the variation of the perturbation amplitude 
with height.‘ This agreement indicates that 
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Fig. 1. Diagram illustrating the dependence of the critical 
wave-length Le (or circumferential wave-number 1c) 
upon the value of the non-dimensional quantity 
R = S/A? (RICHARDSON number) and the latitude 9, 
through four selected isopleths for ne. 


one does not introduce a really essential factor 
in regarding the physical upper boundary 


conditions. 


5. Remarks on zonal currents in tropical 
latitudes 


In the derivation of (3) and of the en- 
suing instability criteria it was not essential 
that the zonal current at the nodal plane be 
directed from west to east and the vertical 
shear be positive. Therefore these relations 
are, at least in principle, applicable to easterly 
winds and/or negative shear. To the north 
and south of the equator there are extensive 
belts of low-tropospheric easterlies in which 
the temperature decreases but slightly north- 
wards. In view of the low latitude this weak 
solenoidal field is associated with a strong 
positive vertical shear, so that the east-west 
circulation is superimposed by a west-east 
circulation. The static stability in this current 
system is on the whole rather low which 


1 CHARNEY’s numerical criterion applies to a prescribed 
(normal) static stability and a representative latitude 45°, 
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Fig. 2. Diagram showing the dependence of the critical 

wave-length Le upon the vertical shear of the zonal wind 

A and the latitude 9, for a prescribed, normal value of 
the static stability S. 


circumstance together with the strong vertical 


shear permits to disregard the term on 


Ho S 
g|4| 
the r.h.s. of (5). We may thus regard the 
elementary criterion (6) to have at least 
the same degree of validity in the latitude 
belt 10°—35° as in the middle latitude belt. 
As is evident (see fig. 1) the critical Richardson 
number is in tropical latitudes much smaller 
than 100 (and would e.g. be reduced to unity 
at latitude 9° for a circumferential wave num- 
ber 3). On the other hand the rates of static 
stability and vertical shear prevailing in “nor- 
mal” conditions in the low latitude belts of 
surface easterlies (trade winds) are such that 
the corresponding value of R exceeds for 
even relatively short waves considerably the 
critical value, though it is, in such condi- 
tions, smaller than in middle latitude westerlies. 
Thus the kind of instability considered here 
is not likely to be realised for any reasonable 
perturbation dimensions, However, in extreme 


conditions the actual value of R may fall 
below the critical value and then instability is 
possible. The present theory may possibly be 
of some use for the understanding of the 
coupling between the storm systems forming 
in or near the trade wind zone and the long 
waves in the near-equatorial upper westerlies 
which has come to light in recent investiga- 
tions (RIEHL 1948). 


It should be emphasized that in applying 
the above criterion to quasi-zonal currents 
with an east-west motion at the ground it is 
essential that the reversal to opposite motion 
(westwind) should take place at rather short 
distance from the ground. In average condi- 
tions, in the regions just to the north and 
south of the equator, this condition is not at 
all satisfield. However, as e.g. the observations 
from the Finnish aerological expedition of 
the year 1939 show (see also VUORELA, 1948), 
the average flow pattern in the belt 5°—15° 
south during the period October—November 
1939 is featured by a rather shallow layer 
with light easterlies components and by uni- 
formily increasing westerly components in the 
troposphere above, the average vertical shear 
amounting to as much as 2.5 m/sec per km. 
The static stability had an average value of 


1.0: 1074. These values bring the number R 
down to a low value of 16. According to our 
criterion this value requires (at 9° latitude) a 
circumferential wave number 12, i.e. a wave- 


length of 30° longitude. 


6. The significance of the absolute circulation 


In order to test the validity of the basic 
equation (5) which had been derived for an 
extremely simplified atmospheric model and 
under many restrictive assumptions, it is of 
importance to find out whether computed 
velocities of propagation (of the perturbation) 
and time-rates of its growth bring out some 


. of the characteristics of middle latitude disturb- 


ances. In the diagram fig. 3 the velocity deficit 
Ho — c is represented as function of the surface 
zonal wind to for different rates of baroclini 

(vertical shear A) at a prescribed latitude (45°) 
and wavelength (51° longitude), the static 
stability being assumed to have the normal 
value. The characteristic time-rate of growth, 
Ti, taken as that interval during which the 
amplitude would increase to e times its initial 
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value, is also given as function of the surface 
zonal wind for different shear values. 

It is seen that for moderate or strong surfaces 
westerlies (Yo > 0), the disturbance moves 
eastwards with a speed that differs little from 
the surface velocity, when the shear is con- 
siderable, and this correspondence is the better 
the closer the surface wind is to its critical 
value (indicated by a circle) where instability 
would set in. This circumstance implies 
that when the westwind circulation is strong 
aloft (and not too weak at the ground) the 
“neutral” and unstable waves move approxi- 
mately with the speed of the surface zonal 
current (a similar result has also been obtained 
by CHARNEY 1947). We find here a verifica- 
tion of a well known empirical rule frequently 
applied in forecasting: a disturbance, which 
does not move quickly “through” the upper 
air flow pattern (long wave), propagates 
eastward with approximately half the velocity 
of the westwind at ca. 5 km height. Now, the 
surface zonal (geostrophic) wind is in such 
cases indeed about one half of the middle 
tropospheric wind. As soon as the initial 
disturbance has assumed semi-mature dimen- 
sions, the results, derived from a theory of 
“small perturbations”, are not applicable. 

For weak surface westerlies or surface 
easterlies (decreasing. upwards) the theory is 
not in agreement with the observations. As 
far as the sense and speed of propagation 
is concerned it breaks completely down. 
However, it brings out a perhaps not un- 
characteristic other feature. It is seen that a 
given shear rate may be sub-critical (not 
involving instability) when the surface circula- 
tion is from west to east, but be super-critical 
when it is from east to west. For example at a 
shear 0.6 m/sec per kilometer elevation, a 
value of %o = + 10 m/sec would predict sta- 
bility whilst with the same shear rate at a 
moderately strong surface east wind (wo = 
= — 10 m/sec) we are well within the region 
of instability. The significant time-rate of 
growth, as defined above, is, in the latter case, 
58 hours. (See fig. 3.) 

The above results show how important for 
the breakdown of zonal motion the absolute 
circulation appears to be. At subpolar latitudes, 
the concurrence of surface easterlies with 
upper westerlies is a common feature of the 
general circulation, The horizontally uniform 
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Fig. 3. Curves for the velocity deficit u, —c (heavy lines) 
and for the characteristic time of amplitude growth 7; 
(thin lines), plotted as functions of the surface zonal wind 
Wo, for different values of vertical shear A. The circles 
at the lower dashed line indicate that the selected shear 
value is critical. — These curves are evaluated for pres- 
cribed values of latitude @ (= 45°), of circumferential 
wave number n(= 7), and of static stability S (= 1.3 - 10-4 
sec—?, normal value). 


zonal current systems we have treated here, 
constitute an integral part of an anticyclonic 
cold polar vortex superimposed by a cyclonic 
vortex. By computing the kinetic and po- 
tential energy of similar current systems, 
RossBy (1949) has demonstrated that an even 
barotropic anticyclonic vortex centered at 
the pole is in a state of dynamic inequilibrium. 
Whether the (cold) vortex as a whole is displac- 
ed equatorwards or whether the transport of 
the cold air towards lower latitudes is realised 
by a breaking up of the zonal motion, seems 
to be dependent upon the dimensional pro- 
perties as well as the shape of the cold dome 
centered at, or near, the pole. 

In order to clarify the matter of the de- 
pendency of the instability on the intensity of 
the absolute circulation, we should add that 
this is, in the present analysis, merely a con- 
sequence of including in the (perturbation) 
vorticity equation the horizontal solenoid 
term. Now, the existence of p, © solenoids 
in the (meridional) vertical plane is absolutely 
necessary to bring about the phenomenon of 
instability considered here, because the ver- 
tical shear (necessarily asscociated with the 
vertical solenoids) controls the inequalities, 
in the vertical, of vorticity transport: to this 
extent this solenoid field controls the defor- 
mation of the flow pattern, in conjunction 
with the effect of the earth’s rotation. On 
the other hand, the horizontal solenoid field 
appears immediately once the perturbation 
has set in, but it only modifies the nature 
of the deformation of the streamlines and 
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is not a necessary condition for amplitude 
growth. 

In conclusion a very brief summary of the 
results of a development of the above ele- 
mentary criteria may be given. Through in- 
cluding into the simplified model of a poly- 
tropic atmosphere the effect of the vertical tilt 
of the pressure wave — assuming in conformity 
with empirical evidence the temperature wave 
to have no such tilt —, one can deal with cer- 
tain types of advection (horizontal mass 
transport). These are associated with critical 
wave-lengths which are longer or shorter 
than those calculated on the basis of the ele- 
mentary criterion. For just the type of ad- 
vection which is most commonly observed 
during the initial stages of a deepening large- 
scale disturbance (see also NYBERG, 1949) the 
critical. wave-length is shorter than that 
evaluated from the elementary criterion. 
The stability curve is then in still better agree- 
ment with the one derived by CHARNEY (1947). 


There is perfect agreement between these 
curves when the height of the polytropic 
atmosphere is replaced by the height of an 
equivalent polytropic atmosphere equal to the 
mean height of the tropopause in the middle 
latitude westerlies. — The inclusion of the 
above mentioned tilt permits finally to evaluate 
the phase lag, between the pressure wave and 
the temperature wave, for which at prescribed 
values of the Richardson number the growth- 
rate of the disturbance has a maximum. This 
phase lag corresponds to an extremely high 
rate of the relevant advection. 


After completion of this investigation a 
paper by Dr. E. T. Eapy has been published in 
this journal (Tellus, 1, 3, p. 33), wherein 
the author assigns to the value of the 
Richardson number a fundamental rôle for 
the behaviour of cyclone waves and long 
waves. 
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Abstract 


J. G. CHARNEY and A. ELrASSEN (1949) have proposed a numerical forecast-method for 
the height profile of the soo mb-surface at a fixed latitude, and they have also calculated 
the constants necessary for 24 hour forecasts at 45° N. From those constants the author 
has derived a 48 hour formula, which has been tried on some simple and regular upper air 
pattern-types, such as sinusoidal profiles of different wave-length, and some types of re- 
gular but assymetric quasi-sinusoidal profiles, the purpose of these tests being to ascertain 
the nature of the circulation changes implied by the Charney-Eliassen method. 


Introduction 


The numerical forecast-method for the height 
profiles of the soo mb-surface, as proposed by 
CHARNEY and ELIASSEN (1949), is a first attempt to 
replace the aerological forecasting procedure of to- 
day, i.e. a mainly qualitative and always somewhat 
subjective argumentation between factors of differ- 
ent kinds, with a quite objective, numerical method. 
The practical possibilities, however, of this pioneer- 
ing way of forecasting are not quite clear from the 
two authors’ mainly theoretical paper. Since the 
meteorological forcasters need an idea of the 
practical qualities of the method, it is necessary 
either to test the method against a series of actual 
situations, or to apply it to some simplified standard 
upper-air patterns. Synoptic tests of the prediction 
method are at present carried out at the University 
of Stockholm. Thus, we will in this paper deal only 
with the application of the forecast-method to some 
schematic height profiles. Thereby we will get a 
general idea of how the method works, and what 
kind of results we may expect in different cases. 

To make the difference between the starting 
profile and the forecast profile more striking than 
it would be with 24 hour forecasts only, we prefer 
studying 48 hour forecasts. In addition to that, 


many of the routine upper-air forecasts of to-day 
are just 48 hour ones. For that reason we will at 
first derive a 48 hour modification of the Charney- 
Eliassen formula. 


A 48 hour forecast formula 


We may write the Charney-Eliassen formula for 
24 hour forecasts in the following manner 


Zag (ie Us) = 2 Bu Zee ch 10 0) (1) 


where z0 = height at t =o, 
Z2, — height at t =24 hours, 


U,, = zonal speed in degrees of longitude 
per 24 hours, 

B, = numerical constants, the 45°-values 
to be found in Table I, 

x = longitude, 

M 0) Euch 2, 


On the right side of (1) we now find for each x, 
a forecast value z,,, which is to be plotted U,, 
degrees east of the point x. That means, that we 
can do the calculation in two steps: the first step is 
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Table I. 


Numerical constants By, in the 24 hour forecast 
formula (1). 


N a? = 10 a? = 14 a = 18 
3 =< 0,070 — 0:05? 10057 
2 — 0.100 — 0.080 — 0.064 
— il — 0.506 == 0.433 — 0.382 

fe) 0.788 0.818 0.839 
I 0.305 0.314 0.316 
2 0.119 0.112 0.102 
3 0.177 0.149 0,123 
4 0.065 0.051 0.038 
5 0.095 0.065 0,047 
6 0.034 0.022 0.014 
7 0.049 0.030 0.016 
8 0.017 

9 0.023 


to compute the z, „„values for each x, as was U = 0; 
the second one is to displace the resulting profile 
just U,, degrees to the east. 

If we now insert the undisplaced values of z,, 
(as calculated for U = 0) in the formula again, we 
get 

2,8 (x) = I By 2,, (x + 102) 
n 
and displacing that profile 2 U,, to the east, we 
have the 48 hour forecast. But (1) gives, when 
U = 0; 


Zo, (x + Ion) = DB„2,(x + Ion + 10m) 
Mm 


Thus 


2,8 (x + 2 U,,) = ZB, By 2 (x + 10 [n + my) 


n,m 


or 


28 +20, )= > C2 (+ ton) (2) 


N 


where C, are numerical constants, easy to compute 
from the B,-values, and to be found in Table II. 


In this paper we have used the constants, which are 
valid when a? = 14.1 If we want a more simple and 


1 According to CHARNEY-ELIASSEN, we have at 45° N, 
that 


where U is the zonal wind in radians per day, and s is 
the number of stationary waves, encircling the globe. 
U is always of the order of 0.3. Then, for instance, 
SS 3) gives a= TA 


Table II. 


Numerical constants Cy in the 48 hour forecast 
formula (2). 


N as 10. | a> 74 ma 78 | Approx. 
— 6 0.005 0.003 0.001 
u Vs) 0.014 0,008 0.004 
dl 0.080 0.051 0.032 
— 323.008 — 0.014 — 0,014 
—r 7 0.056 0.025 0,014 
— I | — 0.876 — 0.770 — 0.688 —- 0.8 
[6] 0.265 0.364 0.438 0.4 
I 0.314 0.388 0.434 0.4 
2 0.075 0.138 0.170 0.1 
3 0.262 0.257 0.232 0.3 
4 0.102 0.125 Q.XE2 0.1 
5 0.186 0.148 0.116 0.2 
6 0.101 0.084 0.065 O.1 
7 05222 0.089 0.054 OT 
8 0.078 0.046 0.025 0.1 
9 0.082 0.020 0.012 
Io 0.049 0.015 0.008 
II 0.024 0.006 0.004 
12 O.021 0.005 0.002 
13 0.010 0.001 0.000 
14 0.008 0.001 0.000 
I5 0.004 
16 0.002 
17 | 0.000 
18 0.001 


less time-consuming formula than the complete 
equation (2), we may use the constants in the last 
column of Table II. That simplification is in fact a 
rather good approximation, as will be shown later 
on. Besides, that is clear from CHARNEY-ELIASSEN’s 
statement, that there is very little difference in the 
forecasts made with a? = 10, a? = 14 and a? = 18. 


Some general qualities of the numerical 
forecast method 


As was already shown, when we deduced the 
48 hour forecast formula, it is possible to split all 
changes due to the Charney-Eliassen method, in 
two components. One is the deformation of the profile, 
not at all depending on the zonal wind Ur. The 
other one is the displacement of the deformed profile 
That displacement depends only on U. 

Of course, the method is not able to change the 
wave-length of a regular periodic profile. 

If the deformation of the profile produces a 


1 As far as the coefficients By, and Cy, do not depend 
on U, which is approximately true, 
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change in amplitude, that change is proportional . 


to the amplitude value. 

A forecast calculated from a profile, which is the 
sum of superimposed profiles, is equal to the sum 
of the forecasts calculated from each one of the 


individual profiles. 


A study of sinusoidal profiles 


Applying the Charney-Eliassen formula to regu- 
lar, schematic profiles, we get a clear idea of what 
changes we may expect in different cases. Besides, 
it will then be much easier to discuss fundamentally, 
if the method is in agreement with experience, than 
by applying it to observed upper-air profiles only. 

For that reason, we have tried the 48 hour formula 
(2) on sinus curves with wave-lengths between 10 
and 180 degrees!. We are then only interested in 
the first or deforming step (U = 0), as the second 
step is nothing but a displacement due to the zonal 
speed. Fig. 1 shows an example, where L = 60 
degrees of longitude. The full line is the starting 
profile, and the broken line gives the 48 hour fore- 
cast. Obviously the deformation step has moved 
the curve about 16 degrees to the west. So if the 
phase velocity, i.e. the speed of propagation of 
ridges or troughs, equals C,, per day, we get 


or 


U2, 7224 


Now other wave-lengths give other values of 
(U,, — C,,), ie. the difference between the zonal 
velocity and the phase velocity. The continuous 
curve in Fig. 2 gives those differences as a function 
of the wave-length L. The dotted line demonstrates 
the well-known formula 


L2 
U—C= B 
4 70 


found by RossBy (1939). 


After studying how that formula fits empirical 
data, Namras & Crapp (1944) found, that it gave 


1 We have done that in spite of the fact, that all those 
wave-lengths are not possible on the earth; to be so, the 
wave-lengths must of course be a factor of 360 degrees. 


SS 


Fig. 1. Diagram showing a sinusoidal profile (full line) 

with the amplitude A and the wave-length L = 60 

degrees of longitude, and the corresponding 48 hour 

forecast (broken line), according to the Charney-Eliassen 

method when U = o. For other values of U, the forecast 

profile is to be moved 2 U;, degrees to the right (to 
the east). 


20 


120 150 


0 0 60 90 180° long. 


Fig. 2. The difference (U24—C24) between the zonal and 

the phase velocity in degrees of longitude per day as a 

function of the wave-length L. The full line shows the 

results of the Charney-Eliassen method. The dotted line 

represents the Rossby-formula, whereas the broken lines 

demonstrate the Namias-Clapp values for certain Un, 
and C24. 


too high values for large wave-lengths; so they 
proposed a corrected formula‘, namely 


pL? 
Creed, — = + 205 


The broken lines in Fig. 2 demonstrate the last- 
mentioned formula, when U,, equals 10 and 20 
degrees per day, and when C,, = 0. 

For short wave-lengths, the numerical forecast- 
method thus is in good. agreement with the Rossby- 
formula. For wave-lengths between 90 and 150 
degrees it just coincides with the Namias-Clapp 
formula for U,, = 15, which is probably a good 
average value of the zonal speed at 45° N. 

In the example of fig. r, where L = 60, we have 


1 It must be mentioned, that this formula was found 
by studying five-day mean maps. So we are not quite 
justified to state, that it is valid as well for the movements 
of individual troughs and ridges, 
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Fig. 3. Percental 48 hour amplitude-changes for wave- 
lengths between 10 and 180 degrees of longitude, accord- 
ing to the Charney-Eliassen formula. 
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Fig. 4a, b. Diagram showing two contrary types of 

askew quasi-sinusoidal profiles (full lines), and the corres- 

ponding 48 hour forecasts (broken lines), valid when 
U= 0. 


got a slight increase of the amplitude (about 1 %). 
Other wave-lengths will give other amplitude 
changes, however. Fig. 3 shows how large those 
changes are. There is mainly an increase in amplitude 
for wave-lengths below 60 degrees and beyond 120 
degrees; for wave-lengths between those values, 
however, we have a decrease of up to 5%. The 
irregularities of the curve as for wave-lengths less 


than 30 degrees probably depend on the 10-degrees 
calculation intervals being too large here, compared 
with the wave-lengths themselves. 


If we use the approximate constants of Table II 
(last column), we will get the same values of 
(U,, — C3) as shown in Fig. 2, but the amplitude 
changes will differ somewhat from those of Fig. 3. 


Quasi-sinusoidal profiles 


We have also tried the forecast-method on regular 
but askew profiles such as given by the continuous 
curves in fig. 4. Looking at the dotted forecast 
profiles, we see that we have here in both cases got 
the contrary type. Having for instance, as in fig. 4 a, 
a trough, where the meridional wind is strongest on 
the back side, we evidently will get the strongest 
gradient on the front side 48 hours later, and vice 
versa. It is also seen, that we have in fig. 4a gota 
deepening and sharpening of the troughs and a 
flattening of the ridges, but in fig. 4b a rise and 
sharpening of the ridges and a filling of the troughs. 

In fig. 5 a we have a series of ridges, sharp and 
flat alternately. Also here we obviously get a fore- 
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Fig. 5a, b. Diagram showing two contrary types of 

quasi-sinusoidal profiles (full lines) with long and short 

wave-lengths alternately, and the corresponding 48 hour 
forecasts (broken lines), valid when U = o, 
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cast, which looks similar to the reflected image 
(fig. 5 b) of the starting profile. Of course that is 
also true in the opposite case. In these two cases we 
do not get any remarkable changes in amplitude. 


Superimposed profiles 


In fig. 6 a, b and c we have given an example of 
two superimposed profiles. The sum of the two 
curves in 6a is given by the continuous curve in 
6b. The dotted curve in that figure is the forecast 
curve. It is equal to the sum of the two individual 
forecasts, shown in 6c. We now see from fig. 6 b, 
that if we have a double trough, the forecast-method 
will make the trough to the west deepen, and the 
eastern one fill. In the double ridge, on the other 
hand, the ridge to the west will rise, whereas the 
easternmost one will decrease extensively. 


Summary 


We have given some schematic examples of the 
consequences of the Charney-Eliassen method. It 
is seen from that, that we can get a trough deepening 
in three ways. Firstly certain wave-lengths always 
give deepening, i.e. an amplitude increase. Secondly, 
if we have a trough, where the meridional compo- 
nent of the wind is stronger on the back side than 
on the front side, there will also be a fall. Thirdly, 
a trough with a small ridge in it, will deepen in its 
western part. 

Similar rules may be stated as to the ridges. 


80 90 100 


O 10 20 30 40 50 6 10 


0 10 


Cc 


Fig. 6 a, b, c. Two sinusoidal profiles of different wave- 

length and amplitude as well as the result of adding them 

(full lines), and the corresponding 48 hour forecasts 
(broken lines), valid when U= o. 
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Abstract 


A solution of the vorticity equation is obtained for a wave pattern changing with 
time. It is shown that such a wave pattern is necessarily associated with ridge and trough 
lines which do not run in the south—north direction. 


1. The long wave pattern of the streamlines in 
the upper troposphere of middle latitudes was inves- 
tigated first by Rosssy (1939). Assuming a purely 
horizontal motion of an inviscid homogeneous 
incompressible atmosphere, he obtained for a 
perturbed zonal current of constant velocity U the 
equation 


2 2 
(eur) (+ SE) + 
ot dx] \ 0x? 0 Ox 

where w is the stream function of the perturbation 
motion, ß is the meridional rate of change of the 
Coriolis parameter, which is assumed to be constant. 
The x and y axes are taken in the east and north di- 
rections respectively. This equation has been discus- 
sed by Rosssy (1939), and later by Haurwitz (1940). 
In the wave models considered by them, the ridge 
and trough lines are oriented in the south—north 
direction. Recently MACHTA (1949) obtained a wave 
solution of equation (1), in which the ridge and 
trough lines are “tilted”, ic. not running along a 
meridian. In all these investigations, however, the 
models contemplated are such that the wave ampli- 
tude neither gains nor loses in size with time. In 
what follows, a more general solution of (1) is 
obtained which gives rise to a wave pattern having 
tilted ridge and trough lines and, moreover, under- 
going development with time. 


2. A solution of equation (1) is sought in the form: 
y =e*[@, (y) cos A(x —ct) + D, (y)sin A (x—ct)] (2) 


where 2x/2 is the wave-length, c the wave velocity 
and oa parameter describing the development of the 
wave pattern with time. Such a form of solution 
would lead to a wave having an amplified (or 
damped) amplitude and such that the ridge and 
trough lines deviate from the S—N direction. 

Substitution of (2) in (1), leads to a system of 
differential equations in ®, and ®,, the general solu- 
tion of which may be written in the form: 


®, (y) = (A cos uy + Bsin wy) e® + 
+ (Ccosuy + D sin wy) e—% (3) 
Da (y) = (A sin Ly B cos Ly) ew + 
+ (D cos wy — C sin py) e—*” (4) 
where A, B, C and D are arbitrary constants, with 


« 2am B 
(22 D u2— @2)2 +4 a2 pe 


(5) 
and 


B (22 + y — 2) 
(22 + u? — a2)? = 4.02 pu? 


U—c= 


(6) 


A SOLUTION OF THE VORTICITY EQUATION $9 


It can be seen from (2), (3) and (4) that the per- 
turbation motion is given by a linear combination 
of the functions 


cos 


et Ak cd) — ny}, 
cos 

— t f \ 

ew +e fits (x — ct) a uy } 


and it is interesting to note that in the wave pattern, 
corresponding to any of these solutions, the ridge 
and trough lines are “tilted” straight lines whenever 
o +0, Le. x +0. This means that for any solution 
involving a wave pattern changing with time, the ridge 
and trough lines are necessarily “tilted”. 

The models considered by RossBy (1939), HAUR- 
wırz (1940) and MAcHTA (1949), can be ob- 
tained from the general formulae (2)—(6) by 
choosing the values of the parameters and con- 
stants in the following manner. 


RossBy: 0e AC On) = B= 6, 
Haurwitz: ¢=0 u$0oA=C=oB+D=o, 
NA Cara Oy seo A = C = -D-= 0 BE 0. 


The solution of Macura (1949) is in fact a 
linear combination of two solutions of the Haur- 
WITZ (1940) type. 

3. In this section the relation between the tilt 
of the ridge and trough lines and the amplification 
(or decay) of the waves is discussed. A ridge or a 
trough line will be called positively tilted when the 
positive angle it makes with the x-axis is less than a 
right angle, otherwise it is described as negatively 
tilted. Also the meridional Aux of west momentum 
and kinetic energy will be examined. These are 
ouv and U ouv (van MIEGHEM 1949), respectiv- 
ely, where o is the density, u and v are respectively 
the z and y-components of the perturbation ve- 
locity, and the bar to denote the average value. 
In the case under consideration o is constant, and 
the average value is to be taken with respect to a 
wave length. 

The streamlines for the total motion, obtained 
by the superimposition of the perturbation and the 
mean motion are given by y, = constant, where 


1 = — Uy + y (7) 


y being given by (2), (3) and (4). To have an idea 
of the orientation of the ridge and trough lines, it 
will be sufficient to find the loci of the points 


0 Je 
where v = 0, i.e. where Sn = o. Differentiating 
x 


the expression so obtained we find that the slope, 
dy/dx, of a ridge or a trough line is given by: 


dx 29, — 9,6, 
dy D? 


A cos? A(x —ct) (8) 


We note also that 


— A 
FFT. ed, PB, — D,'®,] (9) 


which shows that uv is of the same sign as dy/dx. 

It is easy to verify, from (3), (4) and (8), that if 
A =B=0, or C= D =o, the ridge and trough 
lines are straight. It follows that, in this case, uv 
will be of constant sign; and it can be seen that 
west momentum and kinetic energy are trans- 
ported northward for positively, southward for 
negatively, tilted (straight) ridge and trough lines; 
a result pointed out by MACHTA (1949) when o = o. 
The amplification, or decay, of the waves depends, 
by (5), on the sign of (a). It thus becomes evident 
that with straight ridge and trough lines an amplified 
wave pattern is associated with northward increase 
of amplitude if it is positively tilted, and a southward 
increase if it is negatively tilted. 

For a decaying wave, on the other hand, the 
amplitude increases to the south in the case of a 
positive tilt, and to the north in case of a negative 
tilt. The meridional flux of kinetic energy and west 
momentum is towards the region of increasing or 
diminishing amplitude, according as the wave 
amplitude is gaining or losing in size with time. 

We examine now the more general case A? + 
+ B? + o and C? + D? + o. By (8), (3) and (4), 
we find after some reduction that 


dx usinh (2«y + 2k) — «x sin (2uy + 6) 
= (to) 
dy cosh (2«y + 2k) — cos (2 my + 0) 


A? + B2 
where exp ch) = ve 
4 oe, ue BD —AC 
a = (OB AWE 


The denominator of the fraction on the right hand 
side of (10) is always positive and when it vanishes, 
it can be seen that the limit of the fraction is zero 
Thus the sign of dx/dy which indicates the orienta- 
tion of a ridge or a trough line at any point, is 
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determined by the sign of the numerator. Two cases 
are considered: 


(i) Amplified Waves & > 0, (au > ©). 


It is seen from (10) that when y is sufficiently 
large and positive, then dx/dy is positive, and when 
it is sufficiently large and negative dx/dy is negative. 
This means that sufficiently far to the north the 
ridge and trough lines are positively tilted, while 
sufficiently far to the south they are negatively 


tilted. Since the sign of uv is the same as that of 


dx/dy, it follows that uv is positive far enough to 
the north and negative far enough to the south; 
and hence there is an outflow of west momentum 
and kinetic energy from the central part of the 
stream for waves which are gaining in size with 
time. 


(ii) Damped Waves o <o, (au <0). 


Similarly as in (i), it can be seen that the ridge 
and through lines are negatively oriented in the 
northern part of the stream, and positively oriented 


in the southern part. Thus wv is negative far enough 
to the north and positive far enough to the south, 
corresponding to an inflow of kinetic energy and 
west momentum in the central part of the flow. 


It thus appears that a wave pattern changing with 
time is always associated with a tilt of the ridge and trough 
lines. Furthermore, when the configuration of 
these lines is as described in (i) or (ii) above, amplifi- 
cation or decay of the perturbation motion takes 
place respectively, accompanied by an outflow or 
inflow of west momentum and kinetic energy in 
the central part of the stream. This last remark is in 
agreement with a result obtained by VAN MIEGHEM 
(1949). He in fact has shown, from momentum 
and energy considerations, that with the ridges and 
troughs oriented SE—NW to the north of the 
zonal current, and SW—NE to the south of the 
stream, there is accumulation of west momentum 
and kinetic energy which would increase the strength 
of the zonal westerly current. From the above, it 
is seen that this corresponds to the case of a damped 
perturbation, which agrees with a result obtained 
by Kuo (1949). 


Acknowledgement. 


The author wishes to express his gratitude for 
Prof. VAN MIEGHEM for many suggestions and 
valuable discussions. 

Institut Royal Météorologique de Belgique. 
Uccle, July 1949. 


RE RER EINIGES 


Havurwitz, B., 1940: The motion of atmospheric 
disturbances. J. Mar. Res., 3, 1, pp. 35—50. 

Kuo, Hsıao-LAn, 1949: Dynamic instability of two 
dimensional non divergent flow. J. Meteor., 6, 2, 
pp. 105—122. 

MACHTA, L., 1949: Dynamic characteristics of a tilted 
trough model. J. Meteor., 6, 4, pp. 261—265. 

MIECHEM, J. VAN, 1949: Production et redistribution 


de la quantité de mouvement et de l'énergie ciné- 
tique dans l’atmosphère. Application à la circulation 
atmosphérique générale. J. Sc. de la Météor., 1, 2, 
PP 307% 

RossBy, C. G., 1939: Relation between varations in 
the intensity of the zonal circulation of the atmo- 
sphere and the displacements of the semi-permanent 
centres of action. J. Mar. Res., 2, 1, pp. 38-55. 


REG ITES 


The Ionospheric and Radio Wave Propagation Observatory 
ac Kirima,,. 67°50’) N, -20°14.5’ U E 


During the last war the Telecommunications 
Laboratory of the Chalmers University of Technology, 
Gothenburg, started a small ionospheric observatory 
at Askim (s7°14’ N, ı1°59’ E) six miles south of 
the city. Fixed frequency recordings were begun 
in 1941 and a year later, in 1942, a semi auto- 
matic multi-frequency recorder was put in opera- 
tion. Already during the first years of operation 
of the Askim observatory, which unfortunately 
had to be very irregular during the war, a de- 
finite need was felt for a recording sister-observ- 
atory in the very north of Sweden where the 
auroral disturbances are much more frequent than 
at Gothenburg. 

On account of the often very transient con- 
ditions of the “auroral” ionosphere the Research 
Laboratory of Electronics, which now operates the 
Askim observatory, decided to develop a quick- 
sweep high power panoramic recorder for the 
proposed new observatory. It was felt that ex- 
tremely valuable information could be obtained 
by such an instrument with a sweep speed of at 
least 2/3 Mc/s per second.! Preliminary construc- 
tion work was begun in 1945 but was delayed by 
the lack of suitable tubes which were obtained 
from the United States first in 1946. The new 
recording equipment was ready early in 1948 and 
it was decided to move it to a temporary loca- 
tion at Abisko. At the suggestion of Professor 
Rorr M. Sievert, however, this plan was aban- 


doned and the equipment was moved to tempo- — 


rary buildings near Kiruna, the iron ore city of 
the north. The buildings were erected on govern- 
1 The development of such an instrument was dis- 
cussed already in 1939 by WM. O. Reep and O. E. EL. 
Rypseck, both at that time working at HARVARD UNI- 
VERSITY’S CRUFT LABORATORY. 
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ment property where it is hoped that a Swedish 
geophysical research centre will be organized in 
the near future. The temporary buildings, which 
were almost completed by midsummer 1948, 
were placed at the disposal of the Research Labora- 
tory of Electronics through the good offices of the 
Interim Board of the Research Stations in northern 
Sweden of the Royal Swedish Academy of Sciences. It 
is in no small measure due to the positive support 
by Professor Sievert, Chairman of the Board, 
that the Kiruna Ionospheric Observatory was cre- 
ated so rapidly. Regular operation at the new ob- 
servatory was begun on Oct. 1, 1948, and has 
continued practically without interruption from 
that date. 


+ 


The new observatory is casily accessible by 
electric railroads as shown in fig. ı and is situated 
at the south-eastern edge of the zone of maximum 
auroral intensity. The Kiruna observations there- 
fore must serve as very valuable complements to 
those made at Tromsö. The geomagnetic co-latitude 
of the Kiruna observatory is about 24°21’ with a 
magnetic dip of 77°03’ and a total magnetic force 
at ground level of 0.515 CGS. 

The international chain of ionospheric observ- 
atories is steadily being increased as shown in 
fig. 2 which depicts the location of all ionospheric 
observatories known to us to carry out regular 
or irregular observations. In the auroral zone the 
Fairbanks, Churchill, Tromsö and Kiruna observ- 
atories probably will play an important röle. Of 
these observatories only the Kiruna observatory is 
equipped with high power panoramic recording 
equipment (as far as is known to us at the time 
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Fig. 1. The location of the Ionospheric and Radio Wave 
Propagation Observatory at Kiruna. 


(© = approximate geomagnetic co-latitude.) 


of writing) but it is believed that the Churchill 
observatory soon will get equipment of that type. 

A few words should be said in this connexion 
about the research activity and observational pro- 
gramme of the Kiruna observatory. As far as regu- 
lar observations are concerned the results of these 
will be distributed monthly to interested observ- 
atories all over the world. The material is prepared 
in the standard fashion used by the Central Radio 
Propagation Laboratory, Washington, D. C. This 
material will also appear regularly in the Reports 
of the Research Laboratory of Electronics, Chalmers 
University of Technology, Gothenburg, where they 
can be obtained upon request either from the 
Laboratory or from the Chief Librarian of the 
University. 

The ionospheric research at present concerns 
such important and for the practical radio com- 
munication vital questions as the correlation be- 
tween sporadic E-ionization, Polar Black Out and 
Magnetic Bays. Further, if possible, the statistical 


correlation between these phenomena and low- 
latitude fade-out effects caused by solar ultra violet 
light. Auroral radio echo forms and the connec- 
tions between them and the lower ionosphere 
auroral current system are studied regularly. Spe- 
cial attention is being paid to abnormal low height 
reflections, sporadic echoes already having been 
recorded at such low heights as 70 km (special 
report in progress of publication). From a radio 
wave propagation point of view there is especially 
to mention the recording and study of geomag- 
netic triple splitting of the waves reflected from 
the northern ionosphere. The irregular and some- 
times rapid motions of the upper ionosphere, the 
F2-region, is a further object of repeated study. 
Very interesting motion pictures of such motions 
have already been demonstrated by members of 
the Staff of the Research Laboratory of Electronics’. 

Temporary magnetic recordings are carried out 
at the Kiruna ionospheric observatory as an aid to 
the local observer but it is hoped that the Royal 
Hydrographic Office soon will establish high- 
quality permanent recordings in the neighbour- 
hood. 

The observatory is further equipped with an 
auroral spectrograph, to which a photo-electric 
attachment will later be added. In order to study 
the long distance scatter of radio waves from the 
auroral arcs and draperies a fixed frequency record- 
ing equipment with about soo kW pulse sending 
power will be installed this winter. This equip- 
ment is of the same type as that used for the re- 
cording of meteor trails at the Askim observatory 
and will be described in a forthcoming publication 
from the Research Laboratory of Electronics. The 
description of the panoramic equipment will ap- 
pear very shortly. 

It is further planned to open crystal controlled 
long base recordings between Kiruna and the 
mother-observatory in Gothenburg but these plans 
are not likely to be realized this year. 

+ x 
x 

Fig. 3 shows the general lay-out of the two 
main buildings of the observatory. The interested 
reader will note that the observatory has a guest 
room. Visiting scientists are always welcome as 
guests for shorter or longer periods. 

Fig. 4 depicts a corner of the recording room 
showing the 6 kW panoramic recorder which 
sweeps from 1 to 20 Mc/s in 30 seconds or less. 


1 Nature, Vol. 164, Sep. 24, p. 512, 1949. 


D FAIRBANKS 


dPORTAGE LEE PRAIRIE 
ST. JOHNS 


BOSTON Sy @ ALMA ATA 
SAN FRANCISCO) WS WASHINGTON i v ———— 


LANCHOWe  NANKI 


CV 
ITESANDS | q Week 
| | BATON ROUGE ELHI j 0%), @YAMAKAWA 
Vs ) S 1 ene WUCHANG ya 
SE N JUAN, PUERTO = \ ) 
| | N are TRINIDAD 


SUVA, FIJI Is. 
© RAROTONGA T. 


1 
RAOUL, KERMADEC |. @ JOHANNESBURG 


a TOWN) 


FALKLAND Is. 


180 160 Mo 120 100 80 60 40 20 0 20 40 60 80 100 120 4O 160 _ 180 
LONGITUDE WEST LONGITUDE EAST 


Fig. 2. World map of known Ionospheric Recording Stations with regular or irregular activity. 
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Fig. 3. Plan of the two main buildings of Ionospceric and 
Radio Wave Propagation Observatory at Kiruna. 
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Fig. 4. A corner of the recording room showing the 
20 kW panoramic recorder with 16 mm automatic 
film camera. 


Bios: 


General view of the main observatory 
buildings. 


To the right automatic recording 16 mm film 
camera in front of 12GP7 scope tube. 

Fig. 5 is a general view of the main buildings 
with the central mast carrying the aperiodic an- 
tennas for vertical broad band radiation. Fig. 6 
finally is a picture of the observatory with the 
surrounding forests and mountain ranges which 


Fig. 6. The observatory is surrounded by forests and 
long mountain ranges. 


certainly must attract the attention of every nature 
loving scientist. 


+ 


We hope that this new observatory will be a 
valuable link in the international chain of iono- 
spheric observatories and we hope that it will con- 
tribute both to our knowledge of the regular 
ionospheric phenomena and of the complicated and 
fascinating relations between solar, geomagnetic, 
auroral and ionospheric phenomena. 

Our thanks are especially due to the National 
Research Councils for the Technical and Natural 
Sciences which have generously facilitated the 
development of the observatory equipment and 
the erection of the observatory itself. Our thanks 
are further due, and are cordially extended, to 
Professor RoLF M. SıEvErT for his encouraging 
support and to the Wallenberg Foundation for gener- 
al support of our research activity. 


Olof E. H. Rydbeck 


Chalmers University of Technology, 
Gothenburg. 


New publications 


Glaciers and Climate 


Among those scientists who in recent years have 
been particularly concerned with the progress of 
the glaciological science, Professor Hans W:son 
Ahlmann of the University of Stockholm has been 
one of the leaders and pioneers. On the fourteenth 
of November 1949 Ahlmann entered his sixty- 


first year. At this anniversary his colleagues in diffe- 
rent countries, his many friends and grateful pupils 
have wished to honour him and their hommage 
has taken the form of a beautiful volume including 
geophysical and geomorphological articles dealing 
primarily with glaciers and their relation to climate, 
by thirtyfour scientists representing different parts of 
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the world. The »contributions are dedicated to 
Hans W:son Ahlmann as a tribute from the Swe- 
dish Society for Antropology and Geography and 
from a large number of his collegues, friends and 
pupils». The volume was edited by F. Hjulstrém, 
University of Uppsala, L. Kirwan, Royal Geo- 
graphical Society, London, C. M:son Mannerfelt, 
University of Stockholm, C.-G. Rossby, Univer- 
sity of Stockholm and University of Chicago, A. 
Wickman, Litographic Office of the General Staff, 
Sweden, and A. Angstrém, Swedish Meteorological 
and Hydrological Institute, Mannerfelt being the 
chief-editor. It is issued as the first two numbers 
of the »Geografiska Annaler» of 1949 under the 
heading »Glaciers and Climate» and distributed by 
Generalstabens Litografiska Anstalt, Stockholm, 
Sweden. 

The articles can be seperated in different groups 
according to their scientifique content. The volume 
starts with a short note »Glaciers and climatology: 
Hans W:son Ahlmann’s contribution» by L. Kir- 
wan, C. M:son Mannerfelt, C.-G. Rossby and 
V. Schytt, which concludes with some beautiful 
pictures from the glaciated regions Ahlmann 
has visited during different expeditions. Among 
the articles then following there are not less then 
twelve devoted to the matter of climatic varia- 
tions and related problems. Some of them deal with 
purely climatic studies as the following: A. 
Defant, Austria: Konvektion und Eisbereitschaft in 
Polaren Schelfmeeren. — J. Eythorsson, Iceland: 
Temperature variations in Iceland. — B.. Helland- 
Hansen, Norway: Remarks on some variations in 
atmosphere and sea. — H. E. Landsberg, USA: 
Climatic trends in the series of temperature obser- 
vations at New Haven, Connecticut. — G. H. Lilje- 
quist, Sweden: On fluctuations of the summer tem- 
perature in Sweden. — S. Petterssen, Norway: 
Changes in the general circulation associated with 
recent climatic variation. — H. C. Willett, USA: 
Solar variability as a factor in the fluctuations of 
climate during geological time. — A. Ângstrôm, 
Sweden: Atmospheric circulation, climatic varia- 
tions and continentality of climate. They all give 
contributions to the discussion on the recent cli- 
matic variation and its showing up in changes of 
meteorological elements and the general circula- 
tion of the atmosphere. As far as Willett’s article 
is concerned also suggestions for an explanation 
of the causes of climatic changes are given. 

Closely related to these articles are some others 
which primarily deal with glaciological and Quar- 
ternary geological evidences for climatic variations: 


C. E. P. Brooks, Great Britain: Post glacial cli- 
matic changes in the light of recent glaciological 
research. — I. Hustich, Finland: On the correlation 
between growth and the recent climatic fluctua- 
tion. — J. Leighly, USA: On continentality and 
glaciation. — G. Manley, Great Britain: The snow- 
line in Britain. — E. Nilsson, Sweden: The pluvials 
of East Africa. — C. C. Wallén, Sweden: The 
shrinkage of the Karsa Glacier and its probable 
meteorological causes. 

Ahlmann’s scientifique contributions have to a 
large extent been devoted to purely glaciological 
problems. It is therefore reasonable that there should 
also be articles concerned with snow and gla- 
ciers in a more limited meaning: G. J. Klein, 
Canada: Canadian survey of physical characteristics 
of snow-covers. — W. V. Lewis, Great Britain: 
Glacial movement by rotational slipping. — E. 
de Martonne, France: La dyssymetrie des glacia- 
tions regionales et ses facteures. — V. Schytt, Swe- 
den: Re-freezing of the meltwater on the surface 
of glacier ice. — G. Seligman, Great Britain: Rese- 
arch on glacier flow. — C. Troll and K. Wien, 
Germany: Der Lewisgletscher am Mount Kenya. — 
W. Werenskiold, Norway: Glacier measurements 
in the Jotunheim. 

Other branches of physical geography are also 
represented in the publication through some ar- 
ticles. Two of them are particularly dealing with 
Ice-age problems: R. Foster Flint, USA: Pleistocene 
drainage diversions in South Dakota and C. M:son 
Mannerfelt, Sweden: Marginal drainage channels as 
indicators of the gradients of quarternary ice caps. 
Two articles discuss other geological tools in the 
study of climatic changes than the glaciers: S. Thora- 
rinsson, Iceland: Some tephrochronological contri- 
butions to the vulcanology and glaciology of Ice- 
land and A. Noe-Nygaard, Denmark: Samples of 
volcanic rocks from the sea bottom between the 
Faroes and Iceland. Articles which can be said to 
represent the field between Quarternary geology 
and geomorphology are: F. Hjulstrém, Sweden: 
Climatic changes and river patterns. — H. Baulig, 
France: Causalité et finalité en géomorphologie. — 
R. Blanchard, France: Le probleme du Temisca- 
mingue. — G. Lundqvist, Sweden: The orientation 
of the block material in certain species of flow earth. 
— A. Schou, Denmark: Danish coastal cliffs in 
glacial deposits. — J. Sölch: Austria: Ueber die 
Schwemmkegel der Alpen. Of a quite general 
geophysical character is finally: H. U. Sverdrup, 
Norway: Theoretical tools in geophysics. 

C. C. Wallen. 
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Arkiv för geofysik 

Since the beginning of this century the Royal 
Swedish Academy of Science has published »Arkiv 
för matematik, astronomi och fysik», which has 
appeared in 36 volumes during these almost fifty 
years. Because of the increasing number of contri- 
butions during the last years the Academy recently 
decided to divide this series, covering such a large 
field of the natural sciences, into three different 
series, VIZ: 

Arkiv för matematik, 

Arkiv för astronomi, 

Arkiv för fysik 
and to add a new series devoted to geophysical 
contributions called Arkiv for geofysik. 

As before printing and publication is handled by 
»Almqvist och Wicksells Boktryckeri», Uppsala and 


the separate issues will appear at irregular intervals. 


Six issues will form one volume. With the new series 
the format has been increased to 175—250 mm, 
being more convenient, in particular from a geophy- 
sical point of view: 

Vol. I No. 1 of »Arkiv för geofysik» will appear 
in December this year. This issue will contain the 
following papers: 

1) NORINDER, HARALD & SIKSNA, REINHARDS: The 
measurements of the: variation characteristics of 
the density of small-ions. 

2) NORINDER, Haratp & Sarka, Oscars: The 
effect of geological discontinuities on artificial 
lightning discharges as studied in a model test 
arrangement. 

3) Lönnovist, Orov: A general method and a 
simplified formula for calculating effective radiation. 

4) Raas, Lars & Ropsxjer, Nits: A study of the 
accuracy of measurements of the Vaisala radiosonde. 


B. Bolin. 


BOORSREVIEWS 


Franz Baur: Einführung in die Grosswetter- 
kunde, Wiesbaden 1948, 165. pp. Dieterich’sche 
Verlagsbuchhandlung. 


This small volume by Franz Baur is as interesting 
for the author’s philosophical approach to the 
problem of the longer-period variation of the 
weather and its forecasting as for the weather 
facts which are presented. In his introduction Baur 
defines his basic concepts and justifies the recogni- 
tion of the Grosswetterkunde as an important branch, 
in its own right, of the science of meteorology. 

By the term, Grosswetter, Baur refers to the large- 
scale weather pattern over one sector, predomi- 
nantly continental or maritime, of either hemisphere. 
He divides the northern and the southern hemi- 
spheres each into five such sectors, and defines the 
Welt-Grosswetter in terms of the ten individual sector 
weather patterns. This division of the general circu- 
lation pattern into a number of individual sectors in 
either hemisphere is useful as a basis for analog classi- 
fication of the circulation patterns, but Baur does 
not combine the sector patterns into any single in- 
tegrated hemisphere pattern or hemisphere circula- 
tion indices. He deals with the significant features of 
the regional (usually the European sector from 30° W 
to 45° E) weather patterns which are persistent 
(gleichbleibend) over a period of days (Grosswetterlage), 
the Grosswetterlage largely determining the local 
day-to-day sequence of weather in each locality. 

Baur recognizes that the first (and minor) problem 
of Grosswetter research is to establish the fact that 
there actually exists a real Grosswetter in the sense 
that the probability of occurrence of specific weather 
phenomena over long periods of time is determined 
by changing external weather-conditioning causes, 
rather than merely by the superposition of persist- 
ence and internal “chance” occurrences on the 
climatologically normal pattern. An example of an 
internal “chance” occurrence is the deposit of the 
precipitation of one heavy storm just on the snow 
side, rather than on the rain side of the freezing 
temperature, to furnish an extensive snow cover. 
Baur admits that without the presence of external 
(i.c., external to the weather pattern itself, not ne- 


cessarily extraterrestrial) superposed weather con- 
trols, Grosswetter research reduces to a study of the 
comparative morphology of weather patterns which 
has only limited forecasting significance. But once 
the existence of a real Grosswetter is established, then 
the principal problem of Grosswetter research is to 
determine the terrestrial or extra-terrestrial influen- 
ces which condition the time-and-space variations of 
the general circulation pattern. It is primarily this 
problem to which Baur addresses himself. 


In the appendix to this volume, as a subject which 
is too technical for the average lay reader, Baur 
proves to his own satisfaction, by a statistical tech- 
nique which is unfamiliar to this reviewer and too 
elaborate for complete reproduction in this brief re- 
view, that a real Grosswetter does exist. For this pur- 
pose, he selects arbitrarily some nineteen individual 
measurements of weather elements, the occurrence 
or non-occurrence of which can be taken readily each 
day from the weather map. For each element, for 
each of the past so years, he forms a series of the 
daily occurrence or non-occurrence of the selected 
value during a fifty-day period. For eleven of these 
weather elements, the so yearly series run from 
July 1 through August 19, and for eight of them the 
series runs from January 1 through February 19. For 
each of these nineteen sets of so annual fifty-day 
series of values, Baur computes the average “diver- 
gence coefficient”, which is the ratio of the empirical 
scattering to the “Markoff” scattering (the “Mar- 
koff” scattering being the Bernoullian scattering 
corrected to eliminate the effect of serial correlation). 
The so series averages of each of these nineteen di- 
vergence coefficients ranges from 1.03 to 1.40, 
average 1.24. A value of D = 100 is to be expected 
if these series represent only the operation of persist- 
ence and chance occurrences. According to Baur, 
this result lies far beyond the limits of any chance 
occurrence, and proves that there exist superposed 
changeable complexes of conditions by reason of 
which the probability of the occurrence of charac- 
teristic weather features over longer periods of 
time is different from year to year, i.e., there exists a 
real Grosswetter in the sense originally defined. 
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Accepting the existence of a real Grosswetter, 
Baur procedes to investigate its nature and physical 
controls. He discusses first the Grosswetter patterns 
of the European sector of the northern hemisphere, 
indicating how the other weather elements are re- 
lated to the circulation (pressure pattern), both 
current and past. Baur emphasizes in particular the 
relationship of the pressure pattern at the tropopause 
level to the temperature field and sea level pressure 
pattern, and the dominant role of upper level steering 
in the movement of isallobaric centers at sea level. 
He completes this section with a rather elaborate 
classification of sea level circulation patterns in the 
European sector, and some interpretation of the 
climatology or seasonal distribution of these circula- 
tion patterns. There is little in this section of the 
book which will be new or surprising to any meteo- 
rologist who has worked with daily and mean cir- 
culation patterns over large areas of the earth’s 
surface. 

In the second section, Baur extends his discussion 
of Grosswetter patterns to the northern hemisphere. 
He discusses the heat balance of the atmosphere and 
the necessary meridional transport of heat in the 
usual terms (after Baur and Putuipps). He derives the 
form of the general circulation in the widely accep- 
ted pattern of a thermally direct circulation which is 
impressed into zonal cells by the fact of the earth’s 
rotation. The pattern of meridional cells is deter- 
mined in the mean primarily by the seasonal mon- 
soon effects of continents and oceans. The necessary 
meridional heat transport is effected in part by the 
latitudinal austausch of these large semipermanent 
seasonal cells, and in part by that of the irregularly 
migratory pressure centers of middle latitudes. The 
irregular (non-seasonal) fluctuations of the general 
circulation pattern Baur sees essentially as fluctua- 
tions about a normal seasonal heat balance and heat 
transport, any disturbance or unbalance in one direc- 
tion being compensated by a swing off balance in 
the other direction. He offers no comment on nor 
interpretation of the mechanics of the widely varying 
cellular patterns of the general circulation which 
occur from week to week. His following classifica- 
tion of Welt-Grosswetterlagen consists merely of 
setting side by side the five sector pattern designators 
for the northern hemisphere above the correspond- 
ing five designators for the southern hemisphere, 
with the latitude of the heat equator being indicated 
between the two hemispheres. This classification 
of world ‘weather patterns, which. is capable of 
almost infinite variation, is sadly lacking in any uni- 
fying mechanical or dynamic principle, or in any 


unified expression in terms of simple indices of the 
state of the general circulation as a whole. 

With respect to the determination of the terres- 
trial and extra-terrestrial influences that affect his 
Grosswetter patterns, Baur does not have much evi- 
dence to offer which is essentially new or finally con- 
clusive. Among the terrestrial influences which he 
discusses briefly are volcanic outbreaks (Humphrey’s 
theory of the reduction, of insolational heating by 
volcanic ash), the fluctuation of the temperature and 
strength of the ocean currents, ice conditions in 
the Polar seas, and minor oscillations of the inclina- 
tion of the axis of rotation of the earth. The amount 
of this latter variation is utterly small and must be 
insignificant in effect. With respect to the signifi- 
cance of the other three factors, there is statistical 
evidence pro and con, but Baur has no new evi- 
dence to add to previous investigations which have 
led to inconclusive or negative results. Particularly 
with regard to ocean currents and polar ice, the 
statistical and synoptic evidence favors the point of 
view that the ice and current conditions are predo- 
minantly, if not exclusively, results rather than 
causes of the general circulation patterns. Ocean 
temperature or ice conditions at times may be indi- 
cators of future weather patterns, but not in them- 
selves be the physical causes. 

In the fourth section of his book, Baur reviews 
briefly the question of periodicity in the Gross- 
wetter. The daily and seasonal periods find almost an 
infinite variety of expression in weather pheno- 
mena, the seasonal periodicity finding expression in 
the monthly and seasonal frequency distribution of 
Baur’s Grosswetter types. He continues with a dis- 
cussion of the periodic analysis of pressure waves or 
rhythmic variation of pressure, including the compu- 
tation of correlation coefficients for pressure values 
paired on opposite sides of symmetry points as 
defined by Weickmann. Baur says that two years’ 
work spent on statistical analysis of this type at the 
Long-range Weather Forecasting Institute failed to 
produce enough results of forecasting significance to 
justify its continuance. In an interesting section in the 
Appendix he shows by actual experiment that by 
the drawing of numbers in a manner which imposed 
a moderate degree of serial correlation more or 
less corresponding to that in weather data, he ob- 
tains series which contain periodic variations and 
even symmetry points and symmetry correlations 
that are almost identical with those found in series 
of daily pressure observations. In other words, Baur 
attributes the great variety of periodicity and sym- 
metry patterns which may be derived from pressure 
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series to be primarily the incidental result of serial 
correlation in the data, and as such to have no 
greater forecasting significance than the persistence 
itself. With this conclusion this reviewer is in hearty 
agreement! 

In his discussion of the extra-terrestrial influences 
on the Grosswetter patterns (Section V), Baur con- 
siders the effects of the moon and the planets, and 
of solar radiation and sunspots. Except for the mi- 
nute tidal effect ofthe moon on atmospheric pressure, 
he concludes that empirical experience (synoptic 
and statistical) has been able to detect no definite 
evidence of any effect either of the phase of the 
moon or of the planets on terrestrial weather. Most 
meteorologists today are in agreement with this 
conclusion. 

In his discussion of solar variability Baur assumes 
that the final revised monthly mean values of the 
solar constant published by the Smithsonian Insti- 
tution actually express true solar variability, an 
assumption which the statistical work of PAGE! opens 
to serious question. Using monthly mean solar 
constant values from 1919—1937, Baur produces 
some evidence for a double solar constant period 
within the sunspot cycle, with solar constant minima 
occurring at sunspot maximum and at sunspot 
minimum. He accounts for the solar constant mini- 
mum at sunspot maximum by the increased absorp- 
tion of the outgoing solar radiation in the higher 
solar atmosphere at time of maximum, because of 
the increased amount of absorbing material which 
is driven into the higher atmosphere at time of 
maximum solar disturbance. Baur also produces 
some impressive statistical-synoptic evidence (by 
correlation), for the same 1919—1937 period, of in- 
creased zonal westerly circulation over northern 
Europe in midwinter (January-February) with in- 
creased solar constant, but in view of the shortness 
of the period (two sunspot cycles) and the serious 
doubt concerning the reality of the solar constant 
variation, this reviewer would be very hesitant to 
accept these results at their face value. 

Baur agrees with most other investigators when 
he says that the eleven-year sunspot cycle is reflected 
beyond much doubt in meteorological data only 
in the tropics, but even there his evidence is not 
really conclusive. He extends some of the older 
correlations of Körpen and C. E. P. Brooxs to 
more recent data with appreciable deterioration 
of the earlier correlations. He finds statistical evi- 
dence in the analysis of 200 years of Central Euro- 


1 Monthly Weather Review, Supplement 39, pp. 
TES- 120. : 


pean temperature data and 100 years of precipita- 
tion data for a double weather cycle, within the 
sunspot cycle, corresponding to his double cycle of 
the solar constant, and likewise for a corresponding 
cycle of the interseasonal monsoonal pressure con- 
trast in Central Asia. However, the amplitudes and 
regularity of these long-period fluctuations are not 
impressive. Certainly further statistical verification 
of their significance than that given here by Baur 
would be desirable. By and large, Baur does not 
add much to the previous inconclusive evidence as 
to the degree of the dependence of world weather 
patterns on irregular solar activity. He recognizes 
the present unsatisfactory state of this evidence by 
his rather cautious general conclusion on the topic, 
with which this reviewer is in full accord, that at 
least certain portions of the solar radiation are sub- 
ject to fluctuations to which the changes of the 
Grosswetter are definitely related, but that the many 
problems which are posed by this relationship can 
be solved only by close co-ordination of further 
research in atmospheric and solar physics. 

In the final section of his book Baur discusses 
briefly his methods of extended (three to five days) 
and long range (ten-day, monthly and seasonal) 
weather forecasting. The five and ten-day forecasts 
are based on synoptic (the evolution from past to 
current and future daily and average Grosswetter 
patterns, use of analogs) and statistical (by correla- 
tion and periodic analysis) techniques. The longer 
period forecasts are based more on the statistical, 
and less on the synoptic, techniques. Baur does not, 
in this discussion, elaborate on the details of his 
forecasting techniques more than he has in the past, 
nor does he furnish any systematic verification evi- 
dence as to the success or skill of the forecasts. Con- 
sequently there is little point in reviewing this sec- 
tion further. From the impressions which the re- 
viewer has obtained of Baur’s forecasts and fore- 
casting techniques, including the results of an earlier 
attempt at an objective statistical verification of 
two years of his five- and ten-day forecasts*, there 
is no evidence that his results are significantly better 
than those obtained by similar forecasts prepared by 
other agencies. 

In the appendix, Baur summarizes in one brief 
three-point statement the essential conclusions which 
he draws from the material treated in this volume, 
as follows: “Es existiert ein Grosswetter. Dieses 
wird von der Sonne gelenkt, von dem vorausge- 
gangenen Gross- und Weltwetter geformt.” (There 


2 See Supplement No. 39, Monthly Weather Review, 
pp. 88—91. 
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exists a Gross-wetter. This is controlled by the sun, 
moulded by the preceding Gross- and world wea- 
ther.) This reviewer is in complete agreement with 
this statement, and in general with the author’s 
philosophical approach to the forecast problem, 
but he does not feel that other than by some stati- 
stical refinement the author has contributed much 
that is essentially new or significant to the proof 
of these hypotheses. Particularly notable is the lack 
of any integrating principle or unifying approach 
in his synoptic expression of the state or change of 
state of the general circulation pattern (Gross- 
Weltwetterlage). 

On the other hand, there is no question but that 


R. SCHERHAG: Neue Methoden der Wetter- 
analyse und Wetterprognose. Berlin-Göttingen- 
Heidelberg, 1948. 424 pp., 213 figures. Springer- 
Verlag. 


During the 1930's, the publication “Täglicher 
Wetterbericht der Deutschen Seewarte” was an 
important source of information for those who were 
interested in current weather and current aerological 
observations. It was more than that: it actually 
gave, in figures and text, day after day an instructive 
picture of the working methods employed by the 
meteorological service of the Seewarte. The com- 
prehensive use of aerological material gave this 
weather service the character of a pioneer team 
and the descriptive text of the bulletins often wit- 
nessed of the authors’ deep personal interest in the 
problems. The present writer recalls with special 
pleasure the comments by Dr Scherhag which gave 
much valuable information, pointing out significant 
details which easily enough might have been over- 
looked, apostrophizing extreme conditions in a 
language which called things by their right name, 
and, above all, introducing the reader into the 
workshop of those German meteorologists who 
tried to find new laws and rules by studying three- 
dimensional weather. 

Much has happened since these days. For a few 
years the second world war brought an immense 
increase of the aerological material in Germany 
and countries under German control, and an equally 
great increase of the need for detailed and reliable 
forecasts of surface weather and upper air condi- 


Baur’s efforts in the field of longer-range forecasting, 
which form the basis of this volume, constitute 
pioneer work in this field. The difficulties which 
the investigator has to face in this field are evident 
from the present unsatisfactory status of long-range 
forecasting as a science. Baur certainly deserves 
credit for making the effort to present a very subjec- 
tive and technically esoteric yet most important 
branch of meteorology in. a form which is compre- 
hensible to the average intelligent reader. 


H. C. Willett 


Massachussetts Institute of Technology 


tions. We do not know with what feelings the 
German meteorologists followed their orders in 
pursuing their work under these conditions, so 
favourable and so desperate at the same time, but 
we may suppose that many of them, with the in- 
ternational view which is natural to the scientists of 
a truly international science, received some inspira- 
tion from the conviction that their work should 
prove valuable in a future of peaceful, world-wide 
cooperation. If this assumption is correct, it is 
possible to understand that German meteorologists 
who were helped to continue their scientific work 
after the war ended, have been able to produce very 
valuable articles and books in spite of the hard 
times they have experienced. 

The most comprehensive monography which 
has, until now, resulted from this scientific activ- 
ity, is the work of RicHARD SCHERHAG on New 
Methods of Weather Analysis and Weather Fore- 
casting. This book contains 424 pages in big size, 
and 213 figures most of which are original and of a 
very high class. 

The work consists of five parts entitled as follows: 


(t) Development and Fundaments of Synoptic Me- 
teorology, 

(2) General Circulation of Troposphere and Stra- 
tosphere, 

(3) Weather and Weather Analysis, 

(4) Weather Forecasting, 

(s) Long-Range Forecasting. 


In the first part the author. gives an interesting 
survey of the development of synoptic meteorology. 
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The German schools are, for natural reasons, treated 
more completely than others, and little is said e.g. 
about the great development which has taken place 
within this branch of meteorology in the United 
States during the last decade, but on the whole 
Scherhag’s statements seem to be deliberate and as 
objective as one might claim in view of the fact 
that the author has a personal opinion and is ready 
to stand for it. 

In his treatment of the development of synoptic 
aerology, Scherhag gives good arguments for the 
choice of standard isobaric surfaces made by German 
meteorologists, 1.e. the 700, 500, 225, 96 and 41 mb 
levels. It certainly ought to be discussed whether 
this set of standard levels should not be adopted by 
an international convention. 

“The fundaments of synoptic meteorology” is 
the heading of a rather short chapter (13 pp.) which 
contains the greater part of the very few and ele- 
mentary mathematical equations. of Scherhag’s 
work. This brief treatment, of course, adds little to 
the generally accepted dynamical theories but may 
be justified as a recapitulation of theoretical ideas 
which are, in the author’s opinion, essential for the 
synoptician. A theoretical meteorologist would 
probably, after reading this section, feel discouraged 
once more by the contrast between the magnificent 
building of theory itself and the modest, almost 
crippled part of it which is thought to be applicable 
in the weather service. — Although the treatment 
of the problem is very elementary, a few statements 
and formulations are open to criticism. It is too 
plain language to say, for instance, that the coriolis 
acceleration, being the only force acting on a 
straight-moving air mass, must be compensated by 
the gradient force if the wind should not change 
its direction (p. 21); and it is not convincing when 
one of the fundamental findings ofsynoptic aerology: 
high pressure over warm air, low pressure over 
cold air in the upper troposphere, is said to be caused 
by ground friction (p. 29). 

Part 2, dealing with the general circulation, is of 
much greater value. The many tables of new 
aerological mean values (mainly concerning locali- 
ties in Europe), and the extremely instructive 
charts based on these and other (partly unpublished) 
new means, represent no doubt years of hard work, 
especially hard because so many factors act together 
to impede a statistical treatment of aerological 
observations. Scherhag and others who have shared 
the task with him seem to have spared no effort to 
secure the best possible result for the greatest possible 
area in space and time. Nevertheless, the charts must, 


of course, for vast areas be considered merely as an 
indication of the probable general picture, the 
revision of which may be possible and necessary in 
a very short time. In fact, the German weather 
service in Bad Kissingen is already publishing new 
charts showing seasonal mean values which are, 
for some areas, based on a material that was. not 
accessible for the construction of Scherhag’s maps. 

The most outstanding new information that can 
be gathered from these latter maps is perhaps the 
reversal of the pressure gradient which, during 
summer, is taking place in latitudes north of 30° N 
somewhere between the levels of 96 and 41 mbs. 
The strong easterly current at the 300 mb level over 
India in summer is shown to be the lowest-dipping 
part of a tropical stratospherical ringshaped stream 
which during that time of the year winds around 
the earth in latitudes between 10 and 25° N.! This 
is an interesting counterpart to the polar stream 
which in winter does not reach its maximum at the 
tropopause, as generally believed, but, at latitudes 
north of 50°, continues to increase — on the mean 
charts, at least — up to the highest levels from which 
information could be obtained by the aid of radio- 
sondes. 

The last chapter of part 2 gives a survey of the 
so-called stratospheric compensation, that is to say 
the well-known rule that the stratosphere is cold, 
when the troposphere is warm, and vice versa. The 
treatment which Scherhag gives of this question is 
rather comprehensive — 24 p.— and very interesting. 
It is pointed out that although the rule quoted above 
is confirmed with few exceptions by the material 
for any single station and month, it is not valid 
when the annual variation is considered, the strato- 
sphere in middle and, to a much higher degree, in 
high latitudes being considerably warmer in sum- 
mer than in winter. A detailed study of stratospheric 
compensation has revealed that, besides, in winter 
there are variations in the degree of compensation 
from day to day which must be connected somehow 
with different circulation types; the main periods 
seem to be of the order of magnitude of 20 days or 
more. Some outstanding cases of deviations from 
the normal compensation are given, with an inter- 
esting account of relationship to tropospherical 
conditions and the geographical distribution of the 
said anomaly. The idea of taking deviations from 


1 Scherhag suggests that this equatorial ring stream 
prevents any appreciable outflow (“einen stärkeren Ab- 
fluss’’) from the summer to the winter hemisphere, 
but gives no reasons for his suggestion. 
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an empirically determined “normal” compensation 
(which is a function of the season and the geographi- 
cal coordinates) is favourable in several respects; the 
applicability is, of course, limited to such areas and 
seasons where compensation is the normal thing, 
but in middle latitudes this condition is fulfilled. In 
addition to his empirical results, Scherhag gives a 
short outline of the theory concerned. Finally, the 
behaviour of the wind in the compensation layer is 
illustrated by figures which are based on a rather 
great number of radar wind measurements at Berlin 
and show the normal changes of wind velocity 
and wind direction between 11,000 and 16,000 m 
as a function of wind velocity and wind direction, 
respectively, at 11,000 m. It ought to be stressed, 
however, that each of the dots in these figures 
(55—58, p. 106—107) does not represent an indi- 
vidual case but rather a group of cases. One addi- 
tional diagram, at least, showing the scattering of 
observed individual values within a month or a 
season, would have been of considerable interest. 

The third part, which deals with weather and 
weather analysis, is the bulk of the whole book, 
containing nearly 200 pages. It consists of three 
chapters, entitled Elements of the weather map, 
Analysis of the weather map, and Typical weather 
situations (Wetterlagen) in Europe. The disposition 
of the great variety of facts and thoughts given in 
the first of these chapters may at the first glance 
seem puzzling, but a closer study will make the 
reader understand and, quite probably, accept the 
reasons for the untraditional order of presentation. 
The dynamical anticyclone is treated at first, fol- 
lowed by the cold anticyclone, the air-masses, the 
frontal zone, the cyclones and their fronts, the 
hurricanes, high-tropospherical processes, the cold 
pool (“Kaltlufttropfen”) and, at last, stratospheric 
weather processes. 

The general description of air-masses, fronts, 
cyclones, and anticyclones is, of course, on the 
whole in good agreement with the treatment given 
in other text-books of later years. Nevertheless, it 
can be gathered from almost every page how deeply 
the author is personally engaged in the problems 
which he is discussing. Scherhag is equally far from 
a dogmatism which is inclined to reproach Nature 
with its poor fitting to theory, and a nihilism which 
finds everything diffuse and turbid and nebulous. 

For some practical purposes, one might have 
wished that the author should have chosen his 
examples, or at least a great part of them, from 
weather phenomena of moderate intensity, from 
the everyday problems of analysis and forecasting. 


No doubt, a collection of such examples might have 
been useful, too, for students of meteorology, but 
it is not difficult to understand that Scherhag has 
preferred to select situations of special interest, thus 
making his book something of a catalogue of the 
most outstanding weather events in Europe and 
adjacent parts of the Atlantic during more than 
so years. This should not be considered as a kind 
of journalistic record-hunting; in fact, the circum- 
stances under which meteorological extremes occur 
deserve serious consideration and are well worth a 
description, provided this description includes an 
account of the physics of the case that is treated. 

Some of the definitions and statements in this 
section are open to criticism. A front is (p. 248) 
defined as a narrow boundary zone, where air 
masses of different origin and different character 
are moved towards one another; it seems doubtful 
whether this definition accounts for e.g. a fast 
moving cold front. — The rules given by GUILBERT 
in his “Nouvelle Methode de prévision du temps” 
are treated with a remarkable esteem; one might 
suppose that Scherhag feels indebted to Guilbert 
whose work may with some reservation be said to an- 
ticipate the modern prognostic chart and Scherhag’s 
own divergence theory. In the opinion of the pres- 
ent writer, however, Guilbert’s often quoted rule 
of succession of pressure ridges and troughs with 
24 hours interval is in lack of every theoretical 
support and of little or no practical value. — Re- 
ferring to v. FICKER’S so called primary and secondary 
waves, Scherhag writes that a phase difference of 
1/4 wave length between the ground level and the 
upper troposphere can be brought about by a 
horizontal temperature gradient in the layer be- 
tween. This, of course, is true, but if so, the lowest 
mean temperature of the troposphere must be 
found more than !/4 wave length behind the pressure 
minimum at the ground, or, in other words, the 
highs must on an average be colder than the lows, 
and this is, in fact, an exception rather than the rule. 
The phase difference between high-tropospheric and 
low-tropospheric pressure waves is less than the 
phase difference between temperature and low- 
tropospheric pressure waves, and even the latter is 
in most cases less than 1/4 wave length. 

It shall not be concealed that besides these points 
there are some other passages in Scherhag’s book 
where one may have reasons to put a question mark 
on the margin. The book is, on the whole, written 
in such a way that the reader has no chance to fall 
asleep, and some of the weak points may, indirectly, 
be useful as incitements to new research work; with 
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a few exceptions, the weak points are not simply 
errors. 

The main value of Scherhag’s work lies, no doubt, 
in the straightforward and skilful analysis of a great 
number of charts, including upper air charts and 
mean charts (but, curiously enough, no cross- 
section of the type which is so frequent in aerological 
literature), and the accompanying text with its 
multitude of important remarks on general pro- 
blems. Nearly all situations described are of great 
interest, either for illustrating well-known facts or 
for the discussion of new or intricate problems. 
Some of the examples may be briefly mentioned 
here: the great anticyclone over the arctic in the 
middle of Jan. 1940 (p. 124), the famous anticyclone 
over Europe about Jan. 23, 1907 (p. 134), the in- 
structive case of cyclogenesis in a frontal zone over 
eastern North America about Febr. 8, 1933 (p. 
157), the extremely deep cyclones of Dec. 8, 1886 
(p. 179) and Jan. 2, 1933 (p. 187), the unprecedented 
outbreak of Siberean arctic air to Europe about 
Jan. 25, 1942 (p. 233), and the severe gale over the 
Baltic region on July 9, 1931 (p. 280). 

The fourth part, called Weather Forecasting, 
attaches to the third part by continuing the descrip- 
tion of weather situations; the situations described 
in this section are all of them complicated and diffi- 
cult from a prognostic point of view, except one 
(p. 326) which has been inserted to illustrate the 
construction method of prognostic charts. It is 
true that this method stands out quite clearly just 
because of the uncomplicated situation to which 
it is applied but it might have been more natural to 
start with some examples of this type and proceed 
to more intricate prognostic problems. It should 
have been of considerable interest to see a detailed 
description of the construction of prognostic charts 
in a few really difficult cases: one should exemplify 
the striking successes which are sometimes achieved; 
one or two should demonstrate the equally striking 
failures which remind us of the imperfection of 
this (and, probably, any other) approach to a solu- 
tion of the problem of forecasting weather. 

As in the whole book, the casuistic treatment of 
the subject is much fuller and much more valuable 
than the deductive passages. It must be acknow- 
ledged, though, that the author has made an effort 
to compile some rules which, presumably, are of 
prognostic value; these rules, most of which are 
formulated by Scherhag or RopEWALD, no doubt 
deserve the attention of forecasters, but no final 


judgment as to their practical value or theoretical . 


foundation can be given as yet. — In a section con- 
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cerning forecasting of individual weather elements 
the author renders much valuable advice. Some 
principal remarks upon the formulation of forecasts 
are brillantly written and ought to be studied 
thoroughly by all weather forecasters. The warnings 
against “Gummiprognosen” and “a deadly levelling” 
of the wording of a forecast are not only of much 
need but also characteristic for the spirit of the 
author. 

The last few pages of the text deal with the pro- 
blem of long-range forecasting. The method of 
using analogous situations is briefly described, and 
the difficulties in finding sufficiently close analogies 
are discussed. Scherhag’s pessimism on this point, 
as well as concerning the method of multiple corre- 
lations, may be justified, but it stands out in marked 
contrast to the benevolence with which he presents 
such items as symmetry points, reflexion points, 
“singularities”, and sun spot cycles in the atmosphere. 
As far as the so called wave phenomena studied by 
the Leipzic school are concerned, there may be 
some reason to suppose that they are not merely 
accidental in character, but when Dr. Scherhag, 
evidently, believes that long-range forecasts may 
well be based on these theories it is difficult to follow 
him. As for the singularities, Scherhag bases his 
discussion upon the works of ScHMauss, FLOHN 
and other German meteorologists. It is true that 
he warns against an exaggerate use of the singulari- 
ties, and his statement that an increased probability 
(“eine besondere Bereitschaft’’) for a certain weather 
type is present at certain times is not incorrect in 
itself, if the term af certain times is understood to 
mean during periods of the duration of several 
weeks; but one gets the impression that Scherhag 
is ready to accept a good deal of the results which the 
singularity specialists claim to have arrived at. 
After the book was written, professor J. BARTELS 
has given a very clear, sceptical and convincing 
survey of the whole problem from a statistical point 
of view in articles published in “Annalen der 
Meteorologie”, and it may be hoped that, on account 
of this survey, the whole discussion of the singu- 
larities may in a following edition of Scherhag’s book 
be replaced by a mainly historical treatment of the 
singularity epoch in German meteorology. 

The literature list given at the end of the book is 
unusually comprehensive (897 numbers) and seems 
to have been compiled with great care. Among 
the articles mentioned in the list, some have been 
issued in periodicals or series which are not familiar 
to meteorologists in general, a fact which increases 
the value of the list for the research worker. 
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The make-up of the book is of a high standard. 
The figures have been drawn with great care, and 
the reproduction is remarkably fine. But it should 
not be concealed that this high technical standard 
must have been an important cause for the high 
price of Scherhag’s book (75 Reichsmark!). It must 
be hoped that the book will be translated into 
English and French and issued in a number of 
specimens sufficiently great to reduce to price 
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considerably. If such a translation is made, there 
will be an opportunity to correct a few errors, 
some of which have been mentioned above, but 
this will not and should not change the general 
character or the personal note of Scherhag’s great 
work. 

E. Hovmöller 


Swedish Meteorological and 
Hydrological Institute. 
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